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Abstract 


'I  nis  Semiannual  Technical  Summary  describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NOllSAR),  the  Nor¬ 
wegian  Regional  Seismic  Array  (KORKSS)  and  the  Arctic  Regional  Seismic 
Array  (ARCESS).for  the  period  1  October  1990  -  31  March  1991.  Statistics 
are  also  presented  for  additional  seismic  stations,  which  through  cooperative 
agreements  with  institutions  in  the  host  countries  provide  continuous  data 
to  the  NORS.AR  Data  Processing  Center  (NI’DC).  These  stations  comprise 
the  Finnish  K.xperimental  Seismic  Array  (FINKS.A),  the  German  Experimental 
Seismic  Array  (GERESS),  and  two  ,3-component  stations  in  Poland:  Ksiaz  and 
Stary  Folwark.  This  Semiannual  Report  also  presents  statistics  from  operation 
of  the  Intelligent  Monitoring  System  (I.MS).  rhe  l.MS  htus  boon  operated  in  an 
experimental  mode  using  NORESS  and  ARCESS  data,  and  the  performance 
has  been  very  satisfactory. 

The  NORSAR  Detection  Processing  sy:>tem  has  been  operated  through¬ 
out  the  period  with  an  average  uptime  of  9S.(i%  as  compared  to  98.0%  for  the 
previous  reporting  period.  .A  total  of  1837  seismic  events  have  been  reported 
in  the  NORS.AR  monthly  seismic  bulletin.  The  performance  of  the  contin¬ 
uous  alarm  .system  and  the  automatic  bulletin  transfer  by  telex  to  AFTAC 
ha.s  been  satisfactory.  .A  system  for  direct  retrieval  of  NORSAR  waveform 
data  through  an  X.2.3  cotmection  has  been  implemented,  and  has  been  tested 
succe.ssfully  for  actiuiring  such  data  by  .AFT.AC.  J’roressing  of  retjnests  for 
full  NORS.AR/NORESS  data  on  magnetic  tapes  has  progressed  according  to 
established  schedules. 

On-line  detection  proce.ssing  and  data  recording  at  the  NORSAR  Data 
Processing  Center  (.NDPC)  of  NORESS,  ARCESS  and  FINESA  data  have 
been  conducted  throughout  the  period,  with  an  average  uptime  of  99.2%  for 
NORESS,  98.‘1%  for  ARCESS  and  98.9%  for  FINESA.  The  Intelligent  .Mon¬ 
itoring  System  was  installed  at  NORSAR  in  December  1989  and  has  been 
operated  experimentally  since  1  January  1990.  Results  of  the  IMS  analysis  for 
the  reporting  period  are  given. 

There  have  been  no  modilications  made  to  the  NORSAR  data  acquisition 
system.  The  procc.s.s  of  evaluating  and  testing  technical  options  for  refurbish¬ 
ment  of  the  array  is  continuing. 

The  routine  detection  processing  of  NORESS,  ARCESS  and  FINESA  is 
running  satisfactorily  on  each  of  the  arrays’  SUN-3/280  data  acquisition  sys- 
tem.s.  The  routine  proce.s.sing  of  FI.NES.A  data  at  NORSAR  is  similar  to  what 
is  done  in  Helsinki.  GERESS  data  acquisition  and  detection  processing  has 
been  conducted  in  an  experimental  mode  during  the  period,  in  a  way  similar 
to  what  is  done  in  Bochum. 


Maintenance  activities  in  the  period  comprise  preventive/corrective  main¬ 
tenance  in  connection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS. 
In  addition,  the  maintenance  center  has  been  involved  with  occasional  main¬ 
tenance  of  equipment  for  FINESA  and  preparatory  work  in  connection  with 
the  two  stations  in  Poland.  Other  activities  involved  testing  of  the  NORSAR 
communications  systems. 

We  have  continued  our  work  aimed  at  evaluating  the  stability  of  RMS  Lg 
for  yield  estimation  purposes.  We  have  carried  out  a  detailed  RMS  Lg  anal¬ 
ysis  of  NORSAR  recordings  of  Novaya  Zemlya  underground  explosions,  and 
in  addition  conducted  similar  analysis  of  Grafenberg  (GRF)  array  recordings 
(available  after  1976).  The  results  show  a  tight  clustering  in  M(Lg)  of  13 
NORSAR-recorded  explosions  after  1976,  and  GRF  data  of  the  27  Sep  78  ex¬ 
plosion,  for  which  no  NORSAR  data  are  available,  indicate  that  this  event  is 
of  similar  size.  The  correspondence  between  M(Lg)  for  NORSAR  vs  GRF  is 
excellent,  with  an  orthogonal  standard  deviation  of  only  0.02-0.03  magnitude 
units.  This  correspondence  fully  matches  the  excellent  results  previously  ob¬ 
tained  for  Semipalatinsk  explosions,  and  confirms  the  promise  of  RMS  Lg  as 
a  stable  estimator  of  relative  yields. 

The  concept  of  threshold  monitoring,  introduced  by  Ringdal  and  Kvaerna, 
is  a  method  of  monitoring  the  seismic  amplitude  levels  for  the  purpose  of  using 
this  information  to  assess  the  largest  size  of  events  that  might  go  undetected  by 
a  given  network.  In  an  effort  to  demonstrate  the  capabilities  of  this  threshold 
monitoring  technique,  we  have  conducted  a  simulation  experiment,  which  has 
involved  down-scaling  the  recorded  signals  of  the  24  October  1990  Novaya 
Zemlya  explosion  (mi,  =  5,6)  by  a  factor  of  1000  (i.e.,  3  orders  of  magnitude). 
The  resulting  NORESS,  ARCESS  and  FINESA  traces  of  this  “ms  =  2.6  event” 
were  added  at  hourly  intervals  to  the  actual  recordings  for  a  full  data  day,  and 
the  threshold  monitor  was  then  applied.  The  results  demonstrated  that  every 
one  of  these  24  “events”  were  clearly  identifiable  on  the  threshold  trace.  While 
this  clearly  gives  a  very  simplified  picture,  it  serves  to  document  the  excellent 
monitoring  potential  of  these  three  arrays  for  the  Novaya  Zemlya  test  site. 

A  significant  part  of  our  research  effort  has  been  directed  toward  the  fur¬ 
ther  development  of  European  high  frequency  arrays  and  high  quality  three- 
component  stations.  The  aim  has  been  to  provide  for  continuous  transmission 
of  data  (by  satellite  or  land  line)  to  NDPC  and  integration  of  this  data  stream 
into  the  input  base  for  the  Intelligent  Monitoring  System  (IMS).  These  efforts 
have  proceeded  satisfactorily.  During  the  reporting  period,  particular  em¬ 
phasis  has  been  placed  upon  optimizing  the  GERESS  beam  deployment  (in 
cooperation  with  Ruhr  University  scientists)  and  completing  the  integration 
into  the  network  of  the  KSP  and  SFP  stations  in  Poland.  All  of  these  network 
stations  are  scheduled  for  participation  in  the  1991  GSETT-2  experiment. 
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The  multichannel  statistical  data  processing  algorithms  described  in  previ¬ 
ous  NORSAR  Semiannual  Technical  summaries  have  now  been  integrated  into 
the  Event  Processor  system  currently  in  operation  at  NORSAR.  A  description 
of  the  program  and  some  of  the  recent  results  are  provided.  At  pre.sent.  these 
algorithms  are  used  only  experimentally,  in  an  offline  fashion,  but  plans  to 
start  testing  them  in  an  operationcd  environment  have  been  made. 

Results  are  presented  from  a  two-dimensional  finite  difference  approach  to 
modeling  seismic  wave  propagation  in  the  crust.  These  results  are  based  on  a 
cooperative  effort  with  scientists  at  the  HIM  liergen  Scientific  Ontre.  In  the 
generation  of  synthetic  seismogram.s,  a  homogeneous  crust  of  thickness  ,10  km 
and  =  6.5  km/sec  has  been  used  as  a  basis,  with  options  for  pertiirba 
tion  comprising  multilayering,  piece-wise  linear  velocity  gradients  as  well  as 
large-scale  discontinuities.  Randomized  scattering  effects  have  so  far  not  I.een 
considered.  It  appears  that  the  synthetics  generate  all  major  phases,  but  with 
a  relatively  weak  body-wave  coda  generation  compared  to  real  seismograms. 
The  introduction  of  scatterers  is  expected  to  be  of  importance  in  remedying 
this  problem,  and  will  be  the  subject  of  future  investigations. 

An  overview  is  presented  of  results  related  to  determining  the  crustal  thick¬ 
nesses  in  Fennoscandia.  A  detailed  contour  map  (2  km  contour  intervals)  has 
been  developed  for  the  entire  region.  In  view  of  the  extrensive  sediment  de¬ 
posits,  a  map  of  the  crystalline  crustal  thicknesses  is  ahso  presented  for  the  SVV 
part  of  the  area.  In  general,  the  oldest  parts  of  the  Baltic  Shield  exhibit  the 
greatest  crustal  thicknesses  (in  somf-  areas  exceeding  50  km).  In  the  offshore 
Norway  areas,  the  crystalline  thicknesses  are  of  the  order  of  15-20  km.  while 
the  sedimentary  overburdens  can  exceed  10  km.  An  interesting  observation  is 
that  the  Moho  depth  variation  appears  to  have  a  counterpart  in  the  spatial 
distribution  of  earthquakes  in  the  Fennoscandian  region. 

In  the  current  I.MS  implementation  of  the  threshold  monitoring  methe',  a 
limited  number  of  specific  target  sites  are  monitored.  These  sites  include  sev¬ 
eral  mines  in  Scandinavia  and  Western  Russia,  along  with  the  Novav  "emlya 
and  Semipalatinsk  nuclear  test  sites.  We  have  now  initiated  a  s‘  ,  ',y  to  de¬ 
termine  how  this  method  could  be  applied  to  monitoring  more  c  <lensive  geo¬ 
graphical  regions.  The  key  to  achieving  this  is  to  develop  “generic”  relations 
for  attenuation  and  magnitude  correction  factors  for  seismic  phases  of  interest, 
and  to  deploy  a  sufficient  number  of  beams  to  ensure  adequate  coverage.  So 
far,  we  have  developed  preliminary  relations  and  corre,,i,on  factors  for  the  Pn 
and  Lg  phases.  These  “generic”  relations  are  based  upon  systematic  analysis 
of  several  hundred  phase  observations  of  regiona'  events  in  various  geographi¬ 
cal  areas.  The  results  are  applicable  to  Northe:!'  Europe  and  adjacent  regions. 
Factors  that  remain  to  be  assessed  in  detai'  .nclude  the  effect  of  uncertainties 
in  reference  M/,  magnitudes  for  the  evens  in  the  data  base  and  the  effect  of 
signal  loss  in  the  array  beamforming. 
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Using  the  “generic”  amplitude  relations  described  above,  we  have  gen¬ 
erated  a  very  extensive  beam  set  for  a  short  test  interval  using  NORESS, 
ARCESS  and  FINESA  data.  The  results  of  TM  processing  for  this  time  pe¬ 
riod,  which  are  presented  in  the  form  of  regional  threshold  maps,  are  very 
encouraging,  and  indicate  that  this  is  indeed  a  potentially  very  useful  exten¬ 
sion  of  the  concept.  However,  an  “operational”  implementation  might  require 
computer  processing  capacity  of  about  an  order  of  magnitude  greater  than  is 
currently  available  on  a  Sun  Workstation  (type  Sparc  station  1). 
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Summary 


This  Semiannual  Technical  Summary  describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Nor¬ 
wegian  Regional  Seismic  Array  (NORESS)  and  the  Arctic  Regional  Seismic 
Array  (ARCESS)  for  the  period  1  October  1990  -  31  March  1991.  Statistics 
are  also  presented  for  additional  seismic  stations,  which  through  cooperative 
agreements  with  institutions  in  tfie  host  countries  provide  continuous  data 
to  the  NORSAR  Data  Processing  Center  (NPDC).  These  stations  comprise 
the  Finnish  Experimental  Seismic  Array  (FINESA),  the  German  Experimental 
Seismic  Array  (GERESS),  and  two  3  component  stations  in  Poland:  Ksiaz  and 
Stary  Folwark.  This  Semiannual  Report  also  presents  statistics  from  operation 
of  the  Intelligent  Monitoring  System  (IMS).  The  IMS  has  been  operated  in  an 
experimental  mode  using  NORECS  and  .ARCESS  data,  and  the  performance 
has  been  very  satisfactory. 

The  NORSAR  Detection  Proce.ssing  system  has  been  operated  through¬ 
out  the  period  with  an  average  uptime  of  98.6%  as  compared  to  98. 0%,  for  the 
previous  reporting  period.  .A  total  of  1837  seismic  events  have  been  reported 
in  the  NORSAR  monthly  .seismic  bulletin.  The  performance  of  the  contin¬ 
uous  alarm  system  and  the  automatic  bulletin  transfer  by  telex  to  .AFTAC 
has  been  satisfactory.  A  system  for  direct  retrieval  of  NORSAR  waveform 
data  through  an  .X.2')  connection  has  been  implemented,  and  has  been  tested 
successfully  for  acquiring  such  data  by  .AFT.AC.  Processing  of  requests  for 
full  NORSAR/NORESS  data  on  magnetic  tapes  has  progressed  according  to 
established  .schedules. 

On-line  detection  processing  and  data  recordiitg  at  the  NORS.AR  Data 
Processing  Center  (.NDPC)  of  NORESS.  .ARCESS  and  FINESA  data  have 
been  conducted  throughout  the  period,  with  an  average  uptime  of  99.2%i  for 
NORESS,  98.4%  for  ARCE.SS  and  98,9%  for  FINESA.  The  Intelligent  Mon¬ 
itoring  System  was  installed  at  NORS.AR  in  December  1989  and  has  been 
operated  experimentally  since  1  .lanuary  1990.  Results  of  the  IMS  analysis  for 
the  reporting  period  are  given. 

There  have  been  no  modifications  made  to  (fie  NORS.AR  data  acqu  sition 
system.  The  process  of  evaluating  and  testing  technical  options  for  refurbish¬ 
ment  of  the  array  is  continuing. 

The  routine  detection  processing  of  NORESS,  ARCESS  and  FINESA  is 
running  satisfactorily  on  each  of  the  arrays’  StIN-3/280  data  acquisition  sys- 
tem.s.  The  routine  proces.sing  of  FINESA  data  at  NORSAR  is  similar  to  what 
is  done  in  Helsinki.  GERESS  data  acquisition  and  detection  processing  has 
been  conducted  in  an  experimental  mode  during  the  period,  in  a  way  similar 
to  what  is  done  in  Bochum. 
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Maintenance  activities  in  the  period  comprise  preventive/corrective  main¬ 
tenance  in  connection  with  all  the  N’ORSAR  snharrays,  NORESS  and  ARCESS. 
In  addition,  the  maintenance  center  has  been  involved  with  occasional  main¬ 
tenance  of  equipment  for  I'lNESA  and  preparatory  work  in  connection  with 
the  two  stations  in  Poland.  Other  activities  involved  testing  of  the  XORSAR 
communications  systems. 

We  have  continued  our  work  aimed  at  evaluating  the  stability  of  RMS  Lg 
for  yield  estimation  purposes.  We  have  carried  out  a  detailed  RMS  Lg  anal¬ 
ysis  of  NORS.\R  recordings  of  Novaya  Zemlya  underground  explosions,  and 
in  addition  conducted  similar  analysis  of  Grafenberg  (GRF)  array  recordings 
(available  after  1976).  The  results  show  a  tight  clustering  in  M(Lg)  of  13 
NORSAR-recorded  explosions  after  1976,  and  GRF  data  of  the  27  Sep  78  ex¬ 
plosion,  for  which  no  NORS.\R  data  are  available,  indicate  that  this  event  is 
of  similar  size.  The  correspondence  between  .M(Lg)  for  NORSAR  v.s  GRF  is 
excellent,  with  an  orthogonal  standard  deviation  of  only  0.02-0.03  magnitude 
units.  This  correspondence  fully  matches  the  excellent  results  previously  ob¬ 
tained  for  Semipalatinsk  explosions,  and  confirms  the  promise  of  RMS  Lg  as 
a  stable  estimator  of  relative  yields. 

The  concept  of  threshold  monitoring,  introduced  by  Ringdal  and  Kvmrna, 
is  a  method  of  monitoring  the  seismic  amplitude  levehs  for  the  purpo.se  of  u.sing 
this  information  to  assess  the  largest  size  of  events  that  might  go  undetected  by 
a  given  network.  In  an  effort  to  demonstrate  the  capabilities  of  this  threshold 
monitoring  technique,  we  have  conducted  a  simulation  experiment,  which  has 
involved  down-.scaling  the  recorded  signals  of  the  24  October  1990  Novaya 
Zemlya  explosion  {m;,  =  .5.6)  by  a  factor  of  1000  (i.e.,  3  orders  of  magnitude). 
The  resulting  NORESS,  ARCESS  and  FINES  .A  traces  of  this  “ms  =  2.6  event" 
were  added  at  hourly  intervals  to  the  actual  recordings  for  a  full  data  day,  and 
the  threshold  monitor  was  then  applied.  The  results  demonstrated  that  every 
one  of  these  24  “events”  were  clearly  identifiable  on  the  threshold  trace.  While 
this  clearly  gives  a  very  simplified  picture,  it  serves  to  document  the  excellent 
monitoring  potential  of  these  three  arrays  for  the  Novaya  Zemlya  test  site. 

A  .significant  part  of  our  research  effort  has  been  directed  toward  the  fur¬ 
ther  development  of  European  high  frequency  arrays  and  high  quality  three- 
component  stations.  The  aim  has  been  to  provide  for  continuous  transmission 
of  data  (by  satellite  or  land  line)  to  NDPC  and  integration  of  this  data  stream 
into  the  input  base  for  the  Intelligent  Monitoring  System  (IMS).  These  efforts 
have  proceeded  satisfactorily.  During  the  reporting  period,  particular  em¬ 
phasis  has  been  placed  upon  optimizing  the  GERESS  beam  deployment  (in 
cooperation  with  Ruhr  University  scientists)  and  completing  the  integration 
into  the  network  of  the  KSP  and  SFP  stations  in  Poland.  All  of  these  network 
stations  are  scheduled  for  participation  in  the  1991  GSETT-2  experiment. 
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The  multichannel  statistical  data  processing  algorithms  described  in  previ¬ 
ous  NORSAR  Semiannual  Technical  summaries  have  now  been  integrated  into 
the  Event  Processor  system  currently  in  operation  at  NORSAR.  A  description 
of  the  program  and  some  of  the  recent  results  are  provided.  At  ore  -'ut,  these 
algorithms  are  used  only  experimentally,  in  an  offline  fashion,  o”!  plans  to 
start  testing  them  in  an  operational  environment  have  been  made. 

Results  are  presented  from  a  two-di..  nal  finite  difference  approach  to 
modeling  seismic  wave  propagation  in  the  ci.  .  These  results  are  based  on  a 
cooperative  effort  with  scientists  at  the  IBM  Bergen  Scientific  Centre.  In  the 
generation  of  synthetic  seismograms,  a  homogeneous  crust  of  thickness  .10  km 
and  =  6..>  km/sec  ha.s  been  used  as  a  ba.sis,  with  options  for  perturba¬ 
tion  comprising  multilayoring,  piece-wise  linear  velocity  gradients  a-s  well  as 
large-scale  discontinuities.  Ilandomized  scattering  effects  have  so  far  not  been 
considered.  It  appears  that  the  synthetics  generate  all  major  phases,  but  with 
a  relatively  weak  body-wave  coda  generation  compared  to  real  seismograms. 
The  introduction  of  scatterers  is  expected  to  be  of  importance  in  remedying 
this  problem,  and  will  be  the  subject  of  future  investigations. 

.\n  overview  is  presented  of  results  relateil  to  determining  the  crustal  thick¬ 
nesses  in  f'enno.scandia.  A  detailed  contour  map  (2  km  contour  intervals)  has 
lieen  developed  for  the  entire  .egion.  In  view  of  the  extren.sive  sediitient  de- 
posit.s.  a  map  of  the  crystalline  crustal  thicknesses  is  also  presented  for  the  SVV 
part  of  the  area.  In  gein'ral,  the  olde.st  parts  of  the  Baltic  Shield  exhibit  the 
greatest  crustal  thicknesses  (in  some  areas  exceeding  .bO  km).  In  the  off.-;iiore 
Norway  areas,  the  crystalline  thicknesses  are  of  the  order  of  N't  20  km.  wiiile 
tlie  sedimentary  overburdens  can  exceed  10  km.  An  interesting  observation  is 
that  the  Moho  depth  variation  appears  to  have  a  counterpart  in  the  spatial 
distribution  of  earthquakes  in  the  Fentioscaiidian  region. 

In  the  current  IMS  implementation  of  the  threshold  monitoring  method,  a 
limited  number  of  specific  target  sites  are  monitored.  These  sites  include  sev¬ 
en  al  mines  in  Scandinavia  and  Western  Russia,  along  with  the  Novaya  Zemlya 
and  .Semipalatin.sk  nuclear  te.st  site.s.  We  have  now  initiated  a  study  to  de¬ 
termine  how  this  method  could  be  applied  to  monitoring  more  extensive  geo¬ 
graphical  regions.  The  key  to  achieving  this  is  to  develop  “generic”  relations 
for  attenuation  and  magnitude  correction  factors  for  seismic  phases  of  interest, 
and  to  deploy  a  sufficient  number  of  beams  to  ensure  adequate  coverage.  So 
far.  we  have  developed  preliminary  relations  and  correction  factors  for  the  Pn 
and  I.g  phases.  The.se  “generic”  relations  are  based  upon  systematic  analysis 
of  several  hundred  phase  observations  of  regional  events  in  various  geographi- 
•  al  areas.  The  results  are  applicable  to  Northern  Europe  and  adjacent  regions. 
Factors  that  remain  to  be  assessed  in  detail  include  the  effect  of  uncertainties 
in  reference  Mr,  magnitudes  for  the  events  in  the  data  base  and  the  effect  of 
signal  loss  in  the  array  beamforming. 
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Using  the  '‘generic”  amplitude  relations  described  above,  we  have  gen¬ 
erated  a  very  extensive  beam  set  for  a  short  test  interval  using  NORESS. 
ARCESS  and  FINESA  data.  The  results  of  T\1  processing  for  this  time  pe¬ 
riod,  which  are  presented  in  the  form  of  regional  threshold  maps,  are  very 
eiicouragiug.  and  indicate  that  this  is  indeed  a  polentially  very  useful  exten¬ 
sion  of  the  concept.  However,  an  “operational"  implementation  might  require 
coinpiitor  processing  capacity  of  about  an  order  of  magnitude  greater  than  is 
currently  available  on  a  Siin  Workstation  (type  Sparc  station  1). 
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2  NORSAR  Operation 


2.1  Detection  Processor  (DP)  operation 

There  have  hcen  48  breaks  in  tlie  otherwise  conlimious  operation  of  the  .\OR- 
S.\R  online  system  within  the  eurrenl  (i-inonth  reporting  interval.  'The  uptime 
liereentage  for  the  period  is  08.G  as  compared  to  OS.O  for  the  previous  period. 

Fig.  2.1.1  and  the  accompanying  Table  2.1.1  both  show  the  daily  DR  down¬ 
time  for  the  days  between  1  October  1090  and  31  March  1091.  The  monthly 
recording  times  and  percentages  are  given  in  Table  2.1.2. 

The  breaks  can  be  griruped  as  follows; 


a)  Hardware  failure  8 

b)  .Stops  related  to  program  work  or  error  2 

c)  Hardware  maintenance  slops  G 

d)  Power  jumps  and  breaks  0 

e)  TOD  error  correction  0 

f)  Comtminication  lines  32 


The  total  downtime  for  the  period  was  86  hours  and  3  minutes.  The 
mean-titne-between-failures  (MTIIF)  was  4.1  days,  as  compared  to  3.2  for  the 
previous  period. 

J.  Torstveit 
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Fig,  2.1.1  Detection  Processor  downtime  for  October  (top),  November  (mid¬ 
dle)  and  December  (bottom)  l<)t)0. 
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Fig.  2.1.1  Detection  Processor  downtime  for  January  (top),  February  (mid¬ 
dle)  and  March  (bottom)  1991. 
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Date 

Time 

Cause 

S 

Oct 

0032-0146 

Hardware 

failure  NDPC 

6 

Qct 

1656-1758 

Hardware 

failure  NDPC 

9 

Qct 

1216-1252 

Hardware 

failure  NDPC 

10 

Oct 

1304-1314 

Hardware 

failure  NDPC 

11 

Oct 

0042-0144 

Hardware 

failure  NDPC 

11 

Oct 

0754-1006 

Hardware 

maintenance  NDPC 

20 

Oct 

1624- 

Hardware 

failure  NDPC 

21 

Oct 

-1043 

Hardware 

failure  NDPC 

22 

Qct 

2300- 

Hardware 

failure  NDPC 

23 

Oct 

-1245 

Hardware 

failure  NDPC 

16 

Nov 

0855-0903 

Hardware 

maintenance  NDPC 

21 

Nov 

1801-1850 

Hardware 

failure  NDPC 

28 

Nov 

0120-0140 

Operator 

error  NDPC 

29 

Nov 

1515-1534 

Hardware 

failure  NDPC 

4 

Dec 

1212-1223 

Hardware 

maintenance  NDPC 

1 

Jan 

1037-1351 

Software 

work  (new  year) 

3 

Jan 

1017-1023 

Hardware 

maintenance  NDPC 

21 

Feb 

1620- 

Hardware 

failure  NDPC 

22 

Feb 

-0855 

Hardware 

failure  NDPC 

Table  2.1.1  The  major  downtimes  in  the  period  1  October  1990  31  March 

1991. 
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Month 

DP  Uptime 
Hours 

DP  Uptime 
% 

No.  of 

DP  Ur<‘aks 

No.  of  Days 
with  Breaks 

DP  MTBF* 

(days) 

OCT  90 

70.'1.52 

94.83 

12 

10 

2.3 

NOV  90 

717.58 

99.61 

11 

8 

2..5 

DEC  90 

743.55 

99.96 

4 

4 

6.2 

JAN  91 

740.33 

99.51 

7 

7 

3.9 

FEB  91 

655.08 

97.48 

n 

7 

2.3 

MAR  91 

719.88 

99.98 

3 

3 

7.5 

98.56 

48 

39 

4.1 

‘Mean-liine-bolween-failuros  =  lolal  iiptiine/iio.  of  up  inlorvals. 


Table  2.1.2  Online  system  performance,  1  October  1990  31  Marcli  1991. 


2.2  Array  communications 


General 

Table  2. '2.1  roflecls  ihe  performaiire  of  the  coniiminiratioiis  system  through¬ 
out  the  reporting  period. 

The  iiio.st  common  events  which  have  alfected  the  NOKSAR  array  have 
been:  damaged  communications  cables,  loss  of  suharray  power,  irregularities 
in  the  NTA  transmission  systems  and  incorrect  line  level. 

Detailed  summary 

OrlolxT  (Wfxks  40  44).  1.10  4.11.90 

()1.\  was  affecleil  by  a  bad  communications  cable  (week  41).  OlH  was 
affected  by  a  communications  cable  damaged  by  excavation  (weeks  40  41  ). 
.\lso  the  ()'.’(■  coinmiinicalions  (able  was  aliected  for  the  same  r('ason  (week 
l(i).  In  iiddilion  we  had  sync  problems  (week  44)  in  connection  with  0211.  02(' 
and  OOC. 

,V()iy  rnht  r  ('ici  f  k.s  ./.)  TV-  J  14.90 

01. 4  resumed  operation  ti  .Vovetnher.  Dili  was  a/fctcted  appro.vimatejy  0 
hours  week  4-4,  can.se  not  staled.  02('  resumed  o|>eralion  20  Nnveinber  after 
comniiinications  cable  reptiir  by  N'T.\/llamar.  OK'  was  affected  "J.  21  and 
jx  .N'oveinlx'r  in  connection  with  an  irregularity  between  llamar  and  Klverum. 
including  a  change  to  spare  transmission  ('(piipment.  OOC  was  affected  weeks 
4.'').  l-'i  and  47.  reas(.>n  unknown. 

DfcnnUr  (wnks  .j9  :/2).  :t  iO.  14.90 

OdC  and  0G(.'  were  affected  during  this  period;  O.IC  weeks  49  and  ,'>2,  OOC 
weeks  .'if  and  .42.  In  order  to  improve  the  O.IC  performance  .N'T.A/Lillesf com 
raised  the  level  towards  Kjeller  by  d  .'i  dBm.  OfiC  re.sumed  operation  after  a 
Modcomp  restart  27  December. 

January  (trnks  15),  91.12.90  9.1.91 

Week  1  three  .snbarrays  were  affected:  021!  it  .lanuary  in  connection  with  a 
scheduled  niovement  to  another  I ransmis.sion  group.  4  .'>  .January  due  to  power 
failure:  02C  4  January  aud  OitC  4  .4  January,  prrdrably  c.aused  by  interrup'iisn 
related  to  N  I'.A  equipment.  In  Week  2  (X  January)  02C  data  disappeared. 
According  to  NT.A/Lillestrom  a  cable  was  damaged  in  connection  with  road 
work.  On  10  January  the  cable  was  repaired.  In  Week  .3  (IS  January)  power 
to  02B  again  failed.  After  rep.iir  and  replacement  of  the  main  fuses,  (he  .SLEW 
was  finally  reset  23  January.  The  suharray  then  again  delivered  data.  03C 
lost  synchronization  23  January.  .After  a  .Modcomp  restart,  sync  was  again 
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established. 

February  (weeks  6  9),  ). 2-3. 3. 91 

In  this  period  4  subarrays  were  affected:  02B,  02C  week  9;  03C  weeks  7 
and  9;  and  finally  06C  weeks  7,  ft  and  9.  There  wore  short  interruptions  in 
NTA  transmission  equipment/lines,  which  again  resulted  in  sync  problems. 

March  (weeks  10-13).  4-31.3.91 

The  subarrays  02B,  02C  and  OGC  were  affected  1  March,  and  03C  on  3 
March.  ,\ll  were  resynchronized  4  .March  after  a  Modcomp  restart.  A  com¬ 
munications  cable  damaged  by  an  excavation  4  March  resulted  in  loss  of  data 
from  subarray  03C.  5  March  the  cable  was  repaired. 


O.A.  Hansen 


II 


Sub- 

arrays 

Oct  90  (.5) 
1.10-4  11 

Nov  90  (4) 

5  11-2.12 

Dec  90  (4) 
3-30.12 

Jan  91  (5) 
31.12.90-3  2.91 

Feb  91  (4) 
4,2-3  3 

Mar  91  (1) 
4-31.3 

Average 
1/2  year 

1)1  A 

'^'*0.0019 

^♦0,0030 

0  0000 

0.0000 

0.0009 

0.0007 

0.004 

OIU 

^'•O.OOdO 

0.890 

0  001 

0  023 

0  0008 

''*0.002 

0.100 

0211 

*0.304 

a. 025 

0  011 

""*0  not 

'•'"*0.010 

*0.4.58 

0  145 

U2C' 

■9*0,0002 

*'•0,0002 

0,0009 

"'*0  117 

'■"*0  003 

•0.448 

0,1.50 

l)3C 

0  008 

0.004 

"'*0,002 

'-'•1)1)02 

'"*0.027 

*0,448 

0,081 

OdC 

0.003 

'••1  107 

0.001 

0.001 

0.002 

0,001 

0.200 

06C 

*0  011 

•'*2.398 

0,0000 

0.003 

''"*0,0007 

•0  343 

0.1,59 

.AVER 

0  003 

0,040 

0.002 

0.008 

0  000 

0  •213 

0,171 

•Soc  Section  2,2  regarding  figures  (ireceded  by  an  asterisk. 

figures  representing  error  rate  (in  per  centl  preceded  by  a  miml>er  1),  2).  etc  .  are  related  to 
legend  below. 


1  1,1  1  1.13) 

2),-1).7),13),M).1T) 

3l,«),iil,9),l0),15) 

51,16) 


Average  d  weeks  (•10-.13/90).  (l-S/'Jl) 

Average  3  weeks  (d2-M,d6-d8/90),  ((;-8,IO-I3/9I) 
Average  2  weeks  (d l-d2,4r)-d7,d9-r>l/90),  (2,5.8/91) 
Average  1  week  (d8/90),  (0/91) 


Table  2.2.1  Communications  performance.  The  numbers  represent  error  rates 
in  per  cent  based  on  total  transmitted  frames/week  (1  October  1990  31 

March  1991). 
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2.3  Event  detection  operation 

In  Table  2.3.1  some  monthly  statistics  of  the  Detection  and  Event  Proces.sor 
operation  are  given.  The  table  lists  the  total  number  of  detections  (DPX) 
triggered  by  the  on-line  detector,  the  total  number  of  detections  proce.ssed  by 
the  automatic  event  processor  (EPX)  and  the  total  number  of  events  accepted 
after  analyst  review  (telcseismic  phases,  core  phases  and  total). 


Total 

DPX 

Total 

EPX 

/tvTOpted  events 
P-phascs  Coro  phases 

Sum 

Daily 

Oct  90 

11300 

1381 

211 

70 

281 

.9.1 

Nov  90 

12225 

1411 

241 

64 

305 

10.2 

Dec  90 

14158 

1471 

283 

62 

345 

11.1 

Jan  91 

12849 

1354 

20C 

59 

265 

8.5 

Feb  91 

12025 

1497 

204 

43 

247 

8.8 

.Mar  91 

11200 

1409 

.329 

65 

394 

12.7 

1474 

363 

1837 

10.1 

Table  2.3.1.  Detection  and  Event  Processor  statistics,  I  October  1990  -  31 
March  1991. 


B.  Paulsen 
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3  Operation  of  Regional  Arrays 


3.1  Recording  of  NORESS  data  at  NDPC,  Kjeller 

Table  lists  the  main  outage  times  and  reasons,  and  a-s  can  be  seen  the 

main  reasons  for  the  outages  are  hardware  failure  at  the  HI  D  and  software 
failure  al  NDPC. 

The  average  recording  time  was  as  compared  to  9'2.n/{  for  the  pre- 

viou.s  period. 


Date  Time  Cause 


12 

Oct 

0910-0914 

12 

Oct 

1025-1033 

18 

Oct 

0103-0145 

19 

Oct 

0058-0151 

22 

Oct 

1238-1259 

1 

Nov 

0850-0911 

14 

Nov 

0840-0845 

16 

Nov 

0037-0045 

21 

Nov 

2034-2104 

22 

Nov 

0600-0608 

23 

Nov 

2316-2352 

24 

Nov 

0921-0943 

26 

Nov 

2029-2216 

26 

Nov 

2300- 

27 

Nov 

-0621 

28 

Nov 

1159-1354 

28 

Nov 

1430-1711 

3 

Dec 

0810-0824 

8 

Dec 

2200-2205 

18 

Dec 

0855-1029 

24 

Dec 

0534-0549 

31 

Dec 

2035-2045 

9 

Jan 

1045-1112 

11 

Jan 

1225-1228 

Service  at  HUB 
Service  at  HUB 
Timing  failure  HUB 
Timing  failure  HUB 
Service  at  HUB 
Line  failure 
Service  at  HUB 
Line  failure 
Software  failure 
Software  failure 
Software  failure 
Software  failure 
Software  fail'Te 
Software  failure 
Software  failure 
Power  failure  HUB 
Power  failure  HUB 
Hardware  failure  NDPC 
Software  failure 
Line  failure 
Line  failure 
Line  failure 
Service  at  HUB 
Service  at  HUB 


H 


r 


20 

Jan 

2101-2105 

Lino 

failure 

20 

Jan 

2112-2117 

Line 

failure 

20 

Jan 

2132-2135 

Line 

failure 

6 

Feb 

2238-2242 

Line 

failure 

2 

Mar 

0100-1325 

Software  failure 

Table  3.1.1.  Interruptions  in  recording  of  NORESS  data  at  NDPC,  1  October 
1990  -  31  March  1991. 

Monthly  uptimes  for  the  NORESS  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera¬ 
tion)  affecting  this  task  were  as  follows; 


October 

99.8% 

November 

97.8% 

December 

99.6% 

January 

99.9% 

February 

100.0% 

March 

98.4% 

Fig.  3.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  NORESS  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 


J.  Torstveit 
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3.2  Recording  of  ARCESS  data  at  NDPC,  Kjeller 


TIu’  main  roasons  causing  most  of  the  AKCKSS  outage  in  tlie  period  are: 
llardw.ire  failure  at  N'Dl'C  and  line  failuri'.  Outage  intervals  are  listed  in 
1  able  it. 2.  i . 

riie  average  recording  tiim'  was  as  coinparerl  to  fto.tiVi  for  the  pre¬ 

vious  pr'riod. 


Date  Time 


Cause 


5 

Oct 

0413-0438 

S 

Qct 

0456-0523 

5 

Oct 

0629-0723 

5 

Oct 

1321-1620 

5 

Qct 

2038-2236 

7 

Qct 

0743-0756 

10 

Oct 

0810-0911 

10 

Oct 

2138-2220 

13 

Oct 

1410-1511 

14 

Oct 

1120-1232 

15 

Oct 

1100-1133 

19 

Oct 

1603-1644 

22 

Oct 

1333-1534 

23 

Oct 

1355-1406 

25 

Oct 

0650-0834 

25 

Oct 

1555-1606 

9 

Nov 

0034-0107 

10 

Nov 

1900-1929 

13 

Dec 

0315- 

14 

Dec 

-1845 

31 

Dec 

2310- 

1 

Jan 

-0001 

7 

Feb 

0952-1025 

14 

Feb 

2141-2238 

16 

Feb 

0026-0130 

17 

Feb 

1144-1149 

17 

Feb 

1602-1657 

Power  break  HUB 
Power  break  HUB 
Power  break  HUB 
Instaling  UPS  at  HUB 
Instating  UPS  at  HUB 
Software  failure 
Software  failure 
Software  failure 
Software  failure 
Software  failure 
Software  failure 
Software  failure 
Hardware  failure 
Software  failure 
Software  failure 
Software  failure 
Software  failure 
Software  failure 
Line  failure 
Line  failure 
Line  failure 
Line  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 


IS 


17 

Fob 

1708-1738 

18 

Feb 

1259-1311 

19 

Feb 

1031-1038 

19 

Feb 

1234-1238 

19 

Feb 

1451-1507 

19 

Feb 

1556-1759 

19 

Feb 

1839-1910 

21 

Fob 

1725-1741 

21 

Feb 

2300- 

22 

Feb 

-0721 

22 

Feb 

1008-1012 

22 

Feb 

1032-1058 

22 

Feb 

1113-1123 

22 

Feb 

1144-1155 

23 

Fob 

0825-0832 

26 

Feb 

0715-0736 

6 

Mar 

1142-1209 

6 

Mar 

1214-1245 

Hardware  failure 
Hardware  maintenance 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardw2tre  failure 
hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  maintenance 
Software  maintenance 
Software  maintenance 


Table  3.2.1.  Tlie  main  ititerruptioiis  in  recording  of . A  RCES.S  data  at  NDPC, 
1  October  1090  •  31  March  1991. 

Monthly  uptimes  for  the  .4RCESS  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera¬ 
tion)  affecting  this  task  were  as  follow.s; 


October 

98.2% 

November 

99.9% 

December 

94.6% 

January 

100.0% 

February 

98.2% 

March 

99.9% 

Fig.  3.2.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  ARCESS  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 
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Fig.  3.2.1.  ARCESS  data  recording  uptime  for  October  (top),  November 
(middle)  and  December  (bottom)  1990. 
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Fig.  3.2.1.  (cont.)  ARCESS  data  recording  uptime  for  January  (top) 
ary  (middle)  and  March  (bottom)  1991. 


3.3  RecordinR  of  FINESA  data  at  NDPC,  Kjeller 

The  main  reason  for  downtime  ot  the  FINESA  array  in  the  period  arc  line 
failure,  either  between  Norsar  and  Helsinki  or  between  Helsinki  and  the  field 
installation. 

The  average  recording  time  was  98.9%  as  compared  to  .89.0%  for  the  pre¬ 
vious  period. 


Date 

Time 

Cause 

3 

Oct 

0823-0844 

Line 

failure 

3 

Oct 

0902-0907 

Line 

failure 

3 

Oct 

0917-1034 

Line 

failure 

4 

Oct 

1115- 

Line 

failure 

5 

Oct 

-0627 

Line 

failure 

15 

Oct 

1604-1615 

Line 

failure 

27 

Oct 

0241-0340 

Line 

failure 

31 

Oct 

2322- 

Line 

failure 

1 

Nov 

-0530 

Line 

failure 

10 

Nov 

1750-1757 

Line 

failure 

12 

Nov 

1014-1104 

Line 

failure 

13 

Nov 

0753-0757 

Line 

failure 

19 

Nov 

2215-2230 

Line 

failure 

20 

Nov 

0029-0046 

Lina 

failure 

19 

Dec 

2215-2230 

Line 

failure 

20 

Dec 

0029-0045 

Line 

failure 

11 

Jan 

0651-0705 

Line 

failure 

16 

Jan 

1233-1257 

Line 

failure 

19 

Jan 

0810-0817 

Line 

failure 

19 

Jan 

0846-0853 

Line 

failure 

20 

Jan 

2306- 

Line 

failure 

21 

Jan 

-0118 

Line 

failure 

24 

Jan 

0119-0744 

Line 

failure 

13 

Fob 

0832-0849 

Line 

failure 

25 

Feb 

1240-1245 

Line 

failure 

7 

Mar 

2010-2017 

Line 

failure 

12 

Mar 

2232-2238 

Line 

failure 

24 

Mar 

1006-1200 

Hardware  failure 

22 


26  Har 

27  Mar 
27  Mar 


0920-1739 

0624-0757 

1043-1349 


Hardware  failure 
Hardware  failure 
Hardware  maintenance 


Table  3.3.1  The  main  interruptions  in  recording  of  FI.NESA  data  at  NDPC, 
1  October  1990  -  ,31  March  1991, 


Monthly  uptimes  for  the  UNESA  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera 
tion)  affecting  this  task  were  as  follows: 


October 

97.0% 

November 

99..S% 

December 

9f».9% 

January 

9.'<.7%, 

February 

99.9% 

.March 

97.8% 

fig.  3. 1.1  shows  the  uptime  for  the  data  recording  t:i.sk,  or  equivalently, 
the  availability  of  EINESA  data  in  our  tape  archive,  on  a  day  by  day  biusis, 
for  the  reporting  period. 
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Fig.  3.3.1.  FINESA  data  recording  uptime  for  October  (top),  November 
(middle)  and  December  (bottom)  1990. 
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3.4  Event  detection  operation 

This  section  reports  on  results  from  simple  one-array  automatic  processing 
using  signal  processing  recipes  and  the  ’ronap’  recipes  for  the  ep  program,  as 
described  in  NORSAR  Sci.  Rep.  No.  2-88/89. 

IMS  results  are  reported  in  .Section  3.5. 

yORESS  dt'trrlions 

I  he  number  of  detections  (pha.ses)  reported  during  day  274.  1990,  through 
day  090,  1991.  was  30990,  giving  an  average  of  220  detections  per  processed 
day  (182  days  processed). 

Table  3.1.1  shows  daily  and  hourly  distribution  of  detections  for  .NOHESS. 

/Ti'di/.s  automatically  UxraUd  by  .XOllESS 

During  days  274,  1990,  through  090,  1991,  2298  local  and  regional  events 
were  located  by  NORESS.  ba,sed  on  automatic  association  of  P-  and  S-type 
arrivals.  I'his  gives  an  average  of  12.6  events  per  processed  day  (182  days 
processed).  58  %  of  these  events  are  within  300  liin,  and  87  %  of  these  events 
are  within  1000  km. 

.4  It ( 'ESS  detections 

The  number  of  detections  (phases)  reported  during  day  274,  1990.  through 
day  090.  1991,  wa.s  64359.  giving  an  average  of  354  rletcrtions  per  processed 
day  (182  days  processed). 

Table  3.4.2  shows  daily  and  hourly  distribution  of  detections  for  .ARCESS. 

Events  automatically  located  by  ARCESS 

During  days  274,  1991,  through  090,  1991,  3446  local  and  regional  events 
were  located  by  ARCESS,  ba.sed  on  automatic  association  of  P-  and  S-type 
arrivals.  This  gives  an  average  of  18.9  events  per  processed  day  (182  days 
processed).  53  %  of  these  events  are  within  300  km,  and  86  %  of  these  events 
are  within  1000  km. 

EIXESA  detections 

The  number  of  detections  (phases)  reported  during  day  274,  1990,  through 
day  090,  1991,  was  60-547,  giving  an  average  of  333  detections  per  processed 
day  (182  days  processed). 

Table  3.4  .3  shows  daily  and  hourly  distribution  of  detections  for  FINESA. 
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Events  automatically  located  by  FIXESA 

During  days  274,  1990,  through  090,  1991,  U.-CIS  Im  ,ohI  rcKuinai  im’iiis 
were  located  by  FINESA,  based  on  autoinatii  a.s>o(  iaiioii  <if  I’  and  S  ixj..' 
arrivals.  This  gives  an  average  of  21.1  events  per  processed  day  ■  isj 

processed).  72  %  of  these  events  are  within  300  ktn.  and  'X  of  these  .  ..  s’ , 

are  within  1000  km. 

GERESS  detections 

The  number  of  detections  (phases)  reported  during  day  274,  1990,  through 
day  090,  1991,  wa-s  214S0.  giving  an  average  of  185  detections  per  processed 
day  (116  days  processed). 

lable  3.4.1  shows  tlaily  and  he”ri’.  distrilMition  of  det<>clions  for  OERE.SS, 

Evtnls  mitomahctdly  Ux  uUd  by  '.'/-  /f  f  .s.s 

During  days  2!r2.  loOO,  ihront-h  090,  i'e.il.  ;|s;  hn -il  and  regional  euMil^ 
were  located  by  (iKKESS.  I)a,sed  on  .oitoinaiM  .i.-mk  latioj,  of  I’  .md  .^t  ii  pe 
arrivals.  This  givi's  an  average  of  10  7  e\ents  per  pro<e,>ed  da\  TIO  da\- 

processed).  68  V>  of  these  events  are  aitlon  300  kin.  .md  '3  '  of  these  .’wnt- 

are  within  1000  km. 
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3.5  IMS  operation 

The  Intelligent  Monitoring  System  (IMS)  was  installed  at  NORSAR  in  De¬ 
cember  1989  and  has  been  operated  experimentally  since  1  January  1990  for 
automatic  processing  of  multiple-array  data.  The  current  version  of  IMS  pro¬ 
cesses  data  from  the  two-array  network  consisting  of  NORESS  and  ARCESS. 
An  upgrade  of  IMS  scheduled  for  mid-91  will  allow  data  from  additional  arrays 
and  single  stations  to  be  incorporated. 

In  general,  our  routine  operation  of  the  IMS  during  the  reporting  period 
has  progressed  well,  and  the  system  has  proved  to  be  very  powerful  and  flex¬ 
ible.  The  well -developed  automatic  processing  combined  with  very  versatile 
interactive  tools  has  kept  the  analyst  workload  at  a  low  level,  and  in  fact  only 
one  analyst  has  been  required  to  handle  the  regular  processing  at  any  time. 

Since  the  IMS  is  still  in  an  initial  stage,  there  have  naturally  been  some 
problems  of  a  technical  nature,  but  these  have  diminished  considerably  as  the 
system  has  matured. 

Table  3.5.1  presents  an  overview  of  IMS  event  processor  downtimes.  Table 
3.5.2  gives  a  summary  of  pha.se  detections  and  prorcs.sod  regional  events  by 
IMS.  From  top  to  bottom,  the  table  gives  the  total  number  of  detections  by 
the  IMS,  the  detections  that  are  associated  with  regional  events  declared  by 
the  IMS,  the  number  of  detections  that  are  not  associated  with  such  events, 
the  number  of  regional  events  declared  by  the  IMS,  the  number  of  such  events 
rejected  by  the  analyst,  the  number  of  events  .accepted  by  the  analyst,  the 
number  of  events  accepted  by  the  analyst  without  any  changes,  and  finally  the 
number  of  events  accepted  after  some  sort  of  modification  by  the  analyst.  This 
last  category  is  divided  into  two  clas.ses:  Events  where  phases  (not  detected 
by  the  I.MS)  have  been  added  by  the  analyst  and  events  for  which  the  phase 
assignments  tjy  <he  IMS  have  been  changed  or  one  or  more  phase  detections 
have  been  removed. 

From  initial  review  of  the  IMS  bulletin,  it  is  clear  that  the  final  output  is  of 
very  high  quality  from  a  seismological  point  of  view.  More  detailed  evaluation 
will  be  conducted  at  a  later  stage. 


^4 


B.Kr.  Hokland 
U.  Baadshaug 
S.  Mykkeltveit 


1990  Oct  1  00:00 

1990  Dec  11  05:52 

1991  Feb  7  22:00 
1991  Feb  13  12:00 
1991  Feb  14  23:00 


00.000  -  1990  Oct 
10.200  -  1990  Dec 
00.000  -  1991  Feb 
00.000  -  1991  Feb 
00.000  -  1991  Feb 


2  00:00:00.000 
11  14:34:29.600 
8  08:00:00.000 
14  10:00:00.000 
16  03:00:00.000 


Table  3.5.1.  IMS  event  processor  (pipeline)  (lownlimes. 


Oct  90 

Nov  90 

Dec  90 

Jan  91 

Feb  91 

Mar  91 

Total 

Phase  detections 

18127 

19042 

13854 

1UU264 

-  Associated  phases 

2618 

2664 

2325 

2715 

2255 

2943 

15520 

-  Unassociatcd  phases 

15.193 

12914 

15802 

16327 

13397 

10911 

84744 

Events  declared 

991 

1048 

913 

1156 

957 

1126 

6191 

-  Rejected  events 

192 

229 

215 

246 

213 

174 

1269 

■  Accepted  events 

799 

819 

698 

910 

744 

952 

4922 

Unchanged  events 

239 

227 

166 

262 

227 

263 

1384 

Modified  events 

560 

592 

.532 

648 

517 

689 

3538 

Phases  added 

39 

48 

41 

48 

54 

61 

291 

Phases  changed  or  removed 

521 

544 

491 

COO 

463 

R28 

3247 

Table  3.5.2.  I.M.S  pha.se  detections  and  events  summary. 
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4  Improvements  and  Modifications 


4.1  NORSAR 

XOIISA  /?  data  acquifitiuii 

No  modification  lia.s  Ix-i-ii  made  to  llio  NOKS.Nll  data  acquisition  system. 
AO/f.S'.l/?  detection  prorctisiiig 

The  .\ORS.\R  detection  processor  lias  lieeii  running  satisfactorily  on  the 
IBM  during  this  reporting  period. 

Detection  statistics  are  given  in  section  2. 

.\OHS.\R  event  p/wt .s.'iiiy 

'I'liere  are  no  changes  in  the  routine  processing  of  NORS.VR  events,  using 
the  IBM  system. 

\ORSA  li  upgrade 

Testing  of  new  digitizer  systems  has  been  going  on.  It  is  now  clear  that 
it  v.ill  be  possible  to  transmit  sullicient  DC'  power  through  e.visting  cables  to 
operate  ritodern  21  bit  digitizers  at  the  remote  sites. 

The  strategy  for  itpgrading  the  .\'OR.S.\R  .system  is  to  continue  testing  of 
available  2!  bit  digitizer  systettis  to  evalitate  ;lieir  performance,  and  within 
the  coming  yetir  to  actpiire  continuous  data  frotn  one  subarray. 

In  |)arallel  with  hardware  testing,  therr  is  ti  need  to  reanalyze  NORS.Ml 
master  events  to  generate  a  new  datii  ba.se  of  time  debiy  corrections. 


4.2  NORESS/ARCESS/FINESA/GERESS 

De  tection  proee 

The  routine  detection  processing  of  the  array  data  is  running  sati.sfactorily 
on  each  of  the  arrays'  SCN  :t/2.s0  ac(|uisition  systems.  The  same  program 
is  u.sed  for  .\OIiS.\R,  NORl'.SS.  .XRCKSS,  I'l.M'.S.N,  and  GERESS.  but  with 
different  recipes.  The  beam  tables  for  .NORESS  and  ARC'ESS  are  found  In 
N'()RS.AR  Sci.  Rep.  No.  l-8;t/!)0.  I  he  beaut  table  for  FINE.SA/GERESS  is 
found  in  NORSAR  Sci,  Rep.  No.  |.!)U/!I1. 

Defection  stalislic.s  are  givi-n  in  .section  d. 
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During  this  reporting  period,  two  new  three-component  stations  have  been 
installed  in  Poland.  KSP  is  located  at  Ksiaz  (50.8°N,  16.3‘’E),  and  SFP  is 
located  at  Stary  Folwarh  (54.3“N,  23.3“E). 

Data  acquisition  for  the  two  stations  in  Poland  has  been  running  during 
test  periods  along  with  detection  processing.  Detection  (STA/LTA)  processing 
is  performed  on  the  vertical  component  using  13  different  filter  bands.  Incoher¬ 
ent  detection  processing  has  been  tested  on  the  horizontal  components.  The 
subsequent  signal  processing  is  still  in  an  experimental  mode  and  is  focusing 
on  the  difficult  problem  of  phase  identification. 

Event  processing.  Phase  estimation 

This  process  performs  F-K  and  polarization  analysis  for  each  detection  to 
identify  phase  velocity,  azimuth  and  type  of  phase,  and  the  results  are  put 
into  the  ORACLE  detection  data  base  for  use  by  IAS. 

Plot  and  epicenter  determination 

Descriptions  of  single  array  event  processing  are  found  in  NORSAR  Sci. 
Rep.  No.  2-88/89  and  NORSAR  Sci.  Rep.  No.  2-89/90. 
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5 


Maintenance  Activities 


5.1  Activities  in  the  field  and  at  the  Maintenance  Center 

This  section  summarizes  the  activities  in  the  field,  at  the  Maintenance  Center 
(NMC)  Hamar  and  NDPC  activities  related  to  monitoring  and  control  of  the 
NORSAR,  NORESS,  ARCESS,  FINESA  and  GERESS  arrays.  The  recently 
installed  stations  in  Poland  were  not  fully  operational  by  the  end  of  March 
1991,  but  data  from  these  stations  have  been  subjected  to  NDPC  monitoring 
during  selected  time  intervals. 

Activities  involve  preventive  and  corrective  maintenance,  and  in  addition 
installation  and  testing  of  new  equipment  related  to  satellite  communication 
(P.W.  Larsen,  Poland,  October  1990)  and  installation  and  testing  of  a  new 
UPS  system  (ARCESS,  October  1990). 

NORSAR 

This  array  was  visited  in  October  and  November  1990  and  in  January 
and  March  1991.  .“Vctivities  related  to  this  array  have  been  divcr.se,  involv¬ 
ing:  preventive  maintenance,  cable  splicing/location  of  cables,  adjustment  of 
Channel  gain  and  DC  offset,  RA-.u  amplifier  replacements,  work  on  02B  (tel) 
including  channel  VCO  adjustment,  battery  replacement,  SP/LF  instrument 
adjustments,  and  SLEM  reset  after  power  outages.  Finnaly,  there  was  work 
in  connection  with  testing  of  an  RD-3  (Remote  Digitizer)  in  conjunction  with 
the  NORSAR  upgrade  activities. 

NORESS 

The  NORESS  array  was  visited  in  November  1990  and  January  1991.  The 
satellite  clock  was  replaced,  the  UPS  reset  and  a  test  of  the  Megamux  multi¬ 
plexer  was  carried  out. 

ARCESS 

This  array  was  visited  in  October  and  December  1990  in  connection  with 
installation  of  a  new  UPS  (Uninterrupted  Power  Supply),  replacement  of  a 
seismometer  cable,  adjustment  of  fiber  optic  links  between  the  hub  and  remote 
sites,  and  restart  and  check  of  UPS  equipment  (Dec  1990). 

FINESA 

There  were  no  visits  in  the  period. 

Poland 

Satellite  communication  equipment  was  installed  (October  1990). 
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Details  on  all  actitivites  are  presented  in  Table  5.1. 


Subarray/ 

Trisk 

Date 

area 

NORSAR: 

OlA 

Preventive  maintenance  including  adjustment  of 
channel  gain,  DC  offset.  .41so  cable  work  SP02  carried 

out. 

25  Oct 

OIB 

Preventive  maintenance.  Cable  splicing  SP04 

9  Oct 

02B 

Preventive  inaintonanre. 

26  Oct 

02C 

Preventive  maintenance. 

24  Oct 

03C 

Preventive  maintenance. 

26  Oct 

04  C 

Preventive  maintenance. 

29  Oct 

06C 

Preventive  maintenance.  In  addition,  splicing  of 
cable  SP02,  04  05. 

10  Oct 

ARCESS: 

Installed  a  new  UPS.  Replaced  defective  seismometer 
cable  at  site  B5.  Adjusted  all  fiber  links  between 
the  hub  and  remote  sites. 

1-7  Oct 

Poland: 

P.W.  Lar.sen  installed  satellite  communication 
equipment. 

23-30  Oct 

NDPC: 

Daily  status  check  of  all  arrays 

Weekly  calibration  of  NORS.AR  SP/LP  instruments 
Continuous  measurement  of  Ma.ss  Position  and  Free 
Period.  Adjusted  when  outside  tolerances 

Oct 

NORSAR: 

OlA 

Cable  splicing,  SP02 

2  Nov 

Cable  splicing,  SP05 

5  Nov 

06C 

Replaced  RA-5  amplifier  SPOl 

8  Nov 

02C 

V'isited  the  subarray  in  connection  with  communication 
cable  damage 

12  Nov 

02B(tcl) 

Station  05  work,  channel  inoperative 

28  Nov 

NORESS: 

Replaced  satellite  clock.  UPS  re.sel. 

13,19 

Nov 

NDPC: 

Daily  status  clKck  of  all  arrays. 

Weekly  calibration  of  NORS.AR  SP/LP  instruments 
Continuous  measurement  of  .Mass  Position  and  Free 

Period 

Nov 

Table  5.1  Activities  in  the  field  and  the  NORSAR  Maintenance  ("enter,  in¬ 
cluding  NDPC  activities  related  to  the  NORS.Ml.  NORF.SS,  .ARCESS  and 
FINESA  arrays,  I  October  1!)9()  31  March  1991. 
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Subarray/ 

Task 

Date 

area 

ARCESS: 

Restart  of  UPS  equipment 

14,15 

Dec 

NDPC: 

Daily  status  check  of  all  arrays. 

Weekly  calibration  of  NORSAR  SP/LP  instruments 
Continuous  measurement  of  Mass  Position  and  Free 
Period 

.‘Idjustmeni  of  .\Ia.ss  Position  and  Free  Period  u  lien 
outside  tolerances 

Dec 

NORSAU; 

02n(U'i) 

Adjusted  VCO,  ch.  .I 

4  Jan 

OGC 

Adjusted  instruments 

7  Jan 

0'2B 

SLEM  restart  after  power  outage. 

23  Jan 

NORESS; 

Megamux  installed  in  connection  with  a  test. 

9  Jan 

Megamux  disconnected  after  test  completed. 

1 1  Jan 

NDPC; 

Daily  status  check  of  all  arrays. 

Jan 

Weekly  calibration  of  NOR.SAR  SP/LP  instruments 
Continuous  measurement  of  Mass  Position  and  Free 

1991 

Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

NMC: 

Work  in  connection  with  NORSAR  upgrading 

NDPC: 

Daily  status  check  of  all  arrays. 

Weekly  calibration  of  NORSAR  SP/LP  instruments 
Continuous  measuroment  of  Mass  Position  and  Free 
Period 

Adjustment  of  Ma.ss  Position  and  Free  Period  when 
outside  tolerances 

P'eb 

Table  5.1  (font.) 


•lO 


Subarray/  Tjisk  Date 

area 


NORSAR; 

06C 

An  RD-3  (Remote  Digitizer)  was  tested  on 

6,10 

channel  01 

Mar 

06C 

Replaced  RA-5  amplifier  ch  03 

14  Mar 

02B(tel) 

Replaced  channel  4  battery. 

7  Mar 

NDPC: 

Daily  status  check  of  all  arrays. 

Mar 

Weekly  calibration  of  NORSAR  SP/LP  instruments 
Continuous  measurement  of  Mass  Position  and  Free 
Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

Table  5.1  (cont.) 

5.2  Array  status 

As  of  31  March  1991  the  following  N’OR.SAR  channels  deviated  from  toler¬ 
ances. 

Ol.A  01  8  Hz  filter 
02  8  Hz  filter 
04  30  dB  attenuator 


O.A.  Hansen 
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6  Documentation  Developed 

Bannister,  S.C.,  E.S.  Husebye  and  B.O.  Ruud;  Teleseismic  P  coda  analyzed 
by  three-component  and  array  techniques;  Deterministic  location  of  to¬ 
pographic  P-to-Rg  scattering  near  the  NORESS  array,  Bull.  Seism.  Soc. 
.Am.,  80,  1969  1986,  1990. 

Dahle,  A.,  H.  Bungtim  and  I..B.  Kvamme;  Empirically  derived  PSV  spectral 
attenuation  models  for  intraplate  conditions.  Europeau  Earthq.  Eng.. 
submitted  for  publication. 

Havskov,  J.,  L.B.  Kvamme,  R.A.  Hausen,  II.  Hungum  and  C.D.  Lindholm: 
The  Northern  Norway  Seismic  .Network:  Design,  operation  and  results. 
Bull.  Seism.  Soc.  Am.,  in  press. 

Kv^rna.  T.  and  D.J.  Doornbos:  Sc.ittiTing  of  regional  P,,  by  Moho  topogra¬ 
phy,  Geophy.'i.  Res.  Lett.,  in  press. 

Mykkcltveit,  S.  (Ed.):  Bri.sie  Seismohjgical  R(s<arch.  /  October  I9S9  30 
September  1990.  Final  Report.  Kjeller.  .Norway. 

Semiannual  Tech.  Summary.  I  .Apr  -  30  Sep  1990.  NORSAR  Sci.  Rep.  No. 
T90/91,  Kjeller,  Norway. 


L.B.  Loughran 
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7  Summary  of  Technical  Reports  /  Papers  Published 
7.1  RMS  Lg  analysis  of  Novaya  Zemlya  explosion  recordings 
Introduction 

In  recent  years,  much  attention  has  focused  upon  the  use  of  the  seismic  Lg 
phase  to  determine  the  yield  of  underground  nuclear  explosions.  In  the  first  of  a 
number  of  Lg  studies  undertaken  by  the  NORSAR  staff  during  the  19S0s,  Ring- 
dal  (1983)  analyzed  digital  .\OR.SAR  Lg  data  of  selected  .Semipalatinsk  events. 
He  found  that,  when  using  NORSAR  RMS  Lg  instead  of  P  waves  recorded  at 
NORSAR  to  estimate  source  size,  it  was  possible  to  eliminate  effectively  the  mag¬ 
nitude  bias  relative  to  worldwide  mj  observed  at  .NORS.dR  between  Degelen  and 
Shagan  River  explosions.  The  method  consisted  of  averaging  log(RMS)  values 
of  individual  NORS.XR  channels,  filtered  in  a  band  of  0.6  to  3.0  Hz  in  order  to 
enhance  Lg  signal-to-noise  ratio.  RingdaJ  and  llokland  ( 1987)  expanded  the  data 
base  and  introduced  a  noise  compensation  procedure  to  improve  the  reliability  of 
measurement  at  low  SNR  values.  They  were  able  to  identify  a  distinct  P-Lg  bia.s 
between  the  Northeast  and  Southwest  portions  of  the  Shagan  River  Te.st  Site,  a 
feature  that  was  confirmed  by  Ringdal  and  Fyen  (1988)  using  Clriifenberg  array 
data.  Ringdal  and  .Marshall  (1989)  combined  P-  and  Lg  basod  source  size  esti¬ 
mators  to  estimate  the  yields  of  96  Shagan  River  explosions  from  196")  to  1988, 
using  data  on  the  cratering  explosion  of  l.u  January  196.5  as  a  reference  for  the 
yield  calculations. 

Hansen,  Ringdal  and  Richards  ( 1990)  analyzed  available  data  from  stations  in 
China  and  the  Soviet  ITiion.  and  found  that  RMS  Lg  of  Semipalatinsk  explosions 
measured  at  these  stations  showed  excellent  consistency.  They  conchided  that 
for  explosions  at  Semipalatinsk  with  good  signal-to-noise  ratio.  mjlLg)  may  be 
estimated  at  single  stations  with  an  accuracy  (one  standard  deviation)  of  about 
0.03  magnitude  unit.  It  is  noteworthy  that  this  accurary  was  consistently  ob¬ 
tained  for  a  variety  of  stations  at  very  different  azimuths  and  distances,  even 
though  the  b.asic  parameters  remained  exactly  as  originally  proposed  by  Ringdal 
for  NOR.S.MI  recordings  (0.6  3,0  Hz  bandpass  filter.  RMS  window  length  of  2 
minutes,  centered  at  a  time  corresponding  to  a  group  velocity  of  3. .5  km/s). 

In  this  paper  we  apply  Ringdal's  method  to  NORS.XR  and  (Iriifenberg  record 
ings  of  Novaya  Zemlya  explosions.  This  initial  study  focuses  on  the  Northern 
.Novaya  Zemlya  test  site,  .and  we  will  only  consider  explosions  occurring  after 
1976, 

Data 

The  data  base  for  this  study  comprises  seismic  recordings  at  NORS.XR  and 
Grafenberg  for  18  presumed  underground  nuclear  explosions  at  Novaya  Zemlya 
from  1976  through  1990. 


The  NORSAR  array  (Bungum,  Husebye  and  Ringdzd,  1971)  was  established 
in  1970,  and  originally  comprised  22  subarrays,  deployed  over  an  area  of  100  km 
diameter.  Since  1976  the  number  of  operational  subarrays  has  been  7,  comprising 
altogether  42  vertical-component  SP  sensors  (type  HS-10).  In  this  paper,  analysis 
has  been  conducted  using  data  from  these  7  subarrays.  Sampling  rate  for  the 
NORSAR  SP  data  is  20  samples  per  second,  and  all  data  are  recorded  on  digital 
magnetic  tape. 

The  Grafenberg  array  (Harjes  and  Seidl,  1978)  was  establi.shed  in  1976,  and 
today  comprises  13  broadband  seismometer  sites,  three  of  which  are  3-component 
systems.  The  instrument  response  is  flat  to  velocity  from  about  20  second  period 
to  5  Hz.  Sampling  rate  is  20  samples  per  second,  and  the  data  are  recorded  on 
digital  magnetic  tape. 

Fig.  7.1.1  shows  the  Lg  propagation  paths  from  Novaya  Zemlya  to  the  two 
arrays.  The  distance  and  azimuth  are  2200  km  and  2-56  degrees  to  NORSAR. 
compared  to  3300  km  and  213  degrees  to  Grafenberg.  Both  paths  cross  the  Bar¬ 
ents  Sea,  and  as  observed  by  several  authors  (see  Baumgardt,  1990),  this  implies 
significant  Lg  blockage  effects.  The  result  is  particularly  visible  on  NORSAR 
records,  which  show  relatively  weak  Lg  energy  compared  to  the  P  and  Sn  phases. 

Examples  of  VORSAR  recordings  of  one  of  the  explosions  are  shown  in  Fig. 
7.1.2.  We  note  that  this  (as  well  as  most  of  the  other  events  analyzed)  exhibits 
signal  clipping  of  both  P  and  Sn.  This  is  a  result  of  the  very  strong  seismic 
signals  recorded  at  NORSAR  for  Novaya  Zemlya  explosions,  in  combination  with 
the  limited  dynamic  range  of  the  digital  recording  system.  For  this  reason,  we 
have  cho.sen  to  measure  the  RMS  Lg  at  NORSAR  by  selecting  a  2-minute  window- 
in  the  Lg  coda,  starting  at  10  1/2  min  after  the  origin  time  of  the  event  (see  the 
figure).  Previous  studies  of  Semipalatinsk  explosions  have  shown  that  the  RMS 
method  is  not  very  sensitive  to  the  exact  positioning  of  the  time  window,  as  long 
as  it  is  kept  the  same  for  all  events  analyzed. 

On  Fig.  7.1.2,  we  have  also  indicated  a  two-minute  P  coda  window,  which  we 
have  used  to  calculate  NORSAR  P  coda  magnitudes  for  the  explosions,  using  the 
array  RMS  method.  The  P  coda  window  starts  6  minutes  after  the  event  origin 
time. 

In  Fig.  7.1.3  we  show  an  example  of  GRF'  recordings  of  one  of  the  explosions. 
Here,  the  dynamic  range  is  sufficient  to  avoid  any  clipping,  and  we  have  selected 
a  two- minute  window  which  includes  the  main  Lg  energy,  starting  16  minutes 
after  event  origin  time. 

Aualyxis  tr suits 

Applying  the  RMS  measurement  technique  using  our  standard  filter  hand 
(0.6  3.0  Hz)  and  averaging  over  array  elements  as  described  by  Ringdal  and 
ffokland  (1987).  we  arrive  at  results  listed  in  Tables  7.1.1  through  7.1.3. 
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Table  7.1.1  gives  our  results  for  NORSAR  P-coda  magnitudes,  using  the  time 
window  indicated  on  Fig.  7.1.2.  A  constant  correction  factor  has  been  added  to 
the  log(RMS)  values  to  make  these  coda  magnitudes  consistent,  on  the  average, 
with  world-wide  mj. 

Table  7.1.2  covers  the  NORSAR  Lg  results,  and  shows  that  RMS  Lg  can 
be  estimated  for  all  the  events  processed,  including  two  events  below  mj  =  5.0 
(events  2  and  4).  (For  the  27  Sep  78  explosion,  no  NORSAR  data  are  available.) 

Table  7.1.3  gives  corresponding  Lg  results  for  the  Grafenberg  array.  Here,  the 
smallest  of  the  events  (Event  4)  had  too  low  SNR  to  allow  reliable  measurements. 

The  ms  values  listed  in  Tables  7.1.1  through  7.1.3  are  taken  from  Lilwall 
and  Marshall  (1986)  for  events  up  to  1984,  and  have  been  calculated  from  NEIC 
station  reports  for  later  events. 

Fig.  7.1.4  shows  a  comparison  of  world-wide  ms  and  NORSAR  P  coda  mag¬ 
nitudes,  We  note  that  the  correspondence  is  excellent  (orthogonal  standard  de¬ 
viation  is  only  0.027).  Thus  the  NORSAR  recordings  appear  to  provide  a  very 
stable  measure  of  ms  for  events  from  this  test  site. 

Figs.  7.1.5  and  7.1.6  are  scatter  plots  comparing  world-wide  ms  to  NORSAR 
and  Grafenberg  RMS  Lg  magnitudes.  We  note  that  there  is  a  considerable  scatter 
in  both  of  these  plots.  In  particular,  it  appears  that  the  majority  of  events  have 
almost  the  same  M(Lg),  whereas  the  ms  values  vary  from  below  5.7  to  above  6.0 
for  this  group. 

It  is  especially  interesting  to  note  that  NORS.AR  M(Lg)  deviates  significantly 
from  ms,  whereas  NORS.AR  P  coda  corresponds  very  closely  to  ms- 

Fig.  7.1.7  shows  a  scatter  plot  of  Grafenberg  versus  NORSAR  M(Lg).  The 
correspondence  is  excellent,  with  an  orthogonal  standard  deviation  of  only  0.035. 
The  scatter  is  further  reduced  (to  0.025)  if  we  consider  only  events  with  at  least 
5  available  GRF  channels  (Fig.  7.1.8).  Thus,  we  obtain  the  same  close  corre¬ 
spondence  between  Lg  observations  from  these  two  arrays  for  Novaya  Zemlya 
explosions  as  has  previously  been  observed  for  Semipalatinsk  events. 

With  the  current  lack  of  independently  obtained  calibration  data,  it  would 
be  premature  to  draw  any  firm  conclusions  as  to  the  relative  accuracy  of  mi,  and 
M(Lg)  in  estimating  yields  of  these  explosions.  Nevertheless,  it  would  appear 
that  the  close  grouping  in  M(Lg),  especially  seen  for  the  NORSAR  data,  is  un¬ 
likely  to  be  a  coincidence.  It  would  seem  reasonable  to  conclude  that  this  group 
of  explosions  has  very  nearly  the  same  yield,  in  spite  of  the  divergence  in  ni), 
estimates.  However,  additional  analysis,  in  particular  including  available  Lg  data 
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from  Soviet  stations  for  this  event  set,  should  be  performed  in  order  to  further 
test  this  hypothesis. 


F.  Ringdal 
J.  Fyen 
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Ev 

Date  Origin  time 

mb 

NCH 

Noise 

STD 

Pccxla 

STD 

CORR 

m{ Pcoda ) 

01 

29/09/76-273  =  03.00.00.00 

5.7? 

40/42 

1.867 

0.064 

3.133 

0.049 

3.132 

5.732 

02 

20/10/76-294  =  08.00.00.00 

4.89 

40/42 

1.933 

0.066 

2.405 

0.048 

2.  378 

4.978 

03 

01/09/77-244  =  03.00.00.00 

5.71 

34/41 

1.777 

0.062 

3.151 

0.045 

3.150 

5.750 

04 

09/10/77-282=11.00.00.00 

4.51 

40/41 

1.908 

0.066 

2.139 

0.047 

2.046 

4.646 

05 

10/08/78-222:08.00.00.00 

6.04 

33/39 

1.862 

0.066 

3.352 

0.046 

3.352 

5.952 

06 

27/09/78-270=02.05.00.00 

5.68 

*00/00 

0.000 

0.000 

0.  000 

0.  000 

0.000 

0.000 

07 

24/09/79-267=03.30.00.00 

5.80 

38/39 

1.852 

0.059 

3.182 

0.050 

3.182 

5.782 

08 

18/10/79-291=07.10.00.00 

5.85 

28/39 

1.918 

0.050 

3.222 

0.050 

3.221 

5.821 

09 

11/10/80-285=07.10.00.00 

5.80 

32/38 

1.911 

0.060 

3.176 

0.041 

3.176 

5.776 

10 

01/10/81-274=12.15.00.00 

5.91 

20/39 

1.959 

0.063 

3.282 

0.037 

3.282 

5.882 

11 

11/10/82-284=07.15.00.00 

5.52 

27/40 

1.828 

0.085 

2.952 

0.050 

2.951 

5.551 

12 

18/08/83-230=16.10.00.00 

5.84 

25/39 

1.776 

0.054 

3.170 

0.054 

3.169 

5.769 

13 

25/09/83-268=13.10.00.00 

5.71 

24/39 

2.148 

0.055 

3.125 

0.052 

3.123 

5.723 

14 

25/10/84-299  =  06.30.00.00 

5.77 

28/41 

1.932 

0.063 

3.144 

0.066 

3.143 

5.743 

15 

02/08/87-214  =  02.00.00.00 

5.71 

28/40 

1.908 

0.080 

3.170 

0.048 

3.169 

5.769 

16 

07/05/88-128=22.50.00.00 

5.52 

27/40 

1.478 

0.066 

3.014 

0.038 

3.014 

5.614 

17 

04/12/88-339=05.20.00.00 

5.79 

30/40 

2.055 

0.061 

3.223 

0.046 

3.222 

5.822 

18 

24/10/90-297  =  14.58.00.00 

5.60 

38/40 

1.822 

0.070 

3.019 

0.044 

3.018 

5.618 

Table  7.1.1.  .NORS.\R  RMS  P  roRa  magiiitiidos  for  events  in  the  data  ba.se. 
The  table  lists  event  number,  origin  date  and  time,  world  wide  mj,  and  a  list  of 
measurements  made  in  this  study: 


NCH 

Noise 

STD 

Pcoda 

STD 

CORR 

m(Peoda) 


;  Number  of  NORSAR  data  channels  ii.sed.  and  the  total 
number  available 

:  Array  averaged  log  RMS  values  in  a  noise  window 
;  Corresponding  standard  deviation  across  array 
:  Array  averaged  log  RMS  values  in  the  P  coda  window 
:  Corresponding  standard  deviation 
;  Noise-forrecled  log  RMS  values  of  the  P  coda 
:  P  coda  magnitude  derived  by  adding  a  constant  term 
to  the  noise-corrected  values. 


.b" 


Ev  Date  Origin  time  mb  NCH  Noise  STD  Lg  STD  LgCORR  MLg(NAO) 


01  29/09/76-273:03.00.00.00 
02  20/10/76-294:08.00.00.00 
03  01/09/77-244:03.00.00.00 
04  09/10/77-282:11.00.00.00 
05  10/08/78-222:08.00.00.00 
06  27/09/78-270:02.05.00.00 
07  24/09/79-267:07.30.00.00 
08  18/10/79-291:07  10.00.00 
09  11/10/80-285:07  10.00.00 

10  01/10/81-274:12  15.00.00 

11  11/10/82-284:07.15.00.00 

12  18/08/83-230:16. '0.00.00 

13  25/09/83-268:13.10.00.00 

14  25/10/34-299:06.30.00.00 

15  02/08/87-214:02.00.00.00 

16  07/05/88-128:22.50.00.00 

17  04/12/88-339:05.20.00.00 

18  24/10/90-297:14.58.00.00 


5.77 

40/42 

1.867 

0.064 

4.89 

40/42 

1.933 

0.066 

5.71 

34/41 

1.777 

0.062 

4.51 

40/41 

1.906 

0.063 

6.04 

33/39 

1.862 

0.066 

5.68 

*00/00 

0.000 

0.000 

5.80 

38/39 

1.852 

0.059 

5.85 

28/39 

1 . 918 

0.050 

5.80 

32/38 

1.911 

0.060 

5.91 

20/39 

1.959 

0.063 

5.52 

27/40 

1.828 

0.085 

5.84 

25/39 

1.776 

0.054 

5.71 

24/39 

2.148 

0.055 

5.77 

28/41 

1.932 

0.063 

5.71 

28/4  3 

1.908 

0.080 

5.52 

27/40 

1.478 

0.066 

5.79 

.30/40 

2.055 

0.061 

5.60 

38/40 

1.822 

0.070 

3.161 

0.065 

3.160 

5.770 

2.479 

0.065 

2.461 

5.071 

3.147 

0.065 

3.147 

5.757 

2.278 

0.059 

2.235 

4.845 

3.174 

0.057 

3.173 

5.783 

0.000 

0.000 

0.000 

0.000 

3.170 

0.063 

3.169 

5.779 

3.128 

0.060 

3.127 

5.737 

3.175 

0.060 

3.174 

5.784 

3.173 

0.044 

3.172 

5.782 

2.994 

0.074 

2.993 

5.603 

3.197 

0.062 

3.197 

5.807 

3.189 

0.047 

3.187 

5.797 

3.196 

0.075 

3.195 

5.805 

3.197 

0.078 

3.196 

5.806 

3.109 

0.064 

3.109 

5.719 

3.191 

0.053 

3.190 

5.800 

2.996 

0.058 

2.995 

5.605 

NOTE:  The  M(Lg)  values  have  been  obtained  by  adding  a  constant  correction 
term  (2.610)  to  the  noise  corrected  log  RMS  Lg  values.  This  correction  term 
is  preliminary,  and  may  be  subject  to  later  revision. 


Table  7.1.2.  NORS.AR  RMS  l,g  magniliides  for  events  in  the  data  base.  The 
structure  of  the  table  is  analogous  to  Table  7.1.1.  The  rightmost  column  lists  the 
NORS.AR  Mfl.g)  values. 
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Ev  Date  Origin  time  mb  NCH  Noise  STD  Lg  STD  LgCORR  MLg(GRF) 


01  29/09/76-273:03.00.00.00 
02  20/10/76-294:08.00.00.00 
03  01/09/77-244:03.00.00.00 
04  09/10/77-282:11.00.00.00 
05  10/08/78-222:08.00.00.00 
06  27/09/78-270:02.05.00.00 
07  24/09/79-267:03.30.00.00 
08  18/10/79-291:07.10.00.00 
09  11/10/80-285:07.10.00.00 

10  01/10/81-274:12.15.00.00 

11  i:/10/82-284:07. 15.00.00 

12  18/08/83-230:16.10.00.00 

13  25/09/83-268:13.10.00.00 

14  25/10/84-299:06.30.00.00 

15  02/08/87-214:02.00.00.00 

16  07/05/88-128:22.50.00.00 

17  04/12/88-339:05.20.00.00 

18  24/10/90-297:14.58.00.00 


5.77 

02/04 

1.118 

0.086 

4.89 

03/04 

1.318 

0.047 

5.71 

03/04 

1.023 

0.007 

4.51 

*03/04 

1.223 

0.008 

6.04 

05/13 

1.223 

0.069 

5.68 

06/13 

1.270 

0.132 

5.80 

07/13 

1.217 

0.097 

5.85 

10/13 

1.350 

0.100 

5.80 

13/13 

1.350 

0.128 

5.91 

08/13 

1.416 

0.069 

5.52 

13/13 

1.291 

0.121 

5.84 

12/13 

1.214 

0.122 

5.71 

13/13 

1.126 

0.103 

5.77 

13/13 

1.382 

0.131 

5.71 

12/13 

1.033 

0.138 

5.52 

12/13 

1.018 

0.162 

5.79 

13/13 

1.195 

0.147 

5.60 

08/13 

1.452 

0.214 

2.025 

0.035 

2.022 

5.799 

1.435 

0.041 

1.245 

5.022 

2.C97 

0.021 

2.095 

5.872 

1.255 

0.085 

0.000 

0.000 

1.988 

0.102 

1.982 

5.759 

1.896 

0.114 

1.883 

5.660 

2.053 

0.118 

2.048 

5.825 

1.905 

0.116 

1.887 

5.664 

1.968 

0.115 

1.955 

5.732 

2.019 

0.099 

2.006 

5.783 

1.828 

0.120 

1.808 

5.585 

2.066 

0.135 

2.062 

5.739 

2.004 

0.145 

2.000 

5.777 

2.069 

0.124 

2.060 

5.837 

2.035 

0.14  7 

2.033 

5.810 

1.881 

0.148 

1.877 

5.654 

2  038 

0.128 

2.034 

5.811 

1.817 

0.150 

1.773 

5.550 

NOTE:  The  M(Lg)  values  have  been  obtained  by  adding  a  constant  correction 
term  (3.777)  to  the  noise  corrected  log  RMS  Lg  values.  This  correct!  '  term 
is  preliminary,  and  may  be  subject  to  later  revision. 


Table  7.1.3.  Grafenberg  RMS  Lg  magnitudes  for  events  in  the  data  base.  The 
structure  of  the  table  is  analogous  to  Table  7.1.1.  The  rightmost  column  li' ts  the 
Grafenberg  M(Lg)  values. 
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Fig.  7.1.1.  Map  showing  ihp  Lg  propagation  paths  from  the  main  Soviet  test 
site.s  (  Novava  Zenilya  and  Semipalalirisk )  to  the  .NOIiSAl!  ami  Grafenherg  arrays. 
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NORSAR  _ 

No'.cyc  ■■  Sv  2b5-07.'.C.C0  5.30 


Fig.  7.1.2.  K.xaniplp  of  NOHS.AR  rorordings  of  a  iN'ovaya  Zeitilya  explosion  ( 11 
October  1980).  The  center  instrument  of  each  of  the  7  subarrays  is  displayed, 
covering  2.'>  minutes  of  unfiitered  data.  The  positioning  of  lime  windows  used  for 
RMS  f.g,  F^coda  and  noi.se  measurements  is  indicated.  .Mote  the  clipping  of  the 
initial  P  and  that  ahso  the  S  |)h;i.se  is  clo,se  to  exceeding  the  dynamic  range. 


GRAFENBERG 


Fig.  7.1.3.  Kxanipic  of  (Jrafeiiborg  rcroriJing.s  of  a  .N'ovaya  Zemlya  c.xplosion 
(1  Scpirnibor  1977).  TtiP  figure  shows  20  minutes  of  unfiltered  data  from  the 
three  components  of  the  A I  seismometer  site  and  the  vertical-component  A2 
and  A.l  instrument.s.  .Note  that  the  horizontal  components  have  not  been  used 
in  our  analysis.  The  positioning  of  time  windows  used  for  RMS  Lg  and  noise 
measurements  is  indicated. 
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Fig.  7.1.4.  Plot  of  NORSAR  RMS  P  coda  nij  versus  world  wide  mj,  for  events 
in  the  data  ba.se.  The  straight  line  ha.s  been  obtained  by  least-squares  regres¬ 
sion  with  respect  to  the  horizontal  axis,  and  the  stippled  lines  correspond  to 
plus/minus  two  standard  deviations.  The  slope  (S),  intercept  (I),  orthogonal 
standard  deviation  (<t)  and  number  of  data  points  (N)  are  listed  on  the  figure. 
Note  the  remarkably  close  correspondence  between  the  two  estimators. 
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NAO  Lg 


S  =  0.70  1=1.73  ct=0.061  N  =  17 


Fig.  7.1.5.  Plot  of  NOKSAR  RMS  Lg  magnitude  versus  world  wide  ra(,.  Note  the 
much  greater  scatter  in  this  plot  compared  to  Fig.  7.1.4,  Nolational  conventions 
are  as  in  P’ig.  7.1.4. 


S=1.08  |  =  -0.44  ct=0.035  N  =  16 


4.4  4.8  5.2  5.6  6.0 

NAO  Lg 


Fig.  7.1.7.  Plot  of  (irafrnl).  p  ■  vr'-siis  NOHS.VU  KMS  Lg  magnitiuios  for  all 
ronimoii  ovfuts.  Not<’  tlie  rlosr  although  out*  point  in  particular 

(f.vriit  .*)  appears  to  !>«’  an  outlier.  Noiational  conventions  arc  as  in  F'ig.  7.1.4. 


S=1.18  l  =  -1.08  (7=0.025  N  =  13 


NAO  Lg 


Fig,  7.1.8.  Plot  of  Grafenbprg  versus  NOPSAR  RMS  Lg  magnitudes,  using  only 
events  for  which  at  lea.st  GRF  channels  were  available.  Notational  conventions 
are  as  in  Fig.  7.1.4.  Note  tha.  the  orthogonal  standard  deviation  is  as  low  as 
0,02.5.  Also  note  that  in  spite  of  the  very  small  range  of  magnitudes,  the  two 
arrays  show  mutually  consistent  trends. 
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7.2  Threshold  monitoring  of  Novaya  Zemlya: 

A  scaling  experiment 

[ntroduciion 

In  the  previous  NORSAR  Semiannual  Technical  Summary,  Kvaerna  and  Ring- 
dal  (1990)  presented  results  from  a  one-week  experiment  in  continuously  mon¬ 
itoring  the  Northern  Novaya  Zemlya  test  site.  Data  from  the  three  regional 
arrays  NORESS,  ARCESS,  FINESA  were  used  to  calculate  the  thresholds,  using 
the  method  of  Ringdal  and  Kvaerna  (1989).  The  location  of  these  three  arrays 
relative  to  the  test  site  is  shown  in  Fig.  7.2.1. 

In  that  one-week  study,  it  w<is  found  that  the  test  site  could  be  consistently 
monitored  at  a  very  low  magnitude  level  (typically  ms  =  2.5).  In  fact,  every  single 
occurrence  of  the  threshold  exceeding  mj  =  2.5  could  be  explained  as  resulting 
from  an  interfering  event  signal  either  from  leleseismic  or  regional  distance. 

While  these  results  are  very  encouraging,  there  is  clearly  much  work  remaining 
to  be  done  before  the  concept  of  threshold  monitoring  is  sufficiently  well  under- 
st<x)d.  In  this  paper,  we  attempt  to  illuminate  the  concept  further  by  describing 
a  simple  experiment,  involving  down-scaling  of  recorded  signal  traces  of  the  24 
October  1990  explosion  at  Novaya  Zemlya  and  simulating  what  might  have  been 
observed  on  the  threshold  traces  if  such  a  down-scaled  event  had  in  fact  occurred. 

Scaling  of  the  Sf  October  1990  explosion 

The  explosion  of  24  October  1990  had  a  world-wide  ms  =  5.6.  Recorded 
array  traces  of  this  event  are  shown  in  Fig.  7.2.2,  where  also  the  P-wave  SNR 
(STA/LTA)  on  each  filtered  array  beam  is  indicated.  Our  scaling  procedure 
consisted  simply  of  dividing  each  trace  by  a  factor  of  1000  and  adding  these 
down  scaled  traces  to  actually  ob.served  recordings  at  various  points  in  time. 

Two  examples  of  such  “down-.scaled”  signals  superpositioned  on  noise  are 
shown  in  Fig.  7.2.3  and  7.2.4.  The  first  of  these  figures  covers  a  “low  noise" 
interval  (local  night  time),  whereas  the  second  figure  corresponds  to  “high  noise” 
(local  day  time).  In  the  first  rase,  the  P  phase  is  readily  seen  on  all  three  arrays, 
and  the  S  phase  at  ARCES.S  is  also  prominent.  In  the  second  case,  the  phases 
are  far  less  clear,  although  the  ARCESS  P  and  S  still  have  good  SNR. 

Before  proceeding,  we  pause  briefly  to  note  that  a  down-sc2iling  by  a  factor 
of  1000  in  effect  reduces  the  event  ms  by  3  orders  of  magnitude.  In  this  sense, 
the  down  scaled  event  corresponds  to  ms  =  2.6.  We  have  not  attempted  to  apply 
any  source  scaling  law  for  signal  frequency,  partly  in  order  to  maintain  simplicity. 
Furthermore,  such  scaling  laws,  while  certainly  important,  are  not  sufficiently 
well  known  to  apply  with  any  degree  of  Ci.iifidence. 

Moreover,  it  should  be  noted  that  any  shift  toward  higher  signal  frequencies, 
as  would  be  a  natural  consequence  of  applying  frequency  scaling,  would  only  tend 
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to  improve  the  signaJ-to-noise  ratios  of  these  high-frequency  arrays.  Thus,  out 
procedure  can  be  considered  as  conservative  with  respect  to  estimating  detection 
capability. 

Simulation  of  threshold  monitoring 

Turning  now  to  tne  actual  data,  we  selected  a  typical  24-hour  time  period 
(day  104/1991),  and  added  the  down-scaled  signal  at  hourly  intervals  in  order  to 
get  a  picture  of  the  effect  under  different  noise  conditions.  A  total  of  24  identical 
signals  were  thus  added  at  different  times. 

Fig.  7.2.5  shows  the  “actual”  threshold  trace  (day  104)  for  Novaya  Zemlya, 
developed  exactly  as  described  in  detail  by  Kvaerna  and  Ringdal  ( 1990)  for  the 
one-week  monitoring  experiment.  We  note  that  there  is  only  one  peak  signifi¬ 
cantly  exceeding  ms  =  2.5;  this  corresponds  to  a  large  tele.seismic  earthquake  (mj 
=  6.0)  from  the  Ryuku  Islands. 

Fig.  7.2.6  shows  the  resulting  trace  for  that  same  day  after  adding  the  down- 
scaled  signals  and  recomputing  the  threshold  trace.  We  note  that  all  of  the  24 
occurrences  stand  out  clearly  on  the  plot.  Thus,  if  an  explosion  of  ms  =  2.6  had 
indeed  occurred  at  Novaya  Zemlya  that  day,  and  assuming  that  the  scaling  is 
representative,  there  would  have  been  clear  indications  on  the  threshold  trace  of 
such  an  explosion. 

Discussion 

We  emphtisize  that  this  study  is  only  intended  to  give  an  illustration  of  the 
potential  of  the  threshold  monitoring  method,  and  that  clearly  more  data  and 
additional  analysis  is  required  to  assess  the  situation  in  more  detail.  With  our 
procedure  of  scaling  by  a  constant  factor  in  amplitude,  we  have,  for  example, 
not  considered  signal  variance,  which  might  contribute  to  a  greater  variability  in 
the  size  of  the  amplitude  peak,  although  the  effect  is  not  expected  to  be  very 
significant. 

An  interesting  observation  is  the  way  in  which  threshold  monitoring  comple¬ 
ments  the  traditional  detection/location  type  monitoring:  Let  us  for  a  moment 
a.ssume  that  an  ms  =  2.6  explosion  had  in  fart  occurred  at  Novaya  Zemlya,  and 
that  the  resulting  signals  were  similar  to  the  scaled-down  signals  used  here.  It 
might  well  be  that  such  an  explosion  would  not  have  bcxui  detected  and  located  by 
the  regional  array  network.  In  fart,  during  daytime  noise  conditions  (Fig.  7.2.4) 
there  would  very  likely  have  been  only  one  or  two  confident  pha;e  detections  (Pn 
and  possibly  Sn  at  ARCKSS),  and  this  is  not  sufficient  to  locate  in  the  traditional 
network  sense. 

Nevertheless,  as  seen  in  this  paper,  such  an  explosion  would  have  been  clearly 
indicated  on  the  network  threshold  trace.  It  would  not  have  been  po.ssible  to 
explain  this  peak  as  resulting  from  some  “different”  event  (as  was  always  the 
case  for  such  peaks  in  the  Kvaerna  and  Ringdal  (1990)  study).  Thus,  a  peak 
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of  this  type  would  be  a  prime  candidate  for  further  detailed  off-line  anaJysis, 
possibly  implying  efforts  to  acquire  additional  data  in  order  to  further  elucidate 
the  nature  of  the  event. 


T.  Kvaerna 
F.  Ringdal 
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LATITUDE  (DEG  N) 


Novaya  Zemlya 


Fig.  7.2.1,  Location  of  the  target  area  (Novaya  Zemlya)  for  the  threshold  mon¬ 
itoring  experiment.  The  locations  of  the  three  arrays  NORESS  (A  =  2280  km), 
ARCESS  (A  =  1110  km)  and  FINE.SA  (A  =  1780  km)  are  indicated. 
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Novaya  Zemlya  explosion  24  October  1990,  mo  =  5.6  -  Original  recordings 
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Fig.  7.2.2.  P-  and  S-wave  recordings  (filtered  array  beams)  at  ARCESS,  FINESA 
and  NORESS  for  the  Novaya  Zemlya  nuclear  explosion  of  24  October  1990.  The 
SNRs  of  the  detecting  P-beams  are  also  given. 
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Fig.  7.2.3.  “Down-scalrci”  signals  from  (ho  Novaya  Zornlya  nuclear  explosion  of 
24  October  1990  superiinposetl  on  noise  during  a  “low  noise"  intei  'al.  The  origin 
time  of  the  simulated  “event"  is  1991/101/00.30.00. 
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ARCESS 
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Novaya  Zemlya  explosion  scaled  to  mi,  =  2.6  --  Daytime  noise  conditions 


Simulated  recordings  scaled  signal  added  to  data  from  I99>  '04  It. 30 


0/'  "3  ■  1  H  ,1 


Fig.  7.2.4.  Sana’  a.s  Vig.  f  .2.11.  liiil  for  a  "higli  noi.so"  intorvai.  The  origin  time 
of  the  simiilatod  ■•.■vorif  i^  1 1)‘)  1  /  1 01  /  I  1  .ItO.OO. 
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Continuous  threshold  monitoring  of  Novaya  Zemlya  -  Day  104,  1991 


a)  Observed  threshold  traces 
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Fig.  7.2.5.  Threshold  monitoring  of  the  Novaya  Zemlya  lest  site  for  day  1991/lOd 
(M  April  1991).  The  trip  three  traces  represent  thresholds  (upper  90  per  cent 
magnitude  limits)  obtained  from  each  of  the  thr<?e  arrays  (ARCESS,  FINESA, 
NORESS),  whereas  the  bottom  trace  shows  the  combined  network  thresholds. 
.Note  that  for  the  network  trace  there  is  only  one  magnitude  peak  exceeding  2.b. 
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Continuous  threshold  nrtonitoring  of  Novaya  Zemtya  -  Day  104,  1991 
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—  b)  Simulated  threshold  traces:  rriti  -  2.S  explosion  added  at  hourly  Intervals 
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Fig.  7.2.6.  .Same  as  tig.  7.2..5.  but  with  down-.sraied  signal.s  suporimpo.sed  on 
the  data  at  hourly  interval.s.  Nolo  that  all  orcurronce.s  of  the  simulated  04  =  2.6 
event,?  clearly  .stand  out  on  the  combined  network  trace. 
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7.3  Current  status  of  development  of  the  regional  network  asso¬ 
ciated  with  the  NORSAR  Data  Processing  Center 

The  purpose  of  this  contribution  is  to  summarize  the  status  of  development  and 
future  plans  for  the  regional  network  in  northern  Europe  that  contributes  seis¬ 
mic  data  in  real  time  to  the  NORSAR  Data  Processing  Center.  The  network  is 
shown  in  Fig.  7.3.1  and  currently  comprises  the  NORESS  and  ARCESS  arrays  in 
Norway,  the  FINESA  array  in  Finland,  the  GERESS  array  in  Germany,  and  the 
two  3-component  stations  at  Ksiaz  and  Stary  Folwaik  in  Poland.  Also  summa¬ 
rized  in  this  contribution  is  the  current  status  of  development  of  the  Intelligent 
Monitoring  System  (IMS)  and  plans  for  the  near  future. 

The  new  3-coriiponent  stations  in  Poland 

A  description  of  the  two  3-component  stations  at  Ksiaz  and  Stary  Folwark  in 
Poland  is  given  by  Mykkeltveit  and  Paulsen  ( 1990).  The  current  system,  compris¬ 
ing  field  installations  in  Poland  and  associated  telecommunications  arrangements 
for  real  time  transmission  of  data  to  NORSAR,  is  shown  schematically  in  Fig. 
7.3.2.  The  system  is  fully  operational  as  of  April  1991,  and  will  enable  Poland  to 
take  an  active  part  in  the  GSETT-2  (Group  of  Scientific  Experts’  Technical  Test 
number  2)  experiment  during  22  April  -  2  June  1991. 

During  the  fall  of  1991,  the  telecommunication  links  will  be  rearranged  to 
include  also  a  satellite  ground  station  in  Warsaw  for  real  time  reception  of  data 
from  the  two  stations  in  Poland.  A  Sun  Sparcstation-based  data  acquisition  and 
processing  system  is  also  planned  for  installation  at  the  Institute  of  Geophysics 
in  Warsaw.  It  is  expected  that  this  will  effectively  contribute  to  the  broadening 
of  the  scientific  cooperation  between  NORSAR  and  the  Institute  of  Geophysics  in 
Warsaw.  Such  cooperation  is  needed  in  order  to  acquire  relevant  information  on, 
e.g.,  seismicity  and  wave  propagation  characteristics  in  Poland  and  surrounding 
areas,  for  integration  into  the  IMS  knowledge  base. 

The  NORESS,  ARCESS,  FINESA  and  GERESS  arrays 

A  comprehensive  description  of  NORESS  and  ARCESS  is  given  in  Mykkeltveit 
et  o/(1990).  These  array  have  been  in  stable  and  continuous  operation  since  they 
were  installed  in  1984  and  1987,  respectively.  The  uptime  statistics  provided  in 
the  present  and  past  issues  of  the  NORSAR  Semiannual  Technical  Summaries 
testify  to  this.  There  are  no  plans  for  any  significant  modifications  to  these 
arrays. 

The  performance  of  the  somewhat  smaller,  technically  less  sophisticated,  yet 
very  powerful  FINESA  array  in  Finland  has  recently  been  described  by  Uski 
(1990).  Considering  the  simplicity  of  the  FINESA  field  installation,  its  opera¬ 
tional  stability  since  the  upgrade  of  the  data  acquisition  system  in  December  1989 
has  been  remarkable.  There  are  no  immediate  plans  for  modifying  the  FINESA 
system. 
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The  GERESS  array  in  German  Bavaria  has  been  described  by  Harjes  ( 1990). 
Results  from  the  processing  of  GERESS  data  at  NORSAR  have  been  presented 
by  Even  (1990).  Although  the  quality  of  data  received  at  NORSAR  is  not  yet 
entirely  satisfactory,  the  data  are  being  processed  continuously  and  also  used 
experimentally  by  IMS  (see  below).  The  GERFISS  field  system  developer  is  cur¬ 
rently  concentrating  on  solving  remaining  technical  problems.  Cooperatit'e  efforts 
between  NORSAR  personnel  and  scientists  from  the  Ruhr  University  in  Bochum, 
Germany,  currently  focus  on  optimizing  the  GERESS  beam  deployment.  Again, 
active  cooperation  is  needed  for  the  purpose  of  supplementing  the  IMS  knowledge 
base. 

Data  from  all  four  arrays  will  be  contributed  to  the  GSETT-2  experiment, 
along  with  data  from  about  50  other  single  stations  and  arrays  worldwide.  This 
will  proi'ide  another  excellent  opportunity  to  assess  the  capability  of  NORESS- 
type  arrays  for  detection  of  weak  seismic  events  at  both  regional  and  teleseismic 
distances. 

The  Intelligent  MoTiitoring  Hyntevi 

IMS  is  a  system  for  joint  processing  of  data  from  a  regional  network  of  arrays 
and  single  3-component  stations.  IMS  has  bc*en  described  in  detail  by  Bache  el 
al  (1990),  and  initial  results  from  operating  the  system  are  given  by  Bratt  et  al 
( 1990).  IMS  is  distributed  between  NORS.AR  and  the  Center  for  Seismic  Studies 
(CSS)  in  Arlington,  V'irginia,  as  indicated  in  Fig.  7.3.2. 

The  first  version  of  IMS  provides  for  joint  processing  of  data  from  NORESS 
and  ARCESS.  This  version  has  been  in  operation  at  NORSAR  since  January 
1990,  and  event  statistics  are  reported  in  the  Semiannual  Technical  Summaries. 
The  analysis  at  NORSAR  of  regional  events  for  the  GSETT-2  experiment  is 
carried  out  using  IMS  in  its  current  version. 

1  he  IMS  system  developer  S.AIC  is  currently  operating  an  upgraded  version 
of  I.MS  at  CSS.  'I  his  now  version  allows  processing  of  data  from  an  arbitrary 
number  of  arrays  and  single  3-roniponent  stations.  Since  March  1991,  data  from 
NORESS,  ARCESS,  f  lNES.A  and  GERESS  are  jointly  and  experimentally  pro¬ 
cessed  at  CSS.  According  to  current  plans,  this  new  version  of  IMS  will  be  in¬ 
stalled  at  NORS.AR  during  the  summer  of  1991. 


S.  Mykkeltveit 
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Fig.  7.3.1.  1  hp  figure  shows  the  network  of  regional  arrays  and  single  3- 

romponent  stations  in  northern  F.iirope  conti  huting  real  time  data  to  NORS.\R. 
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Fig.  7.3.2.  The  diagram  sliows  llir  two  stations  in  I’oiand  with  a-ssoriatod  data 
commnniration.s  arrangrinonts.  Fho  diagram  also  slutws  how  data  from  thnst' 
stations  art’  mad<‘  available  to  t ho  Intelligent  Monitoritig  System,  both  in  Norway 
and  the  I'.S. 
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7.4  Multichannel  statistical  data  processing  algorithms  in  the 
framework  of  the  NORSAR  event  processing  program  pack¬ 
age 

The  .N'ORSAR  interactive  data  proces.siiig  package  w;i.s  developed  for  the  anal¬ 
ysis  of  small  aperture  seismic  array  observations.  This  package,  called  the 
Event  Processor  ( EP).  has  turned  out  to  be  a  very  convenient  tool  for  the  daily 
production  of  a  seismic  bulletin.  The  small  aperture  seismic  arrays  NORES.S. 
.\RrESS.  FI.\ES.\  and  tlERESS  and  their  associated  data  processing  facili¬ 
ties  are  constructed  for  automatic  recording,  location  and  classification  of  low 
magnitude  regional  events  and  medium  magnitude  leleseismic  events.  Sig¬ 
nal  detection  is  performeil  in  an  online  mode,  whereas  parameter  estimation 
can  be  performed  as  an  oflline  proceduri’  using  i<'Cord<'d  multichannel  seismic 
wavetrains.  .Now,  at  .NORS.M!  the  automatic  system  is  in  full  operation,  pro¬ 
viding  seismic  signal  detection  a.s  well  as  signal  parameter  estimation.  Those 
parameters  are;  onset  time,  a/imiith.  apparent  velocity,  dominant  freqiiencv. 
signal-to-noise  ratio  an<l  so  on. 

This  automatic  system  is  operated  at  a  low  threshold  and  inevitably  jtro- 
duces  numerous  "fal.se  alarms”,  i.e.,  •'events"  c.'iused  by  nf)ise  hursts.  The 
current  seismic  bulletin  is  issued  after  an  interpretation  of  the  detection  list, 
which  may  be  effectively  performed  using  the  interactive  Event  Processing 
(EP)  package.  The  EP  program  has  a  nundrer  of  graphic  routines  and  inter¬ 
faces  for  work  with  NORS.NR  data  basi's.  It  also  comprises  some  sophisticated 
and  rather  time-consuming  data  processing  programs  which  at  present  cannot 
be  used  in  automatic  or  online  processing  moiles.  i.e..  without  human  con¬ 
trol.  In  particular,  adaptive  statistical  multichannel  data  proce.ssing  programs 
have  In-en  installed  recently  in  the  framework  of  the  EP  package.  These  pro¬ 
grams  are  based  on  optimal  methods  of  multidimensional  time  series  analysis 
described  in  Kushnir  <1  ul  (lOHd).  Knshnir  and  Lapshin  (1!)8-1),  Pisarenko 
f/  <il  (inttVl  and  Knshnir  tl  nl  (lOOOa.  1901)1)).  They  comprise  the  fcRowing 
procedures  of  .seismir  array  data  processing: 

1 ,  .Selection  of  the  array  inst  rnmiuils  and  data  lime  interval  to  be  proci-ssed. 

'?  EilliTing  and  ri*sami)liiig  of  tin' data. 

3.  rime  shifting  of  the  )  haniiel  waveforms  with  delays  corresponding  to  a 
given  azimuth  ami  appari'ut  velority  of  a  plane  wave  propagating  across 
the  array. 

1.  Summing  of  shifted  waveforms,  i.e..  beamforming  for  a  given  azimuth 
and  apparent  velocil'. 

.').  Whitening  of  the  lu.rkgrouud  noise  by  arlaptive  filtering  of  the  beam 
output.  This  procedure  provides  sigual-to-noise  ratio  (SNR)  gain  due 


to  differences  in  signal  and  noise  frequency  contents,  but  it  distorts  the 
signal  waveform. 

6.  Adaptive  optimal  group  filtering  (multichannel  Wiener  filtering)  of  the 
array  seismograms.  This  procedure  permits  high  suppression  of  seismic 
noise  due  to  its  coherency.  Theoretically  it  provides  maximum  SNR  gain 
without  any  distortion  of  seismic  signal  waveform  or  frequency  contents. 

7.  Adaptive  detection  of  the  distinct  phases  in  single  channel  seismograms 
or  in  the  output  traces  of  beamforming  and  optimal  group  filtering  pro¬ 
cedures.  The  detection  procedure  takes  into  account  the  difference  be¬ 
tween  noise  on  one  hand,  and  signal  plus  noise  on  the  other,  not  only  in 
amplitude  values  but  also  in  power  spectra. 

8.  Seismic  wave  onset  time  estimation  by  detection  of  the  moment  in  time 
when  the  wavetrain  statistical  features  are  abruptly  changed. 

Array  data  processing  using  the  procedures  listed  above  is  accomplished 
in  the  framework  of  the  EP  system  with  the  aid  of  a  specially  developed  set 
of  commands.  The  major  commands  are  named  GRFADAPT,  GRFFILT, 
ESTDET,  ESTONl,  ESTON3.  Description  of  these  commands  and  examples 
of  performance  are  given  below. 

GRFADAPT  and  GRFFILT  commanils 

The  procedures  1-6  are  carried  out  sequentially  by  each  of  two  EP  sys¬ 
tem  commands:  GRFADAPT  and  GRFFILT.  The  GRF.\DAPT  command, 
unlike  the  GRFFILT  command,  contains  additional  adaptation  algorithms 
which  have  not  been  listed  above.  These  algorithms  are  used  before  the  exe¬ 
cution  of  procedures  5  and  7  and  accomplish  adaptation  to  the  current  noise  of 
the  beam  whitening  filter  and  the  optimal  group  filter.  During  GRFADAPT 
command  execution,  the  processed  data  are  regarded  as  “pure”  seismic  noise 
and  the  autoregressive  (AR)  model  of  the  beam  noise  time  series  and  the  mul¬ 
tidimensional  AR  model  of  the  array  noise  time  series  are  estimated.  Based  on 
these  models,  the  whitening  and  optimal  group  filter  coefficients  are  evaluated. 
They  are  stored  and  used  later  during  GRFFILT  command  execution, 

GRFADAPT  and  GRFFILT  command  execution  produces  seven  output 
traces.  The  first  four  of  them  are  the  main  resulting  traces  and  the  last  three 
-  auxiliary  traces  —  are  needed  to  check  the  adaptation  quality.  These  seven 
output  traces  are  placed  on  top  of  the  EP  system  data  stack  (containing  all 
input  and  output  time  series  during  the  data  processing).  The  main  traces 
are; 


1.  Beam  waveform  composeil  of  filtered  and  resampled  channel  traces  for 
a  given  azimuth  and  apparent  velocity  (OGF  beam). 
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2.  Whitened  beam  waveform  (OGF  Wbeam), 

3.  Optimal  group  filter  waveform  calculated  using  filtered  and  resampled 
channel  traces  as  input  for  a  given  a2imuth  and  apparent  velocity  (OGF 
t-un), 

4.  Whitened  optimal  group  filter  waveform  (OGF  th-w). 

The  GRFADAPT  and  GRFFILT  procedures  also  calculate  the  mean  values 
and  variances  of  the  listed  output  traces  and  of  the  input  channel  traces. 
The  important  result  of  these  calculations  is  the  value  of  the  ratios  of  the 
adapative  optimal  group  filter  ( AOGF  )  output  variance  to  the  beam  waveform 
variance  and  AOGF  output  variance  to  the  averaged  channel  trace  variance. 
After  GRFADAPT  command  execution  the  first  ratio  characterizes  the  relative 
noise  suppression  by  beamforming  and  optimal  group  filtering.  Due  to  the 
signal  undistorting  feature  of  these  procedures  (provided  a  plane  seismic  signal 
wave  is  arriving  with  the  given  azimuth  and  apparent  velocity),  this  ratio  also 
characterizes  the  relative  SNR  gain  due  to  beamforniing  and  the  AOG  filtering. 

Using  the  GRFADAPT  and  GRFFILT  commands,  a  report  file  is  created 
containing  the  input  and  output  numerical  parameter  values  and  description 
of  the  processed  channel  traces.  Particularly,  this  report  contains  the  value  of 
AOGF  and  the  beamforming  SNR  gain  ratio.  .An  example  of  such  a  report  is 
shown  in  Fig.  7.4.1. 

The  format  of  the  GRFADAPT  command  is  given  below: 


grf adapt  vel  [apparent  velocity,  km/sec]  azi  [azimuth,  degrees] 
{filter  typo}  [cutoff  frequencies,  Hz]  factor  [resampling  factor] 

where  {filter  type}  is  one  of  the  three  character  strings:  Ip,  bp,  hp,  which 
means  low-pass,  band-pass  and  high-pass  filter  types. 

Before  the  GRFADAPT  command  is  initiated,  values  of  the  associated 
parameters  have  to  be  assigned.  There  are  additional  numerical  parameters 
which  are  not  specified  with  the  command,  and  which  have  the  following  de¬ 
fault  vaules: 

1.  Filter  frequency  response  decay  factor:  ALPHA  =  10"^, 

2.  Filter  impulse  response  one-side  length:  IRL  =  15; 

3.  Order  of  beam  noise  AR  model  (number  of  beam-whitening  filter  coeffi¬ 
cients):  DARB  =  10, 

4.  Number  of  input  array  data  matrix  autocovariance  coefficients:  LCRC 

=  6, 
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5.  Regularizator  of  matrix  autocovariance  function:  REG  =  10 

6.  Order  of  input  array  data  multidimensional  AR  mode’  (one-side  length 
of  the  optimal  group  filter):  DARGRF  =  6, 

7.  Auxiliary  parameters:  DARARF  =  10,  DMARF  =  20. 

Assigning  of  alternative  values  for  the  parameters  listed  above  can  be  done 
by  the 

EP  command: 


gr  [parameter  name]  [parameter  value] 

To  check  the  current  parameter  setting  before  GRF.ADAPT  and  GRFFILT 
command  execution,  one  should  use  the  EP  command: 


q  gr 

There  is  no  need  to  enter  any  parameter  values  before  the  GRFFILT  com¬ 
mand  execution.  The  computations  are  carried  out  with  numerical  parameter 
values  and  wihitening  and  optimal  group  filter  coeflicients  stored  during  the 
previous  GRFADAPT  command  cxeuction. 

The  numerical  results  of  the  GRF.AD.IPT  command  execution  are  written 
to  the  disk  file  GRFREPORT.OUTPl’T  by  the  command: 

grf report 

The  purpose  of  the  array  data  proce.ssing  using  the  GRFADAPT  and  GRF¬ 
FILT  procedures  is  to  compute  the  adaptive,  statistically  optimal  beam  which 
suppresses  coherent  and  incoherent  array  noise,  thus  providing  the  maximum 
SNR  gain  without  distortion  of  the  signal  waveform.  The.se  procedures  are  es¬ 
pecially  efficient  in  the  case  when  the  signal  and  coherent  noise  power  .spectra 
are  overlapping.  In  this  rase  they  ran  provide  much  '  cr  ''NR  gain  than 
bandpa.ss  filtering  after  conventional  beamforming. 

For  this  reason,  in  the  first  experiments  with  the  (I'RFAD.'tPT  and  GRF¬ 
FILT  procedures  in  the  framework  of  the  EP  system,  we  tried  to  learn  how  the 
program  parameter  values  influence  the  .AOGF  SNR  gain  relative  to  conven¬ 
tional  beamforming  SNR  gain,  when  data  are  processed  in  a  broad  frequency 
band.  The  main  program  parameters  which  inlliience  the  ipiality  of  the  opti 
mal  group  fdter  adaptation  are  I )  the  order  of  the  inpui  data  multidimensional 


AR  model,  2)  the  data  frequency  band  and  3)  regularizator  of  the  data  ma¬ 
trix  covariance  function.  Table  7.4.1  comprises  the  results  of  NORESS  noise 
processing.  The  120  sec  interval  of  array  noise  shown  in  Fig.  7.4.2  has  been 
used  for  optimal  group  filter  adaptation. 

As  one  can  see  from  this  table,  increasing  the  input  data  multidimensional 
-AR  model  order  leads  to  a  strong  increa.se  of  AOGF  SNR  gain.  Neverthe¬ 
less,  it  does  not  seem  worthwhile  to  use  a  multidimensional  AR  model  order 
greater  than  10  12  becaii.se  the  adaptation  procedure  becomes  time  consuming 
and  less  stable  (especially  when  a  large  number  of  array  channels  are  used). 
Choosing  a  higher  frequency  band  leads  to  diinini.shing  of  AOGF  SNR  gain. 
This  can  be  explained  by  the  strong  coherency  of  the  .NORESS  noise  mainly 
at  low  frequencies.  Varying  the  regnlari/.ator  value  from  0  to  10'"*  practically 
does  not  affect  the  AOGF  S.NK  gain. 

The  GRFADAl’T  and  GRFFILT  commands  were  also  tested  by  processing 
some  low  magnitude  local  event  records  from  the  .NORESS.  ARCESS  and 
FINESA  arrays.  In  these  experiments  the  optimal  group  and  whitening  filter 
adaptations  were  made  ttsiitg  the  array  noise  records  preceding  the  seismic 
signal  wavetrairis.  Tlte  duration  of  the  tioise  titne  intervals  used  for  adaptation 
were  front  100  to  l.hO  sec.  Tlte  main  purpose  of  llte.se  experitnents  was  to  learn 
about  possible  differences  in  P.  S  atid  Lg  wave  pha,se  extraction  by  different 
array  data  processing  algorillims  such  as  conventional  beamforming,  beam 
output  noi.se  whitening,  adaptive  optimal  group  filleritig  and  .AOGF  output 
noise  whitening.  What  dislingttishes  these  experiments  from  tliose  described 
in  our  previous  report.?  (Ku.slinir  (t  al,  1000a;  Ktishiiir  cl  nl.  1900b)  is  the 
processing  of  array  data  in  dilTorent  frequency  band.?;  0.2  5  Hz,  0.2-10  Hz. 
atid  0.2  -20  Hz.  Table  7.4.2  comprises  f  lie  ratio.?  of  z\OGF  out  put  SNR  relative 
to  heatnforming  output  SNR  and  AOGF  output  S.NR  relative  to  averaged 
channel  SNR  for  different  phases  being  extracted  from  7  .small  local  event 
records.  Each  ptiase  lias  lieeii  extracted  from  the  noise  by  the  conventional 
beam  and  AOGF  adjusletl  for  tlie  azimnlli  and  apparent  velocity  of  this  phase 
arrival  as  given  in  tlie  NORS.AR  delection  list.  One  can  see  that  in  these 
experiments  the  .AOGF  SNR  gain  relative  to  ilie  cotivent ional  beamforming 
gain  was  betwen  10.2  and  24. -a  dl)  for  tlie  0.2  .7  Hz  frequency  band,  between 
12..'1  and  21.2  dll  for  the  0.2-10  llz  frequency  band  and  around  10-11  dB  for 
the  0.2  20  Hz  frequency  band.  Note  that  the  highest  .AOGF  SNR  gain  of  more 
than  24  dll  was  acliieved  on  the  FINE.NA  records.  Tliis  is  due  to  the  presence 
in  these  records  of  strong,  highly  coherent,  low  frequeiirv  background  noise 
possilily  caused  by  stormy  .seashore  waves.  This  noise  has  been  suppressed  by 
the  AOGF  procedure,  but  not  by  conventional  beamforniing. 

Flxampies  of  GRF.AOAPT  and  GRFFII-T  oiilpiif  traces  as  the  result  of 
event  wavetrain  and  (irerediiig  noise  proce.ssing  are  given  in  Figs.  7,4.3  7.4. .5. 
In  some  of  these  examples,  the  .AOGF  and  whitened  beam  output  wavetiains 
are  similar.  One  may  conclude  that  the  adaptive  whitening  procedure  after 
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conventionaJ  beamforming  can  provide  the  same  results  as  the  AOGF  pro¬ 
cedure  (while  being  less  time  consuming).  But  this  inference  is  true  only  for 
body  waves  of  local  low  magnitude  events,  which  as  a  rule  have  high  frequency 
contents.  Power  spectra  of  surface  waves  and  teleseismic  body  waves  are  often 
overlapping  with  those  of  coherent  seismic  noise.  In  this  case,  the  AOGF  pro¬ 
cedure  has  a  strong  advantage  over  other  filtering  procedures  since  it  r^^ains 
the  signal  undistorted. 

ESTDET  command 

Adaptive  detection  of  distinct  phases  in  the  wavetrain  is  accomplished  by 
the  EP  command  ESTDET.  This  piocedure  is  based  on  the  optimal  statistical 
algorithm  described  in  Pisarenko  el  al  (1987).  Du  ,ng  the  ESTDET  command 
execution  the  time  interval  of  the  data  being  processed  is  divided  into  two 
parts.  The  data  in  the  first  internal  are  regarded  as  “pure”  noise  and  its 
AR  model  is  estimated  (“noise  AR  model”).  The  data  in  the  second  interval 
are  presumed  to  contain  the  seismic  phases.  The  detection  of  these  phases  is 
carried  out  using  a  moving  lime  window.  The  detection  algorithm  consists  of 
calculation  of  the  simplified  Bayesian  test  statistic  for  the  hypothesis:  a)  the 
AR  mode!  of  time  series  inside  the  moving  time  wi  dow  is  the  same  as  the  noise 
AR  model  versus  the  hypothesis:  b)  these  two  models  are  different  (Kushnir  ,  i 
al,  1983).  The  ESTDET  program  takes  as  input  the  wavetrain  at  the  top  of  the 
EP  system  data  stack  and  produces  nine  new  traces  which  in  turn  are  placed 
on  the  top  of  ttiis  data  stack.  Eight  of  these  traces  contain  the  values  versus 
time  of  the  detection  statistics  calculated  using  data  in  a  moving  time  window. 
These  statistics  are  derived  from  four  slightly  different  versions  of  the  simplified 
Bayesian  test  described  above.  The  first  four  traces  arc  the  statistical  '  alues 
in  logarithmic  scale,  the  next  four  are  the  same  values  in  linear  ,'^^^ale.  The 
ninth  output  trace  is  the  auxiliary  trace  for  noise  .AR  modelling  checking.  The 
detection  triggering  of  the  .seismic  phases  is  now  performed  in  an  interactive 
mode  by  comparing  the  detection  .statistic  value  with  the  thre.shold  rho.sen  to 
provide  the  acceptable  false  alarm  rate.  But  it  would  <.  leariy  be  straightforward 
to  develop  a  special  EP  command  for  automatic  phase  detection  triggering. 

The  ESTDET  command  format  is 


estdet  start  [first  point,  sec.]  end  [last  point,  sec.] 
w  [window  length,  sec.]  o  [AR  model  order]  noise 
[noise  interval  length,  sec.] 

where  “start”  and  “end”  are  the  first  and  last  points  of  the  trace  hei;.g  pro¬ 
cessed  (in  sec,  relative  to  the  initial  point  of  this  Iraie);  “w“  is  the  width  of 
the  moving  time  window,  “noise”  is  the  length  of  the  first  part  of  the  data 
time  interval  n.sed  for  the  noise  .Ml  model  estimation. 

Examples  of  E.S'I'DET  command  output  Irace.s  and  input  wavetrains  are 
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given  in  Figs.  7.4.6-7.4.8. 

ESTONl  and  ESTON3  commands 

The  moving  window  detection  procedure  points  out  those  time  intervals 
where  seismic  wave  phases  are  present.  The  next  stage  of  the  signal  processing 
is  an  estimation  of  phase  parameters.  Among  the  most  important  parameters 
needed  for  event  source  location  are  the  wave  phase  onset  times.  For  rough 
estimates  of  these  times,  the  moments  of  detection  triggering  may  be  used.  But 
for  precise  estimation  of  each  phase  onset  time,  special  statistical  procedures 
have  been  developed  (Pisarenko  el  al,  1987).  In  the  framework  of  the  EP 
system,  this  procedure  is  realized  as  two  commands;  ESTONl  and  EST0N.3. 
The  first  command  is  intended  for  single  component  trace  processing.  This 
may  be  the  beam  or  AOGF  output  trace  or  a  “raw”  array  single  channel 
wavetrain.  The  second  command  is  intended  for  S  component  seismogram 
processing  with  the  purpose  of  on.set  time  estimation.  Both  procedures  arc 
based  on  maximum  likelihood  algorithms  for  estimation  of  the  moment  in 
time  when  the  time  series  AR  model  parameter  values  abruptly  change.  The 
ESTONl  procedure  takes  into  account  changing  of  the  time  series  variance 
and  frequency  contents  at  the  moment  in  time  when  the  .seismic  pha.se  arrives. 
The  EST0N3  procedure  also  takes  into  account  changs  in  the  3-dimensional 
time  series  polarization  features  at  this  moment. 

Both  commands  use  the  traces  at  the  top  of  the  EP  system  data  stack  as 
input:  ESTONl  takes  the  upper  trace,  ESTON3  —  the  three  upper  traces. 
Both  commands  produce  one  output  trace  containing  the  onset  time  likelihood 
function  calculated  for  data  inside  a  given  time  interval.  After  command 
execution,  this  trace  is  placed  on  the  top  of  the  EP  data  stack. 

The  ESTONl  and  ESTON3  command  format  is: 


8stonl(3)  start  [first  point,  sec.]  end  [last  point,  sec.] 
w  [ciin,  window  length]  o  [AR  model  order] 


where  “start”  and  “end”  are  the  first  and  last  points  of  the  wavetrain  being 
processed  (in  sec.  relative  to  the  initial  point  of  the  trace);  “w”  is  the  minimum 
width  of  the  data  window  for  the  AR  model  estimation.  The  onset  time 
likelihood  function  is  calculated  for  the  data  inside  the  time  interval  (start  -|- 
w.  end  -  w)  (in  sec.  relative  to  the  trace  initial  point). 
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Examples  of  ESTONl  output  traces  for  different  types  of  seismic  wave 
phases  are  given  in  Figs.  7.4.9-7.4.11.  The  onset  time  values  given  in  these 
figures  as  the  arguments  of  the  likelihood  function  absolute  maximums  coincide 
very  well  with  the  results  of  visual  interactive  analysis  using  the  EP  system’s 
graphic  options. 


A.F.  Kushnir,  Int.  Inst,  of  Earthquake  Predic.  Theory,  Moscow,  USSR 

J.  Fyen 
T.  Kvaerna 
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Frequency  UamI 
(Hz) 


quency  Band 
(Hz) 

Ordor 

C 

of  Multidi 

8 

meii.sioiial  .All  Model 

12 

0.2  -  5 

16.7 

17.9 

18.8 

(21.9) 

(22.4) 

0.2  10 

i;i.2 

- 

16. .7 

(Ui.S) 

(20.2) 

0.2  20 

14.0 

(17.0) 

Frequency  Band 
(Hz) 

Uegularizator  Val 
10-“  I0-* 

0.2  ■ 

-  5 

18.9 

18.8 

18.7 

(22..0 

)  (22.4) 

(22.3) 

0.2  - 

10 

1G..5 

(20.2) 

0.2 

20 

1  1.0 

14.2 

(17..5) 

(17.6) 

Table  7.4.1.  .Adaptive  optimal  )»roup  filtering  S.Ml  gain  relative  to  beam  forming 
for  different  values  of  the  GKFAD.M’T  procedure  main  |)arninelers.  Gain  values 
are  given  in  dB,  values  in  brackets  are  .AOGF  S.NU  gains  relative  to  average 
single  array  channel.  The  results  are  ba.sed  upon  using  120  sec  NORESS  noise 
recordings  shown  in  Fig.  7, .1.2.  The  results  of  the  ujiper  table  are  obtained  with 
a  regiilarizator  '  alue  of  i0"'‘.  whereas  an  .AH  model  order  of  12  has  been  used 

I  for  obtaining  the  results  in  the  lower  table. 

I 

i 

I 
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Event 

Array, 

Phase 

Gain  Relative 

Gain  Relative 

Frequency 

Origin 

Distance 

(AR  Model 

to  Beam 

to  Single  Chan. 

Band 

Time 

(km) 

Order) 

(dB) 

(dB) 

(Hz) 

298:17.51.50 

ARCESS 

508.4 

P  (12) 

16.2 

18.8 

0.2-5 

282:12.04.13 

FINESA 

P  (12) 

21.7 

23.3 

0.2-5 

287.9 

P  (12) 

19.3 

21.0 

0.2-10 

S(12) 

19.6 

21.6 

0.2-5 

S  (12) 

17.2 

19.3 

0.2-10 

Lg(12) 

17.3 

19.4 

0.2-10 

294:09.13.00 

FINESA 

P  (12) 

23.6 

25.7 

0.2-5 

772.3 

S  (12) 

21.0 

23.7 

0.2-5 

Lg(12) 

20.3 

23.1 

0.2-5 

NORESS 

P  (12) 

18.8 

22.9 

0.2-5 

1302.7 

S(12) 

16.6 

21.1 

0.2-5 

Lg(12) 

16.5 

21.3 

0.2-5 

28:09.38.09 

NORESS 

P(6) 

16.0 

19.6 

0.2-5 

1219.0 

P(6) 

12.3 

16.1 

0.2-10 

P(6) 

10.1 

13.8 

0.2  20 

P(8) 

10.9 

14.6 

0.2  20 

FINESA 

P  (12) 

23.1 

24.3 

0.2-5 

1771.0 


Table  7.4.2.  S^!R  gains  of  adaptive  optimal  group  filtering  relative  to 
beamforming  based  on  processing  of  local  event  phases. 
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Event 

Array 

Phase 

Gain  Relative 

Gain  Relative 

Frequency 

Origin 

Distance 

(AR  Model 

to  Beam 

to  Single  Chan. 

Band 

Time 

(km) 

Order) 

(dB) 

(dB) 

(Hz) 

292:12.31.45 

FINESA 

P(6) 

14.6 

17.6 

0.2-1 

164.4 

P(12) 

16..'! 

19.4 

0.2-10 

ARCESS 

390.4 

P  (12) 

16.4 

20.8 

0.2-10 

294:19.32.03 

FINESA 

P  (12) 

24.5 

26.0 

0.2-10 

259.0 

NORESS 

P(12) 

17.9 

21.8 

0.2-5 

564.6 

S  (12) 

16.0 

19.8 

0.2-5 

^g(12) 

1.5.7 

20.0 

0.2-5 

Table  7.4.2.  cont. 
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Output  from  GRFADAPT:  (P  -  phase.  0.2-5  Hz) 

Start  tine;  1990-282:12.02.50.0  Seconds:  120.00  (Data  t.ime  interv.al 
Input  rarmetere: 

Filter  type:-l  (low  pass ( Ip) =~ 1 , band  pasB( bp ) =0 , high  pase(hp)-l) 

FCH  5.25  ( lp,hp-f iltere  cut  frequency,  hz) 

FLH  :  20.00  (up-filter  low  cut  frequency,  hz) 

FHH  :  20.00  (up-filter-high  cut  frequency,  hz) 

IRL  :  1.5  (one-eide  length  of  the  filter  impulse  response) 

ALPHA  0.00010  (decay  factor  of  the  filter  frequency  response) 

AZIMUTH  ;  153.80  (azimuth  of  the  seismic  phase  to  be  extracted' 

VELOCITY  :  7,80  (apparent  velocity  of  this  seismic  phase) 

DAPB  10  (order  of  a  beam  noise  AR-model) 

LCRC  :  6  (number  of  autocovariance  roatrises) 

REG  0 . 000001 ( regularizator  value) 

DARGRF  6  (order  of  the  data  multidimensional  AR-model) 

DARARF  10  (auxiliary  parameter) 

DMARF  :  20  (auxiliary  parameter) 

K  :  4  (resampling  factor) 

Output  parameters: 

NP  1193  (number  of  the  data  samples  being  processed) 


AVCHPOW 
BMEAN 
BPOW 
lERLBLD 
lERMLD 
lERGLD 
GRTUN  (M/P) 
GRTW  (M/P) 
GRAUH  (M/P) 
GRMUN  (M/P) 
Gain  (avch) 
Gain  (beam) ; 


0. 1380 
0.0000 
0.0946 

s  } 

0.0003 

0.0098 

0.0000 

0.0006 

214.33 

146.95 


(averaged  diepereion  of  the  channel  traces) 

(mean  value  of  the  bean  trace,  m.) 

(dispersion  of  the  beam  trace,  ) 

(computation  errors,  error=l.  otherwisc=0) 

0.0006  (m  and  of  the  AOGF  output  trace) 

1.0128  (n  and  of  the  whitened  output  trace) 

0.0007  (m  and  a*  of  the  auxiliary  output  trace) 

0.0006  (m  and  a  of  the  auxiliary  output  trace) 

23.311  (AOGF  SNB  gain  relative  averaged  channel) 
21.672  (AOGF  SHE  gain  relative  beam,  times,  db  ) 


(Description  of  the  array  data  being  processed) 


Channel 

X(E-W) 

X(N-S) 

ELEV 

TDEL 

CHMEAN 

CHPOW 

FIN  AO  sz 

-180.0 

-85.0 

138.0 

0.211 

-0.0070 

0.1287 

FIN  Al“sz 

0.0 

0.0 

138.0 

0.211 

-0.0181 

0.1421 

FIN  A2”sz 

-308'.  0 

-353.0 

162.0 

0.262 

-0.0085 

0.1308 

FIN_Bl”sz 

275.0 

-37.0 

165.0 

0.252 

-0.0074 

0.1264 

FIN_B2_sz 

121.0 

-599.0 

159,0 

0.371 

-0.0053 

0.1280 

FIN_B3_SZ 

-474.0 

-555.0 

176.0 

0.292 

0.0071 

0.1258 

FIN  B4  SZ 

-764.0 

-85.0 

158.0 

0.143 

-0.0167 

0.1171 

FrN_B5_sz 

-436.0 

257.0 

143.0 

0.098 

0.0206 

0.0913 

FIN  B6  sz 

83.0 

277.0 

147.0 

0.153 

-0.0255 

0.1501 

FIN  Cl  SZ 

-1064.0 

-657.0 

158.0 

0.248 

-0.0223 

0.1398 

FIN  C2  sz 

-1110.0 

226.0 

138.0 

0.027 

-0.0460 

0.1368 

FIN  C3  sz 

-162.0 

788.0 

138.0 

0.000 

-0.0100 

0.1507 

FIN  C4  sz 

629.0 

420.0 

138.0 

0.182 

-0.0373 

0.1256 

FIN  C5  sz 

653.0 

-518.0 

138.0 

0.413 

-0.0160 

0.1374 

FIN  C6  sz 

-185.0 

-1108.0 

138.0 

0.460 

-0.0329 

0.2375 

Fig.  7.4.2.  NORESS  noise  recordings  used  for  studying  the  influence  of  GR- 
FADAPT  parameter  values  on  SNR  gain. 
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Fig.  7.4.3.  Extraction  of  the  main  seismic  phases  from  an  event  recorded  at 
FINE.SA  (1990:294:19.,12.03,  distance  2-59.0  km),  (a)  Input  and  main  output 
traces  of  the  GRFFILT  procedure  adjusted  for  Lg-pha,se  extraction;  (b),  (c), 
(d)  main  output  traces  of  the  GRFFILT  procedures  adjusted  for  P,  S  and 
Lg-phase  extraction.  0.2  .5  Hz  band  processing. 
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rig.  r.4.4.  Output  traces  of  GRFFILT  (0.2-5  }lz  band)  procedure  adjusted 
for  P-phase  extraction  from  ARCESS  recordings  of  a  small  regional  event 
(1990:298:17.51..50,  distance  508.4  km).  Note  the  difference  between  P-phase 
waveforms  of  the  AOGF,  Whitened  beam  and  Whitened  AOGF  output  traces 
which  implies  a  shift  of  onset  time  estimates. 


Fig.  7.4.5.  Extraction  of  the  S-phase  from  NORESS  recordings  of  a  regional 
event  (1990. 294.'09, 13. 00,  distance  1302./  kni);  0.2  -5  Hz  band  processing. 
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1  9  90-28.2;  1  2.04.50.000  OGF_th_un 

Fig.  7.4.7.  S  and  I,g  pha.se  rtelerlion  ami  ousel  lime  estimation  using 
AOGF  output  trace  of  the  local  event  described  in  Fig.  7.4.6. 


1  99f.'-292;  1  2.32.0  1  .950  OGF 


Fig.  7.4.8.  I’-pha-sc  lion  and  onset  Uinoostimation  using  ihc  .\OGK  out 
put  irarp  from  FINFS.'V  rprordings  of  a  small  local  event  ( 1990:'2!)2: 12.3 1  .■l.'i.l), 
di.stance  166.4  km). 
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Fig.  7.4.11.  Extraction  of  the  Lg-pha.se  from  the  recording  described  in  Fig. 
7.4.5  and  onset  time  estimation  of  this  phase  using  the  beam  and  the  AOGF 
output  traces. 
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7.5  A  2-dimensional  finite  difference  approach  to  modeling  seis¬ 
mic  wave  propagation  in  the  crust 


Introduction 

It  is  well  known  that  the  direet,  disrreto  solution  of  the  elastic  wave  equation 
constitutes  an  excellent  platfornt  for  synthetic  seisinograin  analysis  as  all  propa¬ 
gation  effects  are  included  in  the  solution  (e.g.,  see  Mooney,  lOSt.')).  ..\  practical 
realization  of  this  approach  has  been  probleinalic  tinlil  recently  due  to  limitations 
imposed  by  currently  available  computers.  This  being  said,  we  will  report  below 
on  '2-dimensional  (2D)  finite  difference  seismogram  synthetic  experiments  which 
have  been  achieved  through  cooperative  efforts  with  scientists  at  IH.M  liergen 
Scientific  Centre  (Bergen,  Norway). 

Elastic  wave  modeling  formulation 

The  basic  equal  ions  governing  wave  propagation  in  a  continuous  elastic  medium 
are  the  momentum  conservation  and  the  stress- strtiin  relation.  Following  ,\(hen 
bach  ( lOTti).  in  the  velocity-stress  formulation.  the-,e  are  given  by 


</  <t 

‘Tt 


jj=  1. 


t)  i) 

I'r  +  7/1  — 
axf  oxj 


jj^  1.,,. 


0  (  0  \  ,  ,  ,  , 

=  +  —  e,j.  . /,  ,/f  C) 

where  Einstein’s  suntmation  convention  is  used.  J  is  the  dimensionality  of  the 
problem,  p  is  density,  and  A  and  p  are  Fame's  parameters.  //  are  body  forces 
and  Vj  and  Ojf  are  velocities  and  stresses,  respectively. 

,V u mrriral  iliscrctizat ion 

Spatial  partial  differentiation  is  achieved  through  cost  optimized,  dispersion 
bounded,  high-order  finite  difference  operators  on  a  staggered  grid.  For  titne 
stepping  a  leap  frog  technique  is  ii.sed.  The  discretization  of  l!ie  elastodynamic 
equations  with  twej  staggered  numerical  space  differentiators,  rr*,  applied  as  in 
I.evander  ( 1988)  to  stresses  and  particle  velocities  leails  to-. 

J 

p+{Vj+(t  -F  A(/2)  -  V)+(f  -  A</2)}  =  A<{E/(<)  +  +  ^  I 

IgJ 

.9,, (I  -F  A/)  -  %(f)  =  AAI  ^  8;v*{l  +  Al/2)  -F  2,,A/fi;  Vs+(l  +  Al/2),  ;, /  =  1 ,  - 
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with 


V;.+  (<)  =  vj{x  +  hi/7,  t),  F;(t)  =  fi{x  +  hi/2,1), 

~  '’'//(’t  +  hj/2  +  h//2,  J), 

pf  =  p{x  +  hi/2),  X  =  X(x),  /i  =  /i(x)  and  /i/,+  =  ,i(x +  hj/2  +  h,/2). 


q(x  +  {(-  l)hj)  -  q(x  -  (hi)  ^  Oq 


(x-h,/2). 


Here  h(j)  is  the  unit  vertor  in  the  jih  dirertion.  A.  p  and  S,j  are  deruied  at  the 
nodes  of  the  Cartesian  mesh,  p'* ,  V*  and  F*  and  defined  at  tlie  links  rotineetins’ 
the  nodes  and  and  an.  defined  at  the  centers  of  the  '■plaqiiette.s''.  rr*  are 
niitneriral  differentiators  of  rooffirionts  q  is  here  velocity  or  stress  and  /.* 

is  the  length  of  the  operator.  For  the  numerical  dispersion  relations,  the  staliilitv 
limit  and  bandwidth  introduced  By  the  discretization,  the  reader  is  referred  to 
Sguazzero  et  a/ (1990). 

Absorbing  and  free  sttrfact  boundary  rnnditiotis 

By  necessity,  the  numerical  modeling  limits  the  medium,  and  to  reduce  ar 
tificial  reflections  from  the  numerical  boundaries,  the  velocities  and  stresses  are 
multiplied  by  c.xponentially  decreasing  terms  near  the  edges.  For  this  procedure 
to  be  efficient,  relatively  large  models  are  required,  that  is,  relatively  large  spatial 
distances  to  the  wedges  and  this  in  ■'ID  modeling  would  be  computationally  very 
demanding.  In  the  latter  ra.se  we  have  experimented  with  boiindarv  operators 
recently  introduced  by  Higdon  ( 1990.1991 ),  which  at  r  =  0  read  like 


which  will  ab.sorb  perfectly  a  plane  wave  t ravelling  towards  the  boundary  at  angle 

lO.a 


Oj  and  speed  Cj.  ni  is  the  order  of  the  operator.  Similar  operators  are  used  on 
the  other  boundaries.  The  condition  for  this  method  to  be  ii.sefnl  is  ibal  the 
number  of  time  steps  is  small  enough  not  to  exceed  a  certain  limit,  after  which 
the  method  will  appear  unstable. 

On  the  to|>  fr<'e  surface,  we  use  the  vanishing  stre.ss  conditions  for  a  free 
hound  ary 

n-T  =  0  (1) 

Mere-  nisi  he  outward  normal  unit  vector  on  the  surface  and  T  is  t  he  si  ress  tensor. 
To  get  computationally  t  ractahle  conditions,  we  assume  the  free  top  sin  face  to  be 
locally  plane.  I  hen  n  =:  k,  where  k  is  the  unit  vector  in  the  vertical  c-directioii. 
s  an<l  >J  are  horizontal  coorclinatos.  (11  then  leads  to 

'Trj-  =  =  <T,.  =  0  (.'■) 

To  increase  the  generality,  one  may  assume  a  topograithic  reKK-f  as  the  free 
surface.  Ily  relaxing  the  re<|iiiremenf  of  the  surface  being  locally  plane,  one 
as.siiiiies  a  given  slope  locally  in  eacb  spatial  direction,  I'lie  resulting  conditions 
on  the  Stresses  become  more  complex,  though  Iraciabie,  as  demonstrated  bv  .lib 
<1  al  ( IbsX). 

;\l  present  rvi'  have  not  incorporated  the  "topography"  free  surface  in  our 
software. 

f 'ru.stn/  irnri  iiroimfinliim  Jl)  Jiiiili  diffi  n  iii  t  .■ttpil/n  lir.- 

The  task  of  "adapting"  the  21)  I’D  software  for  liandling  of  sri.smologli  al 
problems  has  Ikmui  rather  lime  (..n.stiming.  Hence,  only  recently  have  we  been 
able  to  produte  seismic  synthetics  for  crustal  wave  proiiagalion.  W'e  can  al.so 
handle  :(I)  cases,  but  their  seismological  relevance  at  prc.sent  is  limited.  .Anyway, 
in  llie  following  we  will  pre.sent  .some  examples  of  .synthetic  seismogram.s. 

Mfnlf  I  (U  script ititi  iinil  ilafn  tiniilysis 

llasicallv  we  use  a  homogenxuis  crust  of  thickness  .10  km  and  I’,,,'  .-  li..' 
km/sec,  which  l>«*sides  serves  as  a  referem'e  imxlel.  The  options  for  perturbing 
this  model  comprise  multil.iyering,  pit  "ewi.se  linear  velocitv  gradients,  large  sc  ale 
discontinnilies  like  Moho  l)iimp(s),  b\it  so  far  no  randomized  scatter  inclnsions. 
A  schematic  model  illustration  is  shown  in  Fig.  7..">,1,  .Although  the  source  (point 
or  tine  source)  crmhl  b<’  at  any  depth,  the  .^.ecisors  are  always  on  the  free  surface. 
.Any  sensor  configuration  could  be  used,  althoiigb  our  performance  is  for  a  10 
element  line  array  with  0.1  km  sen.sor  interspafing,  which  i.s  rrznvenient  for  velocity 
decomposition  of  the  synthetics.  Occasionally  we  use  a  sensor  spai  ing  of  .A  km  in 
order  to  visualize  the  distance  variability  in  the  records. 

.An  objection  against  21)  solutions  of  the  ela.stic  wave  ecpiation  is  I  hat  all  prop 
agalion  elfects  are  included  and  hence  it  would  be  difficult  to  isolate  the  response 


lOf. 


of  a  specific  body  within  the  synthetic  wavctrain.  To  overcome  this  kind  of  prob 
lems,  we  would  process  the  synthetic  records  in  a  manner  similar  to  that  used 
for  real  recordings.  Principal  techniques  used  are  frequency  wavenumber  (f-k), 
semblance  and  3-component  polarization  analysis  (e.g.,  see  Husebye  and  Ruud, 
1989).  Occasionally  we  would  make  comparisons  with  “ray  tracing”  synthetics 
for  which  more  specific  contribution  effects  arc  specified  a  priori. 

licstilts 

Examples  of  crustal  synthetics  using  the  procedure  outlined  above  are  shown 
in  Eigs.  7. .'1.2,  “..u.S  and  7. ,5.4.  The  following  comments  apply. 

Figs.  7.5.2  and  7.5.3:  Buinp  on  Moho  ranges  160  km  and  210  km 

In  Figs.  7.5.21>  and  7.5.3b  the  homogeneous  cases  are  shown,  while  the  bump 
ra-ses  are  shown  in  Figs.  7.5.2a  and  7.5.3a,  respectively.  A  comparison  here  gives 
that  the  Moho  bump  does  not  strongly  change  the  records,  which  is  also  rather 
obvious  from  a  corresponding  compari.son  of  the  .semblance  plots  in  Figs.  7.5.2c.d 
and  Fig.  7.5..'tc.(l.  respectively.  The  dominant  features  in  the  synthetics  appear 
to  bo  crustal  reverberations  (PmP),  which  are  particularly  abundant  since  the 
signal  source  was  put  at  a  depth  of  10  km. 

Fig.  7.5.4:  Bump  on  Moho  -  ranges  100  200  km 

In  this  case,  we  used  a  linear  velocity  gradient  in  the  crnsl  and  besides  used 
a  sensor  spacing  of  5  km  in  order  to  visualize  distance-dependent  changes  in  the 
records.  As  observed,  the  Pg-pha,se  dominates  the  first  part  of  the  records,  then 
comes  the  corresponding  S-phasos  and  finally  multimode  Rayleigh-type  of  waves, 
V»iih  t!.c  iiiiicl.  I<.:ger  5c:'.f,or  spacing  the  scmblai'ce  resolution  is  very  high,  as 
illustrated  in  Fig.  7.5.4c,d.  We  have  also  tested  the  signal  polarity,  which  further 
adds  weight  to  the  realism  of  the  2D  F'D  synthetics  displayed. 

Ui,<!cu!iriou  and  future  work 

riie  synthetics  generated  seemingly  include  nil  major  phases,  while  iii  com 
parison  to  real  records  the  body  wave  coda  is  weak  to  nonexistent.  .\s  demon 
strated,  long  wavelength  heterogeneities  like  a  bump  on  Moho  do  not  rontrib\ite 
much  in  this  respect.  This  in  turn  implies  that  the  cumulative  propagation  effect 
of  randomly  distributed  scatterers  are  likely  to  be  of  importance. 

.■\  specific  advantage  with  our  technique  for  21)  synthetic  seismogram  calcu¬ 
lations  is  flexibility  in  choosing  model  parameters.  In  our  future  work,  some  sort 
of  a  reference  crustal  model  would  be  established.  'I’heu  we  would  systematically 
change  the  velocity  structure  both  above  and  below  Moho,  .Scatterers  would  be 
introduced  at  various  parts  of  the  travel  path,  and  their  effect  would  be  visualized 
partly  but  taking  tlie  difference  between  “homogeneous”  and  “inhomogeneoris  ' 
synthetics.  Finally,  we  would  naturally  compute  .synthetics  on  the  basis  of  crustal 
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results  presented  in  Section  7,6. 
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30  km 


2D  ELASTIC  MODELLING  -  OSLO  RIFT 


PROCESSING:  F-K;  SEMBLANCE  AND  3C  ANALYSIS 


OBSERVATION:  HORESS  ARRAY  RECORDINGS 


Fig.  7.5.1.  Simple  ono-layorcd  crusfiil  model  used  initially  for  compiitiiig  svn- 
ihefir  spisniogra:iis  l)a,scd  on  finiU'  diirrn'nce  .soliitioii.s  in  2.dim('nsional  (21))  of 
the  olasfir  wave  rqitalions.  'I'hp  poini  sourrp  ).•;  loralod  at  a  deptii  of  10  km: 
Prantal  and  suh-Moho  vploriljps  arp  fl,.a  k:n/spp  anil  8,2  km/sor,  rp.spprt i vt'lv. 
('ortpsponding  dpnsity  values  are  2.8.''.  kg/nv’  and  kg/m‘.  The  Moho 
is  .0  wide  and  2  or  1  kiii  high.  llori/o:ilal  distance  from  source  to  nearest  edge 
of  Moho  hnm|)  is  I  0  km.  So  far  scatter  inclusions  with  contrasts  in  velocity  and 
density  of  the  order  of  2  per  cent  have  not  been  inclmled. 
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Fig.  7.5.6b.  Sf'fnbianro  vriority  (VKSIV\(JF\M  )  analysis  of  tin*  syntl.-'firs  dis 
piaynd  in  Fig.  7. 5. Ha.  J  Ik*  first  part  of  the  syiiUnairs  is  doniinat«‘d  by  rnistal 
Tf’vnrbfrations  (|>bas<*  voloritjos  abovo  S.O  kin/s<*r  hardly  sfHMi ).  Thr  phas»' 
around  .io  snc  pr^ncding  ih*''  S-wavo  at  annin(!  .'tb  see  is  a  coniinonly  ol)snrv<'d 
h'atnrn.  Dispersive  Ravl<’ip,h  waves  are  also  synthesi/«*d.  I'inn'  referenre  is  lied 
tf)  s#’nsor  Ot), 


7.6  Crustal  tliicknesses  in  Fennoscaiidia  —  An  overview 

Barkfjmund 

Crustal  studios  Itoraino  popular  anioiii»  soisiiiolottists  iu  tlio  I'onnosraiidiiia- 
viau  (  oiiutrios  somo  tliroo  docailos  ago,  and  .still  romain  so.  I'lio  nuiiioroiis  .soisniic 
surveys  coiidiictod  within  this  roffion  are  aiinod  at  uiapiiiiii'  crustal  struct iiros  iti 
ovor-incroasiut!;  detail.  We  havr‘  reviewed  the  knowledge  a< cuiiiulated  from  these 
studies  and  itiade  a  new  rrustal  thickness  niap  with  contour  intervals  of  '1  ktn 
for  Feiiuoscandinavitt.  In  sottie  areas  the  sediment  thicknesses  e.xceinl  10  km.  so 
it  is  iiitportant  to  dilferetitiate  hetween  .\Ioho  deptli.s  and  the  crust  (  ry.slalliiie 
thicktiesses.  Hence  for  thesoulh<‘rn  parts  of  I'etitioscatidinavia.  notahly  Detitiiark 
atid  adjacetit  .seas,  ati  additional  tiiap  of  crystallitie  c.ustal  Ihickne.sses  was  made, 
ilelow.  we  will  presmii  the  major  re.siills  from  thi.s  crustal  .study,  while  for  details 
we  refer  to  a  forthcotiiiuK  pap<‘r  hy  Kinck  il  (tl{  1901 ). 

(l(ol<xjirnl  hiiiiu  uxifk 

(leonraphically.  t he  I'luittoscandinavian  pan  of  t he  Halt ic  Shield  comprises  the 
Kola  I’eninsiila  (ititludintt  the  While  Sea),  i'itdaml.  the  ScatidinaviaTi  I’eiiitisiila. 
Dentnark  .and  adjacent  seas  (SkaKerrak.  Katleftai.  the  H.dtic  Sea  atid  parts  of 
tlie  Haretits  SeH)in  Keolo^ical  terms,  this  area  (l  it;.  7.0.1)  e.xhibits  a  variety  of 
difh'rent  tectonic  provinces.  raiiKrtij^  in  from  Anluxtii  to  J’ermian.  1  he  more 
rt'cenf  openini’  of  the  .North  .\llanlic.  commencing,  .some  .'.d  Ma  ago,  affected  only 
peripheral  parts  of  the  shitdil.  that  i^,  the  oiastal  areas  i,d  wi'Steru  and  northern 
.Norway. 

f'rttshil  pri/filinij  Moho  ili  plh  niiipptiin 

The  principal  aim.s  of  crustal  prolirm,i>  surveys  are  crustal  thicknesses  and 
velocity. dejtth  (list rihttl ioiis  ah<jv<‘ and  ludow  Moho.  1  he  former  paratneter  seenis 
well  constrained  In  view  of  small  dilferetices  of  the  order  of  2  d  km  either  hetween 
iidersecting  prolilitig  lines  or  between  rellectioti  and  refraction  lint's.  liegartliiiR 
velocity  depth  distributions  the  reflection  prolilinR  tiata  havt'  poor  resolution. 
The  refraction  profilint’  data  have  relatively  R..;;..!  re.solution  although  the  inver¬ 
sion  schemes  in  general  use  ilo  not  give  uni<|ue  results.  It  suHices  here  to  tnei. 
that  tlilferetit  groups  of  researchers  using  the  same  set  of  observatiotial  data  sel 
dom  produce  the  same  velocity-depth  distribution.  I'lie  inherent  problem  here  is 
that  the  identilication  and  picking  of  secondttry  phase  arrivtils  often  are  iliHicult 
and  hettcf'  the  final  solulion  is  not  well  const rai tier!.  Kinematic  rav  tracing  is  tiot 
too  helpfid  in  this  respect  since  amplitude  inlormation  ami  scattering  contribu 
tions  are  tmrstly  igtiored.  Also,  there  appears  to  he  a  sigiiilicant  improvement 
in  the  published  profiling  results  frrmt  the  mid-.seventie.s  and  onwards,  reflecting 
l.etter  recording  instrumentation  (digital),  denser  sampling  and  the  \i,se  of  more 
soiihislicated  analysis  and  inlerpretational  methods.  I'liese  brief  comments  on  the 


reliability  of  seismic  reflection  and  refraction  profiling  results  should  be  kepi  in 
mind  when  judging  the  major  outcomes  of  our  study  (  Kinrk  (t  <il.  1!MM  ).  namely, 
a  Moho  depth  map  for  lennoscandiiiavia,  tlii<;kn<‘sses  of  the  crystalline  cnisl  in 
the  southern  parts  ()f  the  region  (Denmark  and  adjacent  sea)  plus  a  tahiilation 
of  P-velocity  depth  distributions  for  selected  profiles. 

t'cnnosrundniavian  st  isinir  stinu  ys 

We  liave  carefully  screene<i  the  available  literature  for  |)rofiling  survev.s  within 
Fennoscandinavia,  and  the  onicoim'  of  the.se  efforts  is  tabulated  in  1’ahle  7  O.l 
ami  displaye<i  in  Fig.  7. <>.2.  .\'<ne  that  data  from  soim*  of  tin*  profiling  lim's  have 
been  reanalyzed  and  n’itilerpreicil  and  witli  frvv  exceptions  we  only  refer  to  the 
latest  publication  in  tis  regar<i.  .\  final  remark  here  is  that  imleed  nmdi  c-ffort 
lias  been  invested  in  the  crustal  mapping  of  Fennoscamiinavia. 

lit  Mitliit  (It  fttli  (iinl  rrystulinu  t-ntsltd  lliit  knt  ■''>  /or  I  t  muliiuirid 

In  Fig.  7.<>.:t  the  Moho  <le[)ih  map  is  shown  ami  in  Fig.  7.G.  1  the  crystalline 
crustal  thickness  map  (limile<i  ((»  Denmark  ami  adjacent  s<‘as)  aresjiown.  .A  map 
similar  to  that  in  Fig.  7.<).l  was  atiejn|)le<l  (onstrucied  f‘>r  the  Kola  l*i’niii*.ula 
area,  the  White  Sea  ami  Hie  western  Mar^uits  se.».  bui  at  pres«uit  liere  .,r«’  mu 
euoiigli  data  available  for  smh  a  task.  .Anyway,  the  .Moho  map  in  Fig.  7.fi,;i 
Is  rather  detailed,  in  particular  in  the  areas  olFhore  N«»rway.  as  wr  havr  lu-eii 
alile  to  incorporate  jrci-nt  results  from  marine  seisnii<  retloctiou  snrvevs.  Th«‘ 
crustal  thickening  is  in  general  perp<  iHlicnlar  to  the  coastal  areas  of  soiiih«>rn 
ami  western  Norway,  ami  the  Kola  INuiinsula.  hut  h's>  so  for  th<‘  iiiterpiate  Ualtic 
Sea.  Ill  general,  the  <»ldest  parts  «)r  ilie  lialtic  Sliield  (the  maj<»r  parts  of  th.' 
Fennoscandinaviati  region )  exhibit  tin' greatest  crustal  l hi<  km'v-,<‘s.  This  niav  he 
expres.s<’ii  in  liie  following  form; 

If  I7.:i/ofl(  7')  -  10.2  (  1  i 

where  (I  in  km  is  .Moho  depth  afi<l  7'  l.s  time  in  Ma. 

riio  sediment  t  hicknosse.s  in  i  he  liasin  areas  offshore  Nhirway  are  often  for  mi  da  IF' 
with  corresponding  thicknesses  of  the  crystalline  crust  »>f  tlie  order  of  1."  I'o  km. 
i'liere  is  no  obvious  correlation  with  age  ImU\v«sui  the  crustal  F  veio<iiy  iloptli 
distriliulion.  although  wlietlmr  we  have  piec«'wi.sr  negative,  zero  or  pi>.si!ive  ve 
locily  grarlients  Is  likely  to  afre<t  prohnimliy  seismic  wave  propagation  in  tlte 
crust.  Hegao'F.g  lateral  l‘u  ami  .sn  vehicify  vaii.ilions  within  lVnm>s<  amiin.o  ia , 
this  jirohlein  lias  ln'<‘n  slmii^al  1)V  lotno^rapliic  lerlitii<|ues  using  lo<al  seisniulogi 
<al  iuilletin  data  (e.g..  .-ee  Ifannister  tl  al,  1001  ].  I  heir  ma  joi  li;.  !;•  "<  are  Hiat 
proiioiinced  low  v<‘hMily  ar^-as  are  associated  with  tlie  ('ahdoni<le  im'unt.tins  nf 
wesO'rii  Norway  ami  the  rift  and  hasin  areas  olfsliore  Norwa\ .  The  ( «'nt ral  part s  cf 
the  shield  are  rather  liomogetieoiis  in  tliis  respect.  .\  corresp<uiding  t(iim'gr.iphi< 
'•Indy  nf  crustal  ve|oeity  variations  f|*g  ami  Sg  pha.ses)  w.t.s  not  attempted  '.ime 
t  lu’  I’g/Sg  ray  paths  t  annot  unnjUely  he  determined. 
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( to  many  oUht  rontiiMMital  tin*  n'snll.s  display**!]  in  1  i^. 

r.^i.d  arf<l  7Ji.  1  ajv  ind^vd  vnrv  d4‘lail4*»l.  On  tin*  other  hand,  structural  <h'tails 
aj!’  very  poorly  resolved,  whicli  in  turn  r(*tl**cts  the  liata  at  hand;  niaitdy,  r** 
tra(  lion  proliline  results.  .Mtiiou^h  these  result-'  are  no!  a«ie<piaU'  for  rest  rai 
hy  jiot  !lese^  on  f  In*  f e*  r(»td<’  evolnt  ioti  af  I  ennoscatidinav  ia .  we  do  t  hink  l  hat  I  h«'S<' 
result  V  may  la*  iiist  ruiti<*nt  a!  in  juovidin^  a  Im*i  f  er  und!Tsfandifi^,  of  seismic  recorris 
.ii  lof  al  amt  re|>ional  <lislam«*s  tljr!!!!!*!!  synthetic  '<*isniou,raiu  analysis.  In  tlii'- 
!e>pe«i  We  CiHismler  tin*  '20  litiiteil  <liirer<'fK<‘  lechniqj;'  presenter]  jji  Si'ction  7.' 
•  o  1'*'  most  '•uiiahh*  sin(‘<‘  \v»'  could  inrorptrrai!*  a  filling,  Nt!>h!)  loeriln'r  with  anv 
ktmi  of  veha  ity  «:ia!li!'nt  al'ove  aml/or  helow  Moho. 

A  liiial  retnai  k  is  that  t In*  Moho  <lept h  variai it)n  app!’ars  to  hava-  a  <ount«  rparl 
iti  the  sjvatial  dist ribut i«>i,  of  earlh<piakes  within  this  region.  As  is  W!'li  known. 
i!ie  '■ei-nii' ity  i'-  by  far  laiuest  in  t!'*-  <-oasial  ar'*a‘- of  Norwav.  whert*  tlM*  (rust  is 
e\-cep;ionally  il/in.  r»jrtlienm»r*-.  all  iIm*  larm"'’!  «*.ti t  h!piak‘’s,  incliiditij*  the  his 
liuital  <»nes.  have  iak‘‘n  pla(!*  in  ar!*as  wln'ie  ili!*  <  rn>t  is  thin.  In  olli‘‘r  words, 
'tress  iuctimnlal  itiiis  within  lennos!  amlin.o  ia  appe.;r  'o  be  insirlfif  ient  for  crack 
me  or  cansine  fna  }<*r  »*arf hqtiakf's  in  ar‘‘as  with  lhi<  k  <  rust  ( //  >  U)  km  ),  which 
imt nr.dl V  is  st  loneer  than  the  thin  <  rust  in  tin*  <  <uist al  areas.  Nat  nrally.  t luT!’  ar<' 
imiir.  ar!-a'>.  includinu  l)!-nmark.  with  thin  <  iust  btH  seismically  !jniesci’nt. 
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7.7  Initial  development  of  generic  relations  for  regional  threshold 
monitoring 

Introduction 

In  earlier  reports  (Kvatriia  and  Ringdal,  1030a.  Kvceiiia  and  Ringdai,  iOOOb) 
we  have  demonstrated  applications  of  the  threshold  monitoring  (TM)  technique 
to  regions  of  limited  areal  extent  like  mines  and  nuclear  test  sites.  This  method 
has  proven  to  provide  a  simple  and  very  effective  tool  in  day-to-day  monitoring 
of  areas  of  particular  interest.  One  of  the  ba.sic  underlying  assumptions  has 
been  that  each  target  region  should  be  defined  such  that  all  events  within  the 
region  show  similar  propagation  characteristics.  This  has  enabled  us  to  get  the 
necessary  magnitude  calibration  factors  from  processing  previous  events  with 
"known"  magnitude,  using  the  relation 

huj  =  'i'v  -  log(.S',,j)  (i  =  I.-  •.  A  ;  y  =  1 .  •  •  •. /,)  (1) 

whore  b,  j  is  our  estimate  of  the  magnitmle  correction  factor  for  phase  i  and  event 
j.  rhj  is  the  estimate  of  the  magnitude  for  event  j  ( l>a.se(l  on  inde|>endent  networks 
or  knowledge  about  the  exjilosive  charge)  and  S,^j  is  our  I'stimate  of  the  .signal 
level  at  the  predicted  arrival  time  of  phase  i  forevent  j.  A’  is  the  number  of  phases 
considered  (there  might  be  .several  stations  and  .several  pha.ses  per  station),  and 
I.  is  the  number  of  events. 

file  magnitude  correction  factor  to  bo  insed  for  phase  i  is  then  given  by 

6,  =  a;  <  i>,  j  >  CA) 

where  A'  denotes  statistical  expectation.  Parameters  like  window  lengths  for 
signal  level  estimation,  travel-times  of  the  different  phases.  fre(|uency  filters  and 
steering  delays  for  array  beamforming  are  taken  from  processing  of  the  calibration 
events. 

Extension  of  the  TM  method  to  regions  where  no  calibration  events  are  avail¬ 
able,  requires  that  we  have  generic  formulas  for  all  variables  describing  the  pro 
cessing.  Such  relations  will  make  it  possible  to  monitor  new  and  larger  geograph¬ 
ical  regions,  and  will  in  addition  enable  us  to  get  a  more  thorough  understanding 
on  how  events  originating  in  one  region  influence  the  threshold  in  other  regions. 
.Applying  such  generic  relations  will  of  course  involve  a  tradeoff  where  a  wider 
geographical  coverage  is  achieved  at  some  expense  with  regard  to  optimized  mon¬ 
itoring  of  limited  target  area.  Thus  it  should  be  seen  as  a  supplement,  and  not  a 
replacement  of,  the  target-specific  threshold  monitoring. 

In  the  following  we  present  results  from  a  preliminary  study  on  methods 
for  obtaining  such  generic  relations,  with  special  application  to  the  regional 
Eennoscandian  array  network. 


Phases  to  consider  and  their  travel-times 

A  standard  method  of  estimating  the  magnitude  of  local  and  regional  events 
is  based  on  a  iiteasurement  of  the  amplitude  of  the  maximum  peak  in  the  S- 
wavetrain  (Richter,  ItWo;  Bath,  1981;  Alsaker  it  al.  19!H)).  I'lie  NORSAR  record¬ 
ings  of  Fig,  7,6,1  siiow.  that  the  position  of  the  maximum  peak  vary  stronglv 
from  one  region  to  anrjther  Fveiit.s  originating  within  iiie  I'ennoscandian  Sh'eld 
(event  1  and  event  "i)  will  usually  have  the  maximum  energy  a,s.soriated  with  the 
h,  phase  (group  velocity  ,9.,'>  km/s).  On  the  other  hand,  events  wit'i  propagation 
paths  crossing  the  .North  Sea  grabeti  structures  (event  !  i  or  events  originating  in 
oceanic  regions  (event  6)  will  have  the  their  iiiaximntn  energy  a.s.sociatetl  with  the 
.S„  arrival  (group  velocity  about  l.-h  ktii/s),  Iti  ailditi.in,  Kv;erna  and  Mykkeltvit 
(  1985)  have  shown  that  the  regiotis  iti  which  the  /.j  arrit.il  is  the  liominant  plia.se 
are  dependent  on  the  fre(|Uencies  cotisidered.  l.e..  the  S„  phase  betnini's  more 
domitiatit  a,s  the  fretpiencies  increa.se. 

'1  ht'  I’M  ini't hod  re<|nire  that  tin*  travt*l  limes  ol  tin'  coiisidt'red  phases  ,ire 
given  a  priori  for  all  target  areas.  For  optimntti  performatice.  one  phase  should 
be  associated  with  the  energy  maximum  of  the  wavi  irtiin.  f  rom  the  cttmplexitie.s 
described  above,  it  is  obvious  that  we  cannot  obtiiin  generii  formula.s  for  the 
travel-time  of  this  amplitude  peak  without  extensive  data  amilysis  and  regioiial 
mapping.  For  NORFSS  recordings,  we  have  from  the  study  i>f  Kvierna  and 
.Vl  vkkeltvf-it  (198.5)  an  idea  of  th.e  geiwraphii  al  regions  for  w  hich  .V.,  or  /. ,  is 
the  dominant  phase,  but  similar  information  is  currently  not  .nail.ibh'  for  other 
seismic  arrays  and  single  stations. 

From  several  years  of  experience  with  seismic  dat.i  from  local  and  region.il 
events,  we  know  that  the  energy  a.s.sociated  with  the  |>  phase  often  exhibits  its 
amplitude  maximum  several  seconds  after  the  initial  1’  onset.  'I'his  feature  is 
partly  illustrated  in  Fig.  7.6.1.  For  optimum  TM  computations,  it  is  ahso  ben- 
oficial  to  make  use  of  the  phases  for  which  the  travel-time  difference  is  as  largt- 
as  possible.  W'e  will  therefore  in  the  following  proceeil  with  the  first  arriving 
F-pha,se  {P„  or  and  the  pha.se  in  the  FM  analysis,  using  the  stamlani 
Fennoscandian  travel-time  tables  a.s  the  generic  formulas.  To  compensate  for  the 
uncertainties  in  the  positioning  of  the  maximum  amplitudes  of  the  wavetrain,  we 
will  introduce  so-called  time  toh'ranct's.  'Fhis  conct'pt  w'ill  be  outlined  in  one  of 
the  following  sections. 


Frequenry  hands 

To  ensure  optimum  performance  of  the  I'M  method,  we  introduce  bandpa.ss 
filtering  of  the  data  in  the  band  where  the  considered  phase  is  expected  to  have 
the  largest  SNR.  These  bands  are  however  dillicull  to  predict  as  large  variations 
occur  regarding  attenuation  properties  of  the  different  propagation  paths,  source 
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spectra  and  noise  conditions. 


In  the  context  of  monitoring  regions  witliin  local  and  regional  distai.ci'S,  the 
work  of  Sereno  (1991)  gives  an  excellent  picture  of  the  average  properties  of  re¬ 
gional  phase  attenuation,  .source  spectral  scaling  and  backgoiitui  noise  conditions, 
f  rom  an  ajisiimplion  on  the  event  magnitude  .\//_  aiul  the  epicentral  distance,  we 
could  use  his  results  to  predict  the  best  S.VIl  frequency  band  of  a  pha,se. 

We  will,  however,  in  this  preliminary  study  base  our  selectioti  of  filter  bands 
on  statistics  from  the  detection  proce.ssing  of  the  regional  arrays  NORFiSS  and 
.■MK’F.SS.  The  l.\S/I.MS  system  (Itache  il  ol.  1990)  is  used  for  routine  analysis 
of  data  from  these  arrays,  and  all  information  concerning  the  detected  seismic 
pha.ses  are  .stored  in  a  large  data  base.  The  statistics  on  the  dominant  fre(|uency, 
i.e.,  the  frequency  with  the  largest  SNR,  give  us  ati  idea  on  how  the  o|)timum  fre¬ 
quency  band  varies  as  a  function  of  epicent ral  distance.  I'he  statistics  cover  both 
NORKSS  and  .MIC'KSS  data  from  the  time  interval  1990/01/29  to  1991/04/2!t. 

The  //,  (I'j)  results  ar<'  givtui  in  'I'able  7,0. i.  and  show  large  variability,  es 
pecially  at  distances  below  r,00  km.  .\i  hirger  dislajices  the  fre(|iienry  hand  9 
to  ')  llz  cover  tlie  vast  majority  of  the  occurrem es.  To  retain  simplicity  in  this 
(ireliminary  study,  we  have  dnisen  to  use  the  9  to  llz  fre(|uen<y  band  for  the 
first  arriving  I’-phtise  at  .ill  distances.  For  larger  distance.-i  this  is  also  in  genera! 
agreement  with  i>re9.'.riions  ha.sed  on  the  results  of  Seri'iio  (  1991  ). 

The  /,,,  results  given  in  I  able  7. (>.2  also  show  large  variability  at  distances  be 
low  .700  km.  It  should  be  noted  that  tlie  doiniiiaiil  frequencies  for  I.,,  are  relative 
to  the  preceding  .S',,  (oda.  and  not  relative  to  liackgoiind  noise  conditions,  as  was 
the  situation  for  l'„.  We  want  optimum  performance  relative  to  hackgonnd  noise 
cot  ditions.  so  till'  /.,  st.ilistics  shoulil  be  interpreleil  with  some  caution.  On  the 
other  hand,  nnmerons  studies  of  /.,.  propagation  charactoristies  (a.o,.  liaumgardt , 
1990;  Soretio,  1991;  Kvti’rna  and  Mykkeitveit.  IftSfi)  eonlirm  the  “low-frequenry” 
nature  of  at  distances  above  .700  km.  .Mso  in  this  case  we  will  make  a  compro. 
miso  and  use  the  1  ..7  to  .9.7  Hz  freipiencv  band  for  /,.,  at  all  distances.  This  will 
give  close  to  o|ilimnm  performance  for  at  longer  distances,  which  is  considered 
the  most  important  for  the  overall  threshold  monitoring  capability. 


drtd  (U  finitions  and  itmf  toU  rdnri  i 

riireshold  monitoring  of  a  larger  geographii  al  r<-gion  implies  that  each  target 
jioint  have  to  represent  a  finite  surrounding  area.  If  we  divide  the  region  to  be 
monitored  into  a  grid,  as  sliown  in  Fig.  7.0.2.  tlie  area  siirroiinding  tin'  target 
point  i.s  given  by  a  rectan.gle  as  indicated  on  the  same  figure. 

Fhe  travel  time  of  the  considered  plia.se  is  given  by  l\.  where  denote  the 
distance  frotii  the  station  .1  to  the  target  point  .1/.  Let  be  the  minimum 
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travrl-lime  from  any  within  the  ret:tanglo,  e.g.,  M\,  and  lei  be  the 

fnaxinuJm  travel  time  from  any  point  within  the  rectanp^le.  e.g.,  \f2. 

If  the  density  of  the  grid  is  such  tliat  the  magnitude  calibration  factors  do 
n<)t  vary  .signifirantly  witliin  the  rectangle  surronndtug  each  grid  point  (see  Fig. 
we  use  the  following  procedure  for  monitoring: 

I. el  .S(/)  denote  th(‘  >ignal  level  observed  at  time  t.  Instead  of  meL'.-^Mriug  the 
signal  level  at  time  l\  as  predicted  from  the  p<»sition  of  tin*  target  point,  we 
intro(hice  time  tolerances  such  that 

~  max(.VtMi  (d) 

where  /  6  /  a>1-  I  herein*  the  ('siimated  signal  level  can  Ih*  sai<l  to  represent 

an  ii{)per  limit  for  any  M)urceh  witliin  the  rectangle.  The  lime  tolerames  can 
also  be  used  to  compensate  for  uncertainties  in  the  position  of  the  maximum 
amplitmle  of  the  svavetrain.  but  we  note  iliat  the  resolution  of  the  IM  method 
will  be  deteriorated  if  tlie  liiiH’  toh'rances  l)e((mies  too  Iarg‘’. 


S  l'A  h 

In  determining  the  (»{>(imum  SIX  window  length,  ne  need  to  t.ike  thre»’ 
facior.s  into  arronni; 

■  Average  Si  \  during  noise  <onilitions. 

•  \'aiial)ility  of  S I  \  during  tioise  e«>nditi'*ns. 

•  M<iximum  S I  .\  value  when  the  vinnal  ox«ijis. 

In  practici*.  it  is  de^irable  to  have  n  signal-io  none  ratio  a>  Kirue  poNsUd,'. 
measiire<{  relative  to  mulfiple.s  of  the  noiv-  standard  <hni;itjon.  Our  appro. uh 
toward  solving  tldn  problem  outline<i  in  the  following. 

In  this  initial  study,  we  havi^  cIm>s»ui  t»>  sample  ilie  ,lat.t  hv  1  ^etoMi  'lii.n 
1  erm  averages  ( .S' /  .t )  sampled  at  1  Ms  omi  interval'-  I  his  tien-ion  i-  b.tsed  .>i; 
coinpromise  between  <lata  resolution  and  mali.tgable  d.ita  volhines 

Intuitively,  an  iinl ant ane<*us  jdia^e  with  short  tlniatioii  te.g..  ‘  shoidtj  I.* 
represente<i  !)y  an  .S  /  \  average<l  of  a  short  time  vsindovs .  u  here. is  i  he  ,niip!ii  ude 
level  of  an  eineige-n!  jdiase  with  long  duration  I...f  slnadd  be  represented 

by  a  longer  time  window  ]  he  initial  data  sampling  I  1  se<  .  si  t  values i ,  .dlosK 
us  to  use  any  int«'ger  multiple  of  1  second  as  windt*v.  lengths  for  the  tonsideped 
jdiases. 

I.el  <lenole  ilie  average  of  I  he  logi.s/'  \)  untb'r  nois«‘  i  imdii  loiis.  .oib 

let  rrlA/)  b»'  tlie  associated  standard  ileviation.  A/  refers  to  a  partnail.ir  >  i A 
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window  length.  Let  y{At)  be  what  we  consider  the  “worst  <a.se’'  noise  situation 
given  by 

y(At)  =  I(At)  +  rrT(At)  (4) 

Let  .SlAt)  be  the  ina.ximum  of  log(5T.t)  for  the  signal.  We  introduce  the 
term  “noise  damping",  ^(At)  by  the  formula: 

tr(Af)  =  .S-(A/)-y(An  (5) 

The  “noise  damping"  is  then  a  measure  of  the  “effective"  signal-to-nois.’  ratio, 
i.e.,  how  much  the  signal  e.xceeds  the  “worst  case"  noi.se  situation.  The  optimum 
SI  A  window  length.  At,  is  the  arguif.ent  for  which  the  noise  damping  ;lAt) 
attains  its  ma.xiimim. 

To  assess  the  optimum  S'i A  window  lengths  for  /'„  and  L,  and  to  reveal  any 
distance  dependency,  we  computed  ina.ximum  signal  SI'A  tallies  with  different 
window  lengths  for  events  at  various  epiceiitral  distances. 

Lsing  the  z-component  of  the  center  instrument  of  .N'ORKSS,  .VKC'KSS  or 
I'lNF.S.A,  the  /'„  data  were  filtered  in  the  ;i-.')  11/  pa-ssband.  The  starting  jioint 
of  the  ST  A  windows  were  at  the  predicted  arrival  time  of  the  p. pha.se.  and  to 
accomodate  for  uncertainties  in  the  positioning  of  the  amplitude  maximum  of  the 
F’  wavotrain,  we  introduced  a  time  tolerance  of  ±  seconds.  Information  on  the 
!\  data  are  given  in  Table  .'1.  The  interpolated  curves  of  Fig.  7. 6. .3  give  S(Al)  for 
several  events  for  a  set  of  different  window  lenghts.  For  this  study,  the  absolute 
scale  of  S(At)  is  without  any  significance,  .so  for  display  purpo.ses,  an  offset  was 
added  to  each  of  the  curves.  .As  expected,  the  shortest  window  length  ( 1  second) 
gave  the  largest  S{At).  but  there  is  a  distinct  difference  in  the  slopes  for  events 
above  and  below  300  km  epiceiitral  distance.  We  will  therefore  in  the  following 
proceed  with  two  average  signal  curves,  one  for  all  events  within  300  ktn  of  the 
stations,  and  another  for  for  the  rest. 

The  noi.se  charai  leristics  for  the  3..'i  H/.  freiiuency  band  wins  obtained  from 
analysis  of  six  30  minute  noi.se  intervals.  Information  on  the  noi.se  intervals  are 
given  in  Table  7.0,1.  For  consistency  with  the  /’„  analysis,  a  time  tolerance  of  .3 
seconds  wa.s  used.  Values  of  .'IfAf)  for  all  noi.se  samples  are  given  in  Fig.  7.0.1. 
together  with  the  average  over  all  six  samples.  Similar  curves  for  of  Af )  are  given 
in  F  ig.  7,0..'). 

Now  turning  to  the  noise  damping  of  the  phase  for  events  within  300  km 
of  the  station.  Fig.  7.0.0  give  the  noise  damping  c(Af)  for  a  set  of  confidence 

levels  j  ■  irfAf)  (j  =  1.2 . .')),  and  show  that  for  any  choice  of  confidence  level, 

a  1  second  window  length  will  do  the  best.  For  events  more  distant  than  300 
km  from  the  station,  we  get  the  same  conclusion  a,s  inferred  from  the  results  of 
Fig.  7.0  7.  It  is  clearly  possible  that  a  shorter  time  window  than  I  second  might 
further  improve  the  /’„  pha.se.  but  we  have  not  so  far  investigated  this  possibility. 


I  he  dofinitioii  of  tho  “worst  rase  "  situation  is  soiii«'wliat  arbitrary,  but  in 
conjiirtion  with  tho  lota)  number  of  samples  per  day  (HO-tOO),  the  '.\a  level  is  a 
res(uiable  praetiral  ( i^mproinise.  This  nu'ans  that  «>f  the  data  will  be  below 

this  liniil.  We  also  see  that  for  all  confidence  levels  up  to  'tfr.  the  conclusion  on 
lh('  b<-st  window  iehg.th  for  will  remain  the  same. 

Sirfiilar  analysis  was  ron<lmter)  fr>r  the  /.^  pljasi-.  The  data  were  bandpass 
tillered  bet  we*u>  1  .'>  .md  d..',  H/.  and  the  cenO'r  point  of  i  he  si*;!^!  anal>  si'-  window 
w.is  set  at  the  e\pi‘cie<i  ampliitide  maximum  of  the  /.  .  phase  (i-*’..  tt  a  {^roup 
velocitv  of  3."»  kni/si.  lo  rtcco.iiodate  for  uncertainties  in  the  posit!i»iiin  ;  of  lh(' 
ampniiide  maximum,  we  used  a  lifue  tolerance  of  i.r>  seconds.  Details  on  the 
l.^  phases  are  e.iven  in  I'able  7. (>.3.  and  tlie  valti<‘s  of  -^1  AO  h)r  (‘vents  at  various 
disianres  arc*  siiown  in  l  ie;.  T.ti.s.  .Mso  in  tliis  case  evenis  above  and  below  300 
km  show  different  slopes,  and  we  will  in  the  following*  proce/>d  with  (he  averages 
for  these  two  populations. 

I’he  data  interval.s  of  I’able  7.0.3  were  also  to  assess  llie  noise  charac- 
t(‘ristics  of  the  l.r>-3.r>  ID  fretpiency  baml.  The  estimated  c\irv(»s  for  .tfAO  are 
given  in  Fig.  7.0.0,  and  the  corresponding  rr-valu<*s  are  given  in  Fig.  7. (>.10. 

The  noise  damping,  computed  from  “an  average"  signal  within  !i00  km 
opiccntral  di.stance  ami  from  “average"  tioi.se  conditions,  is  given  in  Fig.  7.0.1 1. 
When  considering  the  levels  '.\a  and  higher,  all  window  lengths  of  5  seconds  ox  less 
siHun  It)  do  almost  e<jually  well.  The  corresponding  t  virve.s  for  events  exceeding 
300  km  epirentral  distance  are  shown  in  Fig.  7.0.12.  They  indicate  that  an  ST  A 
window'  lengtfi  of  ftl'  s4*conds  will  be  close  (<>r>ptiniiim  for  all  coriridojire  levels  up 
to  O/T. 

Our  preliminary  as>esMm*nt  i.s  that  a  -3  second  window  length  should  b(‘  used 
for  i.j  phases  originating  from  eveiiis  within  300  km  epjcentral  distance,  whereas 
a  10  second  window  >hould  lx-  us«*d  for  4‘vmii.x  4‘x<«M‘ding  300  km. 

An  increa.M'  in  tli‘‘  tinu'  iob‘raiic<‘s  will  imr«‘as»'  the  values  r»f  .dfAt)*  whereas 
rr(  At)  will  derrxave.  Fig.  7.(*.13  illustrates  this  for  a  noise  sample  in  the  l.o-H.') 
11/  fr<M)uency  band  using  a  It)  second  ST  A  wimlow  huigth.  We  see  tliat  the  value 
of  .1(  A  / )  -f  3  A/ )  remain  almo.st  <  oust  ant  for  atiy  time  tt/lerance,  implying  that 
the  results  we  obtained  with  a  time  tolerance  of  i  .7  secoiuis.  also  siHun  to  he 
valiil  for  other  <  hoice;-  of  tiin«‘  lob*ranr<*s. 


SiftTvig  (if lays  aud  tjjfcfs  iff  ons-sifirxnti 

One  of  the  main  h'aliir^vs  of  seismic  arrays  is  the  ability  to  imj)rove  the  signal 
to  noise  ratio  (SN'H)  by  In^amfonning.  Instead  of  rompiitiiig  tin'  SlA'i^  from 
bandpa.ss  filtered  single  component  sensors,  we  steer  lieams  towards  each  target 
point,  filter  them  in  ilie  at>propriate  fr«'quen<  y  hands,  and  finally  compute  the 


ST  A  values.  In  this  way,  we  significantly  reduce  the  noise  levels  (for  uncorrelated 
noise,  by  a  factor  of  v/iV,  where  N  is  the  iiuiiiber  of  sensors).  Kvaerna  ( 1989)  have 
estimated  the  SNR  gain,  the  noise  suppression  and  the  signal  loss  for  P-phases, 
using  data  from  the  NORESS  array.  In  the  .'I-o  Hz  frequency  band,  appropriate 
for  it  was  found  that  an  SNR  gain  of  12  dH  cosild  he  achiev(‘d  with  optimum 
plane-wave  steering  delays.  It  was  also  found  that  even  though  the  array  was 
steered  with  optimum  st*'ering  parameters,  the  signal  am})!itiides  were  reduced 
by  the  beamforming,  due  to  lack  of  coherency. 


As  shown  in  Fig.  7.G.1-I,  the  steering  delay.s  (apparent  velocity  and  azinnith) 
appropriate  for  the  target  point,  will  not  lx*  optimum  for  the  rest  of  the  points 
witliin  the  surrounding  rectangle.  We  will  in  the  fidlowing  consider  the  “worst 
case"  situation,  and  account  for  the  maximum  signal  loss  for  any  points  within 
the  rectangle.  If  we  assutm*  that  lh<‘  expected  slownesses  of  all  pinnts  within 
the  rectangle  is  identical,  which  is  resonable  for  and  I.,,,  the  mis-steering  will 
primarily  he  caused  l)y  d<‘viating  azimuths,  as  shown  in  Fig.  T.G.l  t. 

Fig.  T.ti.lo  illustrate  ih<*  Ios.h  of  tlie  niaximnm  STA  a.s  a  fiinciion  of  mi.^- 
stwring,  for  NORf.SS  and  .AKCKSS  P  heams  filtered  between  .'i.O  and  5.1)  11/., 
Information  on  the  events  are  given  In  'I'ahle  7. (>.5.  The  app.irent  vi'lociry  of 
eacfi  phase  is  taken  from  broadband  f  k  analysis,  the  S  T  A  huigtii  is  one  second, 
and  tiie  time  tolc*ranc(‘  is  i-5  scM'oitds.  I  ho  mis  st<‘<Ming  is  ini ro<|ijced  as  azimulli 
deviations  normalized  relative  to  an  apj)arent  velocity  of  s.()  km/s.  Let  denote 
the  azirniitlt  deviation  relative  tc)  an  apparent  velocity  of  s.o  km/s  and  let  c, 
denote  the  apparent  velocity  of  the  incoming  wave.  If  is  t  he  azimuth  deviat  iott 
relative  to  Vj,.  we  gel  the  following  relation: 


=  2  arcsinf  —  sin  —  ) 

s.O  2 


(hi 


fig.  7.b.l.5shows  that  the  sigttal  loss  is  about  I  dH  fc>r  a  itormaliztMl  azimuth 
mis-stoering  of  20  dc'grees.  if  otir  gri<l  is  conslrticte<l  in  >uch  a  way  that  the 
maximum  allowed  azimuth  deviation  is  within  20  dc'grees  (sc'o  Fig.  7.0.1  1).  the' 
Pn  signal  loss  at  NORESS  and  AROESS  will  bp  within  I  dIE  For  arrays  with 
smaller  radius  (e.g.,  FINF.SA).  the  signal  loss  will  he  loss. 

We  have  not  .so  far  inv(‘stigatc‘d  the  signal  loss  due  to  azitmitli  mi, "-steering  of 
the  i,j  pha.ses.  I’he  appanuit  velocity  is  lowc'r  than  fc^r  which  indicate  highc’r 
signal  loss,  but  fh<*  lowc’r  frecpiency  tillc'r  iiscmI  for  I.,  (1.5-:i..5  Hz  vc'rsns  3. 0-5.0 
Hz)  works  in  the  opposite  dircTiion, 

Due  to  the  large  regional  variations  in  prcipagalion  charactc’ristics.  it  is  nsuallv 
diflicult  t(»  preclict  tlx*  aj>parent  velocitic's,  givem  the  co(>rdinates  of  tlic’  targc't 
point.  lahle  7.0.0  givc's  file  l•stiInaled  appanui!  velocity  of  the  first  arriving  P- 
pha.sc’  { Pg  or  /’„ )  as  a  function  of  epicentral  distance*.  I'he.sc*  statistics  are  taken 
from  lh»’  I.\S  data  base*,  and  contain  liofh  NORESJS  and  ARC  ESS  observations. 


Similar  statistics  on  tin*  I.,j  phase  arc  given  in  lahlo  iioth  tables  show  a 

large  scatter,  illustrating  the  <liirirulty  in  pnMlicting  the  apparent  velocity  given 
the  epicentral  (iistance.  Another  complicating  factor  is  the  dispersion  of  the 
wave  train,  itnplying  that  the  estimates  of  apparent  velocity  will  he  a  function  of 
l)oth  the  frequency  band  an<l  the  positioning  of  the  analysis  window. 

We  have  initially  not  attempted  to  do  any  systematic  regionalization  ot  the 
apparent  velocity  (d»ei vat it>ns.  In  the  mean  time,  we  use  an  appaieiil  velocity  of 
S.O  kni/s  when  forming  beams  sleer<*d  towards  t argei  points  more  di>oiiri  ti.  /. 
■JoO  km.  At  clo>'er  distanci's.  we  tise  (>.5  km/s.  l  or  b<Mms.  an  apparent  velocit  \ 
(jf  l..'t  km/.-'  i-N  assumed  for  target  point.',  at  all  distames.  I  hese  paraiiK'ters  are 
cuiTeutly  list  d  for  all  arrays  (NOKKSS.  AIK’KSS  and  I'lNI’SA). 

1  he  signal  loss  will  also  l)e  depiuolent  on  the  array  geometry,  but  this  ha.- 
-o  far  not  [xvn  studied  in  c(mu<M-(iou  witli  mis-st-'ering  <»f  tin*  beam.s.  A  natural 
iiext  -lep  will  be  to  <‘valual(‘  all  the  elfecls  of  beamfoniiitig  array  geotnet  ciew 
and  mis  ste<‘ring  iu  (he  context  cjf  thr^'shold  :muii(oring.  Hut  in  this  f>rehmi},;)ry 
St  udy,  1  h<‘  signal  loss  i.s  ac<  <>niited  for  by  ailding  a  const  ant  term  of  0.2  (  1  d  H  )  1 1 . 
t  111'  observed  [ogf  S I  ,  1 )  values  for  /'i  {  P-j).  and  O.-'l  { 0  d  B)  to  t  In’  logf  .S' 7  A  )  values 
for  I-., 


corrir///i?»  fii<  for>  (nit{  rnnafif'( 

We  are  now  in  t  lie  po.-^it ion  to  computi' t  lie  generic  relati<)n.s  for  i  he  inagnit tide 
rorreiiiou  factor.-,  a-  the  other  I'M  variables  have  lu'en  preliminary  assesse(l, 
.\Uaker  I  f  til  i  1000).  i-oli<M  !i‘(l  a  larg('  event  dat.i  base  w  hen  est  imatiiig  forniula,s 
for  a  .\//.  scale  in  Norway,  and  they  siibsetjur'iitly  (omjiuteil  neiwairk  averaged 
estimat.'-  for  all  evniis.  We  wj))  in  the  following  use  their  ilata  ba.e  arid 
inagiiitudes  a-  a  liasi-  for  c<»mputing  tlie  geneiir  relations  for  the  magn  Hide 
correct  ion  factors. 

The  data  base  coiji ains  observai ion-s  from  21  ditferent  s» at iotis  (see  Fig.  7.0. 1  (i 
most  of  whi<  h  with  differefit  Inst r»in»etit  re.-potrse  fur>'  mojis.  In  order  to  compare 
the  .S/'.\  values  at  the  n’spectivf'  stations.  w(‘  ne< d  to  jind  a  common  basis  for 
com fiari-nn.  .A.-  f  he  Individual  amplit  ude  response  fijiici ions  show  only  -mall  v.iri 
ations  within  tin*  relativi'ly  narrow  passbands  considered  for  ami  we  can 
in  a  ft  .tp()roxiftiaf  e  wav  t  ransfer  t  he  .S' 7  A  values  into  units  of  ii  in  or  ii  injs  simpl  v 
by  multiplying  by  the  dis|)laceinen(  or  vrlocilv  respons<‘  at  the  center  fre(piemy 
of  the  passbaiifl.  smli  tijat 

5;  )1  (7) 

wlirre  .S’r.t,..,  is  till'  iil)si'rviMl  SI  A  in  i|uanliini  units.  |,A,(|  is  tho  (lisplaci'itic'nl 
aiiipiit  iiilc  rnsponsr,  and  llic  rentrr  frcfpipnry  w  hern  w'|  and  ..'2  arc 

I  111'  low  and  IhrIi  i  nloffs  of  I  hr  passhand.  A  similar  lypr  of  riinalion  r;  a  hr  nsrd 


if  we  instead  convert  the  ST  A  "alues  to  ground  velocity. 

In  accordance  with  earlier  regression  analysis  of  magnitude  relations  {.Msaker 
"'e  the  following  parameterizition: 

A/i  =  log  .S'7’.-l,  +  Cl  +  C2  •  log  A,  +  f '.'t  •  A,  (;=1, {S) 

where  .V  is  the  nnniher  of  observations.  M,  is  the  network  magnitude  of  the  event, 
STA,  is  the  instrument  corrected  and  A,  is  the  epicentral  distance. 

The  data  ba.se  of  Alsaker  el  al  (1990)  contains  711  observations  distributed 
among  lO-A  events  (.see  Kig.  7.6.17).  To  ensure  good  SNl?  in  the  r„  and 
frequency  bands,  all  data  were  visually  inspected.  .-After  rejecting  data  with 
insufficient  S.N’Ii  or  with  other  data  quality  problems.  <d)Scrvatioris  remained 
for  analysis  and  .72.'<  for  A.,.  The  STA  value.s  wer/-  rom|)uted  using  the  recipes 
outlined  in  the  preceding  sections,  and  the  results  from  the  regres,sion  analy.ses  are 
given  in  fable  7.6,H.  Fislimates  of  the  standard  deviation  are  also  given,  and  .show 
a  a  value  of  0. 19  for  1.^.  The  /'„  data  show  a  much  larger  scatter,  and  we  obtained 
a  f7  value  of  0.36.  Compared  to  site  specific  moiiiloring.  these  tr  estimates  are 
significanlly  higher,  as  the  typical  rr  values  for  site  speciPc  monitoring  are  less 
than  0.2  for  P,,,  and  le.ss  than  0.1  for  If  dilTerenl  filters,  travel-time  modefs  or 
other  parameters  were  to  be  used  in  the  TM  analysis,  new  magnitude  rorrection 
factors  would  have  to  be  obtained  from  reanalysis  of  the  calibration  events,  using 
the  new  reripies. 

•As  the  T.\I  method  computes  upper  magnitude  limits  from  a  cumulative  dis¬ 
tribution  with  a  given  mean  and  .standard  deviation,  we  have  the  ojition  of  bal¬ 
ancing  the  term  <"'1  against  the  standard  deviation  n.  This  implies  that  we  can 
reduce  n  if  ("1  is  increa.sed.  Our  philosophy  ftehind  the  I'M  computations  ha.s 
been  to  make  conservative  estimates  of  the  upper  tnagnitude  limits,  in  order  not 
to  overestimate  the  capabilities.  In  this  way,  wc  can  add  a  constant  term  to  Cl 
or  increase  rt  if  some  of  the  attenuation  relations  or  other  underlying  parameler 
estimates  of  the  TM  method  are  considered  particularly  nncortain. 


/ti.srn.s.sioTi 

fhe  results  presented  in  this  study  give  us  a  tiieans  of  testing  the  concept 
of  threshold  monitoring  applierl  to  large  geographical  regions.  It  enalde.s  us  to 
extend  the  original  “sile-specifie”  threshold  monitoring  to  what  we  might  call 
"regional  threshold  moiiiloring”.  I’sing  these  initial  generic  relations,  Ringdal 
and  Kvmrna  ( 1991.  this  issue)  have  already  shown  how  colour  compuler  displays 
can  he  applied  to  interpret  the  results  from  TM  analysis,  fhey  also  indicate 
new  applications  of  the  regional  Ihre.shold  monitoring  coni  ejit  which  should  lie 
investigated  in  parallel  with  improvements  of  the  generic  relations. 

I  he  data  base  used  for  obtaining  the  magnitude  calibration  factors  consists  of 


events  fron)  Fennoscaiuiia  and  adjacent  areas,  making  the  results  representative 
for  this  kind  of  j»**o!ogi('al  environm<‘nt .  If  we  want  U»  extend  the  1  M  analysis 
to  other  types  of  gmlo^ira)  regions,  exhibiting  different  wave  propagations  char 
act  eristics,  new  generic  relations  hav«*  to  be  found.  Another  uncertain  factf^r, 
/'op/erniijc  Hx*  *  [jiagfiit ude  caJjbrat-'oji  is  efferl  c'f  tins 

particular  data  base  f(»r  regression,  as  the  same  dat.i  ba>e  wa.^  u^ed  for  obtaining 
t  h(‘  A//,  scale  for  .N’orvvay  (Alsaker  tt  a/, 

1  lie  effect  of  signal  loss  due  to  iiii>.-.steeriiig  ol  the  an.iy>.  'hould  be  more 
tlioroughly  investigated.  I'lte  signal  loss  is  a  fniii  tion  *»f  several  variables,  among 
otfuTs:  fiha.se  type,  signal  coh<'renrv.  frequemy.  degree*  <if  mi.s  steering  and  array 
geometry.  This  also  imjilies  that  when  m-w  arrays,  with  dill'erent  array  geome¬ 
tries,  are  introdma.'d  in  the  i\t  < om()ufaf ion.",  new  nnaJels  for  .signal  h^.s.s  have  t(» 
be  a.ssessed. 

We  are  also  investigating  tin*  }>osMl>ilify  of  using  'ev«’ral  fiber  bands  when 
ref)resenting  the  amplitude  i(*vel  of  a  phase.  The  <  urrent  n-odej  of  a  fixed  fre 
‘pjenry  band  for  /*„  and  /.^  is  clearly  not  ofitimal.  Hut  iii  order  to  make  such 
irnprovenients.  new  generic  relations  have  to  fie  obtained  for  a  set  of  different 
filter  bands. 

Regionali/alion  of  the  travel-time  models  for  the  maximum  amplitude  peaks 
in  the  wavelrain  will  opi.ittii/e  the  I’M  comfiul  at  Ions.  I'he  data  base  of  .\lsakor 
(/  n/ ( 1990)  contain.^  several  rec  ordings  at  NORiiS.'^  and  .MUTSS  wliich  can  be 
u.sed  to  legionali/e  tin*  trave)-iiim‘  models  at  tlu'se  two  stations,  iiut  lor  the 
other  stations  currcmt Iv  providing  digital  data  to  .\'()1<S.\  U  f  I'l N  KS A.  ( iKU  I .SS. 
Ksiaz  and  .Sfary  Kolwarkf,  a  new  event  <lata  f/a,se  will  have  to  }>»'  fsillecied.  If 
iiulefiendetit  network  averaged  magnitude's  can  Ive  provided  for  lh«‘s<‘  c'vcmts.  the 
generic  rc'lations  for  magnitude*  calibration  can  also  fie*  tmfirove'd. 

In  conrlusion,  the  ke*y  for  further  improvements  of  th*'  generic  re'latienis  for 
re-gjona.l  uiesfiold  monitoring  is  ea.sy  Hcce*ss  to  a  large  event  data  base  includ¬ 
ing  recordings  at  all  relevant  .staf ie>n.s.  Network  locations  an<)  network  averaged 
magnitude's  should  be*  av<diable  for  all  e*venls.  With  this  at  hand,  we  have  the 
pos.sibiiity  to  investigaf*-  regional  be'havionr  atid  the»  <‘fr<*ci  of  elifb'rent  fiarame- 
t»'r  settings,  in  oreier  to  furtfM*r  imf>rove  the*  fx'rtormaue e  r*'gional  thresimld 
monitoring. 


Tormod  Kvaerna 
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Origin  time 

Lat. 

Long. 

Distance 

Ml 

Station 

Pn  data 

Lq  data 

J990-101;ll. 51.5.5.'! 

62.8 

170  8 

FINKSA 

ye« 

ye.H 

1990-101:13.46.07.0 

60.9 

29.3 

183.8 

FINFSA 

yes 

yes 

1989-167:11.23.26.0 

69.-4 

30.6 

200.0 

3.0 

ARCESS 

yes 

yes 

1989-076:1 1.48. 53.0 

69.4 

30.6 

200.0 

2.9 

ARCESS 

yes 

ves 

1990-101:10.21.21.0 

59.5 

25  0 

224.8 

FINESA 

yes 

ves 

1989-033:18.28.55.0 

67.1 

20.6 

338.5 

2.5 

ARCESS 

yes 

ves 

1989-0.59:18.3«>.4.5.0 

67.1 

20.6 

3.38.5 

2.5 

ArK;(:ss 

yes 

ves 

1989-105:08. .50.53.0 

68.1 

.332 

348.9 

2.7 

ARCESS 

ves 

yes 

1989-133:08:!8.49.0 

68.1 

33.2 

.348.9 

2.7 

ARCESS 

yes 

ves 

1988-258:08.59.58.0 

64.7 

:i0-7 

.584.0 

2.9 

ARCESS 

yes 

ves 

1988-141:09,54.24.0 

.59.5 

2.5.0 

760.0 

2-7 

NORESS 

no 

yes 

1989-0.51:13.19.57.0 

.59.5 

25.0 

760.0 

2-5 

NORESS 

no 

yes 

1989-108:13.41.15.0 

.59.5 

26.5 

841  0 

2.8 

NORESS 

no 

yes 

1988-075:11.52.22.0 

61-9 

30.6 

882  2 

2.8 

ARCESS 

yes 

ves 

1990-10.3:10.18.55.0 

59.2 

28.1 

9.37.0 

3.1 

NORESS 

ves 

no 

I989-00.5;l0.tXJ.07.0 

61.9 

-30.6 

1024.3 

2..5 

NORICSS 

no 

ves 

1988-258:08.59  -58.0 

64-7 

.30.7 

1069.4 

2.9 

NORESS 

no 

ves 

1990-103:10,28,41.0 

64.6 

31.2 

109.3.7 

3.0 

NORESS 

ves 

no 

1989-090:12.16.17.0 

.59.5 

26.5 

1119.8 

3.0 

ARCESS 

no 

yes 

101X)-103:08..37.0R.O 

67.6 

33.5 

1302.7 

2.8 

NORESS 

yes 

no 

1989-167:1 1.23.26.0 

69.4 

:K).6 

1307.3 

.3.0 

NORESS 

no 

ves 

1989-108:08.59,23.0 

68.1 

33.2 

1314.3 

2.9 

NORESS 

n4> 

yes 

Table  7.7.3.  Inforination  on  Iheovcnts  used  for  roniputation  of  maxinmtn  signal 
ainplitndps.  (ionott'd  .9(A<). 
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Start  tiit»<* 

Stat  ion 

1990-()9<KiZ.S0.00 

NOKrSsS  1 

l99O-0’^v2:}.00.UO 

AHCKSS 

Noar.ss 

AliCKSS  1 

NOKKSS 

09f»;0;K(KJ  00 

AHCKSS 

Table  7.7.4.  Stan  limos  of  noi.'io  intervals  used  for  assessing  average  noise 
eliarartorislii  s.  The  letigtli  of  all  intervals  were  .'iO  ininnte.s 
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■'ii  1 
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i^>*)irjo;(n  12.^4. 1 
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>5.0 

aa.ii  ll 

1S9M20;0'J  19  -'T  O 

*•7.0 

:h..9 

AfO'T:SS 

1001  U'O  fr.i  2(1  .'I'v.t 

120  4 

7  S 

I2S  )  ; 

lWM20.ll.VJ.2.’1.0 

:V2.0 

ATK’KSS 

loOMiO  12-00  4'.  H 

Kri,,; 

X.7 

sri 

lOlM  120:12. .M.4r>.0 

♦iO  1 

HI  0 

AHCKSS 

l<*l*M2(l,12.  l-.i  7 

•M.  1 

7  7 

28'<.4  !! 

Table  7.7.5.  List  <'venls  userl  for  the  preliminary  assessmi’nl  of  signal  loss  due 
to  mi.s-.steering  of  the  P-beaiiis.  Ihi-  event  locations  are  the  antotnatir  network 
solntions  from  the  eeirerali/.ed  heainformine  nieihod,  see  Ringdal  amt  Kvterna 

( iftsrt). 


Table  7.7.6,  I  iiis  ial)l<‘  an  ova-rvicw  of  tin'  ^'stiinatiMl  aj>{>arfMil  v<’I(jrity  f)f 
{ h('  first  arrivih*^  P-phasp  (  /\,  or  l\j).  I'arh  olrnuMii  of  i liis  table.  t he  nninlxT 
of  observat  ions  r)f  t  Ik*  apparerP  v«‘lo<  ily  for  a  ftiven  appar‘*iil  V(‘lociiy  and  anre 
range.  The  data  are  taken  from  routine  d<‘Uation  proia'ssing  of  the  I.\S  sysitnn. 
and  tfie  ^taf ( over  both  NORKSS  ami  .\K(TiSS  data  from  the  time  interval 
l<)00/0l/2:t  to  Ail  events  were  below  .U/,  .'TO. 


Table  7.7.7.  Same  as  Table  7.7.t).  but  for  the  I.j  [)hase. 
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Cl  Ci 
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nobs  { 
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i.uafi 
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OIZ 

•  U.'J.VC  U.Ki'IJ 

i.aisHK-o.i 

0.1U2 

r,2H 

Tabic  7.7.8.  Results  froin  n‘Rressi(Mi  a||aly^is  nf  tfje  <lala  us<Mi  by  Alsaker  ft 
(il  (Ib'iO).  riie  reores^ioIl  eoelfirien!  s  ;ui<l  I  lie  (J  value>  for  f\,  (/’j)  and  f,.,  were 
()l>taine(|  froin  proressing  tlie  data  with  the  'CM  recipie^  outlined  in  the  precedint; 
>ert  ioti^. 
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Fig.  7.7.1.  Illustration  of  variation  of  rrlativo  iinportanco  of  thr  phasos  and 
l.g.  riio  standard  group  voloritios  of  I.T*  and  kiti/s,  (oiunionlv  assignofl  to 
-St,  and  Lg,  respoctivoly,  arc  markcil  by  da-shed  linos,  riif^  upper  throe  traces 
co’.er  tfio  distance  interval  d.^0--5-'»0  kni,  while  the  throe  lower  traces  correspond 
to  opicentral  distances  in  the  range  1225-i;{2«  km.  The  location  of  the  NORSAR 
array  is  denoted  by  a  ring  on  the  map,  and  tlie  traces  are  from  the  NORSAR 
seismometer  02B01.  1  he  data  are  bandpass  filtered  1  to  S  H/.  The  reduction 
velocity  is  8.0  km/sec. 
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Fig.  7,7.2.  This  illnsl  l)ic  nort'ssily  t)f  iisinj'  tinir  toloranrrs.  TIh'  plus 

siRfis  indicat*’  tarj;«’i  points,  hihI  a  rr<  taiigl<»  surrounding  oik*  of  tho  targ»’t  points 
(A/)  is  als(j  givnn.  1  h<’  point  within  llin  n’riangU'  with  thn  ininimuin  travrl-tiinr 
is  dcnotrd  A/j,  whi'rnas  itin  point  with  tho  inaxiinum  travel  ti mo  is  donotod  A/j 


Relative 


statistics 


STA  leneth  (sec) 

Fig.  7.7.3.  I  hr  as  I  orisks  of  this  sIk)W  «»l)Sf*r\al  ions  <»f  maxinnuii  lo£;{  .S'  /  .1 1 

(dmotod  .S'(^0  )  for  forasrt  of  .S /M  hs.  l  lio  ol>.s»'rva1i(»ns  rorrospondiiig 
to  til#'  sain#'  {)has#'  ar<*  intrrpolat#*#!  bv  <lHsh#'d  or  <lott#M|  lin<*s.  and  thr  #'])ir«'nt ral 
dislaiKO  of  oarli  «'V#'iil  is  in<ii<*at#'d.  Infornialion  on  tlio  cv#>nts  aio  givon  in  'Fabb' 
T.T.d.  For  display  [)ur|M)s#'s  an  ojfs#*!  was  add#'d  to  oarh  of  tho  curvrs,  as  thr 
absobitn  seal#’  is  without  aiiv  signiliranro.  Not#*  th#'  difr#’r«'n<<'  in  thr  sUipps  for 
f’Vcnts  abovp  and  bolow  dOO  kin  f’pir#'iilra|  distance. 

H9 


Fig.  7.7.4.  F'.jfig  from  Um'  r*'ji>*'r  n?  AO  of  l)ot}i  ihr  NOKl 

.iml  t)ii‘  AlU  rSS  .irray.  llii>  litiutf'  >-lio\v>  tin-  a\ 'T.i*;*' ol  joii(  >  /  A)  uml'-r  in)i>*’ 
‘•oii'liiions  .OAM  )  for  a  "''I  I'f  >' /  . \  fis.  1  h<- d.ila  w<‘r<' lilNT‘’(j  in 

» li»->  pa^''l)ajnl  A  't  11/  and  a  f  iiii*-  >ol«TaiM  of  j  ">  \va>  ux'd .  I n formal  hum >ti 

t  luMiai a  jni‘‘rv<il>  u.ivon  in  Tal>l«'  7.7.1.  I  In-  a\ar.i‘;<M>t  !  Im  six  noisi-  .d)srrvat  iori'*. 

for  snl>s«Mjn»ml  analysis  ot  ftojsr  clampiiiL’,.  i.s  indii  ai'^i  iiv  filU'd  s^juar*  <  ami 
a  solid  lino. 


Fig.  7.7.5.  lids  litrnr'*  j’ivos  1  in*  standar<l  do\i,tiion  of  tho  nois#'  ofis{'rv;ttions 
a naly /od  in  l  it;.  7.(1.  I  ho  avora^o  of  t ho  st aiolard  doviaiion  (  n r  \ os .  nsod  for 
sul)•^onl|^-nl  anal>sis  of  nois«>  lampinji,.  is  indi<  alo<l  l»v  fillod  stpiaros  ami  a  solid 
lino. 
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Fig.  7.7.6.  Tho  ri'lalivo  iidiso  daiiipiiiR  given  in  this  figure  is  rompnted  from  av¬ 
erage  /’„  signal  hetiavior  for  events  within  .'iOOkin  eplrentral  distance  (taken  from 
Fig.  7.7,:!).  averagi’  noise  .-omlitions  (taken  from  Fig.  7.7.4 )  and  average  values  of 
noise  standard  deviation  (taken  from  Fig.  7. 7. .a).  Tli.’  relative  noise  damping  for 
a  set  of  confidence  levels  is  shown  and  the  3(7  level  use<!  to  characterize  a  "worst 
ca.se"  situation  is  given  by  fillet!  .stpiares  and  a  solid  line. 
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Fig.  7.7.7.  Same  a  Fig.  7.7.r).  bnl  representing  events  with  epicentra)  distanres 
exceeding  300  kin. 


.stat  istics  ( i  .5 


3  f)  Wv.) 


2  -1  ' 

-  ^ 

m  -  - 

AVERAGE 

K 

ARCt>;h' 

< 

^.3 

*1  ' . 

o 

NORPS."' 

c£ 

C 

2  2 

L  vit 

o 

^  ^ 

^  « 

-  -  ^  -  -# 

2  1 

-  0 

B 

—a 

C 

2,0 

'  G 

0 

-  - 

> 

o 

o 

o 

o 

1.9 

o 

o 

o 

1,8 

f 

, 

) 

5 

!0 

If) 

20 

STA  length  (sec) 


Fig.  7.7.9.  Same  a  Kig.  7.7.1.  hut  tlio  noise  inl<Tvals  were  analyzed  in  tlie 
Hz  filter  hand. 
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Fig.  7.7. iO.  Ihis  ligiire  gives  the  slan<lar<l  deviation  of  the  noise  observations 
analyzed  in  f'ig.  7.7.0.  ’Die  average  of  the  standard  deviation  curves,  used  for 
subsequent  analysis  of  noise  (lamping  is  indicated  by  filled  squares  and  a  solid 
line. 
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Fig.  7.7.11.  rin’  r‘'!;uivc  !icmvi.  iliunpiuii  <;iv*Mi  iti  tliis  ligun’  is  romputod  lr<Mn 
avt’rau*'  /.  ^  siuii.d  l.«‘liav i')ur  foi  within  .{0(1  kin  ral  disfann'  (  fak^'r; 

Innu  I  iti.  7.T.S1.  ava-rajif-  nois*-  • « nulii ioi!>>  (t.ak^'n  tiom  1  ii*.  7. 7. in  and  avt-rajif' 
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Fig.  7.7.13.  Th  is  figure  gives  the  meat,  of  log(.S'7’.\)  togei  lier  with  the  mean-*  firr 
level  for  a  set  of  time  tolerances.  The  first  noi.se  segment  of  Table  7.7.1  was 
bandpass  filtered  between  1..5  and  3.7)  Hz.  and  the  estimates  were  obtained  using 
an  ST  A  length  of  10  seconds. 
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7.8  Continuous  threshold  monitoring  using  “regional  threshold 
displays” 

Introduction 

Continuous  threshold  monitoring  (Ringdal  and  Kva-rna,  1989)  is  a  method 
of  monitoring  seismic  amplitude  levels  for  the  purpose  of  assessing  the  largest 
size  of  events  in  a  given  target  region  that  might  go  undetected  by  a  monitor¬ 
ing  network.  The  method  has  recently  been  implemented  v.ithin  the  Intelligent 
Monitoring  System  (IMS)  (Bache  et  at,  1990).  In  previous  Semiannual  Technical 
Summaries,  as  well  as  in  the  present  issue,  several  examples  of  application  have 
been  presented.  In  particular,  Kvaernaand  Ringdal  (1990)  conducted  a  one-week 
monitoring  experiment  of  the  Novaya  Zemlya  test  site  using  the  Fennoscandian  re¬ 
gional  array  network,  and  concluded  that  continuous  threshold  monitoring  down 
to  event  size  as  low  as  mj,  =  2.5  appeared  feasible  for  this  site. 

Regional  threshold  monitoring 

In  the  current  IMS  implementation  of  the  TM  technique,  a  limited  number 
of  specific  target  sites  are  monitored.  These  sites  include  several  mines  in  Scan¬ 
dinavia  and  Western  Russia,  along  with  the  Novaya  Zemlya  and  Semipalatinsk 
nuclear  test  sites.  For  each  of  these  sites,  a  number  of  calibration  events  are 
available,  and  thus  it  has  been  possible  to  fine  tune  the  parameters  in  order  to 
obtain  close  to  optimum  monitoring  performance. 

“Regional  threshold  monitoring”  is  defined  as  an  extension  of  the  original 
“site-specific”  threshold  monitoring  concept.  It  entails  using  the  same  basic  prin¬ 
ciples  to  obtain  wide  geographical  coverage,  including  coverage  of  regions  for 
which  no  calibration  events  are  available.  The  key  to  achieving  this  is  to  de¬ 
velop  “generic”  relations  for  attenuation  and  magnitude  corrections  of  seismic 
phases  of  interest,  and  to  deploy  a  sufficient  number  of  beams  to  ensure  adequate 
geographical  coverage. 

Kvaerna  (1991,  this  issue)  has  developed  initial  such  generic  relations  for  the 
Pn  and  Lg  phases  of  NORESS,  ARCESS  and  FINESA.  His  relations  are  appli¬ 
cable  to  Northern  Europe  and  adjacent  regions,  and  are  based  on  a  systematic 
analysis  of  several  hundred  phase  observations  of  regional  events  in  various  geo¬ 
graphical  areas.  Kvaerna’s  results  form  the  ba.sis  for  the  study  presented  in  this 
paper. 

Threshold  maps 

The  regional  threshold  monitoring  approach  lends  itself  naturally  to  displays 
in  the  form  of  contoured  geographical  maps.  By  using  a  spatial  grid  covering 
the  area  of  interest,  interpolation  can  be  applied  to  get  a  visual  representation  of 
threshold  variations  over  an  extended  geographical  region,  and  examples  will  be 
given  later. 
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These  contour  maps  are  in  many  ways  similar  to  the  standard  network  capa¬ 
bility  maps  traditionally  used  in  seismic  monitoring  studies  ( N'etworth,  Snap/1), 
etc  ).  However,  there  are  some  fundamental  dilferences; 

-  Standard  capability  maps  use  as  a  ba.sis  stali.slical  models  of  signal  and 
noi.se  rharacferi.stir.s;  in  particular  a  signal  variance  and  a  noise  variance  is 
assumed  to  compensate  for  statistical  fluctuations.  In  contrast ,  the  regional 
I'M  maps  give  ‘'‘snapshots”  of  the  capability  as  actually  observed  at  a  given 
point  in  time. 

With  standard  maps,  no  allowance  is  made  for  iiiMisiiai  conditions,  sm  h 
as,  e.g.,  the  occurrence  of  a  large  earlhquak<'  or  aii  aftershock  secjuence 
which  may  cau.se  the  network  capability  to  deterlorati'  for  hours.  With  the 
IM  aitproach,  the  tictual  variation  in  iletection  capablily  is  iiiimeiiiately 
tipparent . 

•  Sttindard  ca|)abibty  majts  require  .i.s.Mim|)tioiis,  e.g..  with  regard  to  "SN’K 
threshohi  required  for  detection"  and  ‘■minimum  number  of  stations  re 
(jiiired  to  locate".  The  I'M  itiaps  nniuire  no  such  assuitiplions  sime  they 
are  not  tir'd  to  "detecting  anil  locating"  seismic  events,  but  rather  describe 
directly  the  observeil  '‘seismic  field"  at  any  point  in  lime. 


W'e  will  briefly  comment  further  oti  the  last  itetii  itientioiied  above-.  I'lie  re. 
(|uiremeiii  of  mult isttil ion  detectioti  with  the  standard  tiiethod  will  sometimes 
resiill  iti  unrealisiictdly  high  thresholds,  e.g.,  jn  areas  near  a  station  of  the  mon 
iloring  network.  The  imiliistation  reiptiri'inent  ahso  imjdies  that  the  method  is 
not  able  to  adequately  represent  the  possibility  of  particularly  favoralde  source 
station  paths.  ca.so  in  point  is  the  out  standing  capahiliiy  of  l  lie  .NOliKSS  array 
in  detecting  exiilosions  ;it  .Shagan  River.  Thus,  if  .NORK.SS  has  no  detection,  it  is 
highly  unlikely  that  any  explosion  at  llml  site  of  lUf,  >  :i  has  occurred,  whereas  a 
capability  map  based  on  l-slalion  detection  reipiireiueni  may  well  show  a  thresh¬ 
old  an  order  of  maguilude  liigher. 

The  threshold  moiiitoriiig approach  will  avoid  (he.se  inconsistencies.  I  hiis.  un 
der  noritial  noise  conditions,  the  thresholds  will  be  very  low  witltin  a  few  hundred 
km  of  each  network  station.  I'lirlhermore.  .since  the  I  .M  thresholds  are  dominated 
by  the  "best"  station  of  t  he  network,  part icularly  favorable  source/ receiver  pat  hs 
may  be  accommodated,  althoiigli  this  would  require  ,i  combitiatiori  of  regional 
and  slte-s()eritic  riiouiloring. 

Diaptny  I  iai}ipli  s 

I'siiig  the  generic  relations  developed  by  Kva-rna  (1111)1),  we  romputed  a 
tlireshold  monitoring  grid  of  20  x  '20  geograpliical  aiming  points  for  a  40  mimite 
time  interval.  Data  from  the  three-  arrays  NORT.SS,  .-MICKSS  and  TI.N1.S.\ 
were-  used,  ('onionring  ma|i.s  wi-ri-  develop<-il  iiy  interpolation  in  this  grid,  and 


displayed  in  the  form  of  color  maps  where  the  color  scale  is  tied  to  the  actual 
threshold. 

Figs.  7.8.1  and  7.8.2  show  two  repre.sentative  examples  of  output  from  this 
procedure. 

Fig.  7.8.1  shows  the  “absolute”  TM  threshold  levels  (with  m;,  units  indicated 
on  the  color  template)  at  a  specific  time  during  a  typically  “quiet”  period  (i.e.,  no 
seismic  event  occurring).  We  note  that  the  areas  immediately  surrounding  each 
array  (deep  or  light  blue)  show  the  lowest  thresholds  (below  mj  =  0.5),  whereas 
most  of  the  remaining  area  at  regional  distances  has  a  green  color,  indicating 
thresholds  in  the  range  m(,  =  0.5  1.5.  The  yellow  color  seen  further  away  from 
the  network  stations  indicates  thresholds  of  1.5  to  2.5. 

Fig.  7.8.2  shows  a  typical  map  at  a  time  corresponding  to  a  mining  explosion 
(magnitude  2.2)  at  the  Apatity  mine  in  the  Kola  Peninsula.  In  contrast  to  Fig. 
7.8.1,  we  have  here  chosen  to  display  relative  thresholds  (i.e.,  thresholds  relative 
to  the  average  thresholds  during  noise  conditions  at  each  geographical  point). 
This  is  done  to  emphasize  more  clearly  the  effects  of  the  .seismic  event  in  causing 
threshold  increases  outside  'he  source  area.  We  note  that,  naturally,  the  area 
surrounding  the  mining  site  has  the  highest  relative  threshold  (red),  whereas  the 
“side  lobe”  effect  causes  significant  threshold  increase  ahso  in  other  regions,  some 
of  which  quite  far  apart  from  the  mine. 

The  computer  displays  shown  in  Figs.  7.8.1  and  7.8.2  also  include  fields  for 
displaying  threshold  traces  and  selecting  various  plotting  options.  At  the  present 
time,  however,  these  features  have  not  been  operationally  implemented. 

Perspectives 

We  consider  that  the  regional  approach  to  threshold  monitoring  would  imply 
a  significant  enhancement  of  practical  monitoring  of  underground  nuclear  explo¬ 
sions.  In  particular,  a  graphics  display  system  could  be  developed  to  provide 
the  analyst  with  very  useful  interactive  tools.  Among  features  that  might  be 
desirable  are; 


“Snapshots”  of  regional  threshold  maps  taken  at  times  when  a  peak  occurs 
on  a  threshold  monitoring  trace.  For  example,  if  a  peak  is  obcC.ved  on 
the  threshold  trace  used  to  monitor  Novaya  Zemlya,  such  a  snapshot  could 
immediately  reveal  that  this  peak  might,  e.g.,  be  a  side  lobe  effect  from  a 
remote  earthquake. 

Threshold  displays  taken  during  the  coda  of  very  large  earthquakes,  indi¬ 
cating  the  resulting  effects  on  detectability  in  various  regions. 

“Cumulative”  displays  showing  the  largest  possible  events  that  might  have 
occurred  during  a  given  time  period  (e.g.,  24  hours). 
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-  ('ombiii;itioiis  i)f  tlirosliold  displays  and  conventioiiaJ  ('picciiter  maps  (if  d('- 
tcctcd  ('vents. 


.\it  e.xl ri'inely  interesting  ap.pliration  woidd  ho  a  real  lime  "vidod’'  display  of 
luiw  the  tliri'shol  l  situation  flurtnatos  with  time.  When  a  seismic  ovt'iit  occurs,  a 
real  time  display  of  this  typo  would  illuslralo  how  the  threshold  first  increases  at 
"side  lohe”  loca  ions,  with  suh.sequenl  focusing  upon  the  actual  epicentral  area. 
.Snell  a  vid(Hi  option  could  of  course  just  as  easily  he  implemented  for  off-line 
(retroactive)  display  of  lime  periods  of  interest. 

In  order  to  'take  efh'ctive  use  of  the  regional  threshold  monitoring  approach 
and  the  associated  display  options,  a  workstation  with  powerful  computational 
and  graphii  al  <  ijiahilities  will  he  required,  and  we  are  currently  evaluating  possi 
hili ties  in  this  r-  gard.  We  are  also  continuingour  re.search  aimed  at  integrating  the 
“regional"  and  *site-specifir"  threshold  monitoring  methods,  whirh  we  ronsidi'r 
to  have  a  rom  ined  potential  of  heroming  a  basic  tool  in  practical  monitoring 
applicat  ions. 


F.  Ringdal 
T.  Kvaerna 
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by  Kolf  M.  .-len  of  N'OKS.Mf,  using  the  “NOGR.A"  graphics  software  system. 
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Fig.  7.8.1.  Example  of  regional  threshold  display  of  “absolute"  threshold  levels,  at  a  typi¬ 
cal  "quiet"  period.  See  text  for  detailed  explanation. 
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Fig.  7.8.2.  Example  of  regional  threshold  display  of  "relative"  threshold  levels  at  a  time 
when  a  mining  explosion  iKetirrcd  in  the  Kola  Peninsula.  .Sec  text  for  detailed  inlormation. 
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